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FOREWORD 


This  report  presents  the  final  results  of  one  of  the  46  projects  comprising  the  military -effect 
program  of  Operation  Plumbbob,  which  include  24  test  detonations  at  the  Nevada  Test  Site  in 
1957. 

For  overall  Plumbbob  military-effects  information,  the  reader  is  referred  to  the  “Sum¬ 
mary  Report  of  the  Director,  DOD  Test  Group  (Programs  1-9),”  ITR-1445,  which  Includes: 

(1)  a  description  of  each  detonation,  including  yield,  aero-point  location  and  environment,  type 
of  device,  ambient  atmospheric  conditions,  etc.;  (2)  a  discussion  of  project  results;  (3)  a  sum¬ 
mary  of  the  objectives  and  results  of  each  project;  and  (4)  a  listing  of  project  reports  for  the 
military-effect  program. 
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Project  5.3  participated  in  Operation  Plumbbob  to:  (1)  measure  thermal  and  blast  response  of 
the  FJ-4  airplane  to  nuclear  explosion  effects;  (2)  correlate  experimental  response  data  with 
the  analytical  predictions  to  confirm  the  delivery  capability  of  the  FJ-4  aircraft;  and  (3)  obtain 
data  to  improve  the  methods  of  predicting  the  blast  response  of  swept-wing  aircraft. 

Two  aircraft,  a  primary  and  a  standby,  were  flown  during  the  operation.  Provisions  were 
included  for  identical  instrumentation  in  the  two  aircraft.  This  Instrumentation  provided  for 
the  acquisition  of  data  on  thermal  input  and  response  of  thin-skin  and  honeycomb  structures, 
nuclear  effects,  overpressures,  dynamic  response,  engine  response,  and  chord-wise  pressure 
distribution  over  the  swept  wing. 

Instrumentation  included  strain  gages,  thermocouples,  calorimeters,  radiometers,  pres¬ 
sure  transducers,  film  badges,  and  pitch-rate  and  attitude  gyros. 

The  aircraft  participated  in  seven  shots,  ranging  in  yields  from  10,3  to  74.1  kilotons.  The 
only  damage  sustained  was  the  blistering  of  the  neoprene  rain-erosion  coating  on  the  nose 
radome  during  Shot  Hood. 

From  die  final  analyses  of  the  test  data  the  following  conclusions  were  made:  (1)  correla¬ 
tion  of  radiant  exposure  data  indicated  that  shielding  affects  the  thermal  output  of  the  nuclear 
device;  (2)  the  methods  for  predicting  the  thermal  response  of  both  thin-skin  and  honeycomb 
panels  to  known  radiant  exposure  and  irradiance  from  a  nuclear  detonation  are  satisfactory; 

(3)  the  critical  dynamic  response  of  the  aircraft  structure  was  determined  to  be  bending  at 
Wing  Station  17.5;  (4)  the  maximum  stress  levels  at  Wing  Station  17.5  and  the  maximum  am¬ 
plitudes  of  airplane  rigid-body  pitching  motion  were  accurately  predicted,  with  the  former  hav¬ 
ing  an  average  error  of  only  2.55  percent  of  limit-allowable  stress;  (5)  for  the  high-blast  in¬ 
cidence  angles  experienced  in  this  operation,  the  blast  effect  and  rigid-body  response  were 
both  nearly  proportional  to  free-stream  overpressure  and  blast-incidence  angle;  (6)  there  was 
a  discernible  diffraction  pulse  of  force  on  the  FJ-4  wing  during  blast  impingement,  however, 
this  pulse  had  a  negligible  effect  on  the  critical  stresses  of  the  FJ-4  wing;  (7)  correlation  of 
gamma  doses  was  satisfactory  when  shielding  was  taken  into  consideration;  (8)  the  correlation 
of  postshot  free-stream  overpressures  and  times-of-shock  arrival  with  measured  values  was 
excellent;  and  (9)  no  power-plant  structural  damage  or  adverse  operational  effects  occurred. 
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Chapter  1 

INTRODUCTION 


1.1  OBJECTIVES 

Two  FJ-4  aircraft  (Figure  1.1)  participated  in  Operation  Plumbbob  to:  (1)  measure  ther¬ 
mal  and  blast  response  of  the  FJ-4  airplane  to  nuclear  explosion  effects;  (2)  correlate  experi¬ 
mental  response  data  with  analytical  predictions  to  confirm  the  delivery  capability  of  the  FJ-4 
aircraft;  and  (3)  obtain  data  to  improve  the  methods  of  predicting  the  blast  response  of  swept  - 
wing  aircraft. 


1.2  BACKGROUND 

With  the  advent  of  nuclear  weapons  in  modern  warfare,  the  problem  of  accurately  defining 
the  capabilities  of  aircraft  to  safely  deliver  such  weapons  became  acute.  The  development  of 
methods  for  predicting  weapon  effects  and  aircraft  responses  to  these  effects  was  severely 
limited  by  a  lack  of  quantitative  data  defining  nuclear  weapon  explosion  phenomena.  Laboratory 
simulation  of  the  problem  was  not  able  to  yield  the  required  data.  Consequently,  full-scale  nu¬ 
clear  test  programs  were  established  in  order  to  provide  the  basic  data  necessary  to  develop 
weapon  effects  prediction  methods  and  delivery  criteria.  Towards  this  end,  numerous  full- 
scale  nuclear  tests,  beginning  with  Operation  Crossroads  in  1946  and  extending  through  Opera¬ 
tion  Redwing  in  1956,  were  conducted  both  at  the  NTS  and  EPG.  These  tests,  in  turn,  defined 
problem  areas  to  be  investigated  in  subsequent  tests  and  confirmed,  corrected,  and  consider¬ 
ably  extended  the  knowledge  of  weapon  phenomena. 

Commensurate  with  progressive  refinement  of  weapon  effects  prediction  methods  it  be¬ 
came  necessary  to  comprehensively  Investigate  and  refine  aircraft  response  prediction  tech¬ 
niques  in  order  to  fully  describe  and  exploit  the  delivery  capabilities  of  modern  aircraft.  In 
addition  to  confirming  the  delivery  capabilities  of  FJ-4/4B  aircraft  with  nuclear  weapons,  the 
participation  of  Project  5.3  in  Operation  Plumbbob  was  specifically  planned  to  evaluate  re¬ 
sponse  theories,  and  in  particular  to  obtain  blast  response  data  for  swept -wing  aircraft  and 
thermal-response  data  for  thin-skinned  aluminum  honeycomb  sandwich  materials. 

The  FJ-4  and  FJ-4B  aircraft  both  have  a  delivery  capability  with  low  yield  nuclear  weap¬ 
ons.  In  addition  the  FJ-4B  has  a  delivery  capability  for  weapons  in  the  megaton-yield  range. 
Since  the  FJ-4  and  FJ-4B  have  Identical  response  characteristics  and  identical  positioning 
limits  with  the  exception  of  the  maximum-allowable  temperature  in  honeycomb  control  sur¬ 
faces,  the  confirmation  of  delivery  capability  using  FJ-4  test  aircraft  also  implied  confirma¬ 
tion  of  the  FJ-4B  delivery  capability.  Therefore,  the  FJ-4  test  aircraft  were  modified  by 
replacing  all  honeycomb  control  surfaces  with  FJ-4B  honeycomb  control  surfaces,  i.e.,  the 
thicknesses  of  face  plates  were  changed  from  0.012  inch  to  0.016  inch  and  the  bond  was  changed 
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from  low  temperature  FM-47  to  high  temperature  HT-424.  This  not  only  permitted  extension 
of  the  data  to  the  FJ-4B  but  also  permitted  the  use  of  a  higher  temperature  limit  for  positioning. 
Operation  Plumbbob  was  a  low-yield  test  series;  however,  it  was  desirable  that  the  FJ-4  par¬ 
ticipate  to  confirm  the  low-yield  delivery  capability  of  the  aircraft  and  to  provide  enough  re¬ 
sponse  data  to  reliably  predict  escape  criteria  and  delivery  capability  for  the  FJ-4B  with 
high-yield  weapons. 


Figure  1.1  FJ-4  airplane  instrumented  for  Operation  Plumbbob. 


1.3  THEORY 

Two  basic  classes  of  data  are  obtained  from  full-scale  nuclear  tests:  (1)  data  that  will 
define  previously  uninvestigated  areas  of  weapon  effects  and  aircraft  response;  and  (2)  data 
that  will  provide  verification  or  modification  of  existing  analytical  prediction  methods  or  sug¬ 
gest  the  development  of  new  methods.  Weapon  effects  which  were  specifically  investigated  or 
which  imposed  critical  limitations  in  the  positioning  of  the  FJ-4  aircraft  in  Operation  Plumb¬ 
bob  were  radiant  exposure,  nuclear  radiation,  overpressures  associated  with  the  shock  wave, 
and  material  velocity  of  the  air  behind  the  shock  wave. 

It  should  be  noted  that  the  problems  of  dictating  safe-escape  requirements  for  aircraft  in 
the  combat  delivery  of  nuclear  weapons  are  basically  the  same  as  the  problems  associated 
with  the  safe  positioning  of  aircraft  in  nuclear  test  operations,  although  the  margins  for  safety 
requirements  may  change. 

1.3,1  Thermal  Effects  and  Response.  The  Special  Weapons  Effects  Branch,  Bureau  of 
Aeronautics,  recommended  the  use  of  a  total  temperature  of  350F  as  a  tentative  limit  for 
aluminum-bonded  honeycomb  sandwich  material  for  positioning  in  Operation  Plumbbob.  Labo¬ 
ratory  tests  and  theoretical  analyses  demonstrated  that  all  structural  components,  Including 
the  honeycomb  material,  were  capable  of  withstanding  total  temperatures  up  to  350F  without 
failure  or  permanent  deformation.  Positioning  calculations  for  thermal  effects  in  Operation 
Plumbbob  were  based  on  the  methods  given  in  Appendix  A.  A  comparison  was  made  of  the  heat¬ 
sink  characteristics  of  honeycomb  and  conventional  skin-stringer  constructions  and  also  of  the 
effects  of  both  flyaway  and  convective  cooling  on  the  temperature  rise  in  each  of  these  struc¬ 
tures.  These  comparisons  showed  that  the  thermal  responses  of  these  structures  were  differ¬ 
ent  and  could  not  be  determined  by  either  a  unified  method  or  by  considering  one  structure  as 
having  a  constant  fraction  of  the  response  of  the  other.  For  these  reasons,  the  thermal  prob- 
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lem  was  divided  into  two  phases,  radiant  exposure  and  thermal  response.  Since  the  test  air¬ 
craft  were  instrumented  for  the  collection  of  both  radiant  exposure  and  temperature  data, 
comparisons  of  the  measured  and  calculated  values  of  both  radiant  exposure  and  thermal  re¬ 
sponse  were  made. 

The  methods  used  in  predicting  the  magnitudes  of  radiant  exposure  and  temperature  re¬ 
sponses  for  honeycomb  and  conventional  structure  are  presented  in  Appendix  A. 

1.3.2  Dynamic  Response.  A  study  was  initiated  to  determine,  by  a  comprehensive  analy¬ 
sis,  the  actual  detailed  motions  of  the  FJ-4  airplane  in  response  to  a  complex  blast  input.  This 
included  not  only  the  motions  of  the  airplane  in  space  but  also  the  transient  vibratory  motions 

of  the  structure.  These  latter  motions  served  as  the  basic  means  of  evaluating  the  stress  levels 
existing  throughout  the  structure.  A  summary  of  the  methods  used  in  the  establishment  of  the 
dynamic -response  positioning  criteria  is  included  as  Appendix  B. 

1.3.3  Nuclear  Radiation.  The  maximum-permissible  nuclear  dose  for  Naval  personnel 
engaged  in  Operation  Plumbbob  was  established  by  the  Chief  of  Naval  Operations  and  concurred 
in  by  Commander,  Field  Command,  AFSWP,  as  follows:  (1)  30  rem  over  2  consecutive  years, 

(2)  15  rem  over  3  consecutive  months,  and  (3)  5  rem  per  event.  Envelopes  were  then  computed, 
delineating  the  spatial  regions  around  ground  zero  in  which  the  pilot  would  receive  a  total  ex¬ 
posure  of  5  rem  or  more.  The  method  used  for  computing  these  envelopes  for  positioning  was 
based  on  weapon-effects  data  obtained  from  Reference  1.  The  methods  used  for  postshot  data 
analysis  were  developed  from  References  2  and  3  and  are  presented  in  Appendix  C. 

1.3.4  Static  Overpressure.  A  study  was  performed  to  determine  the  capability  of  the  FJ-4 
to  withstand  the  crushing  overpressure  associated  with  immersion  in  the  shock  front.  The  study 
consisted  of  a  review  of  previous  FJ-4  pressure  analyses,  an  investigation  of  the  application 

of  impact  loads  (Reference  4),  a  fuselage  frame  analysis,  and  a  skin-panel  and  access-door 
pressure  analysis.  The  results  of  this  study  determined  that  the  canopy  was  the  most  critical 
component,  because  of  the  impact  nature  of  the  applied  force.  The  maximum-allowable  over¬ 
pressure  was  calculated  to  be  3.10  psl.  Envelopes  were  then  computed  defining  the  spatial 
regions  around  ground  zero  in  which  the  canopy  would  fail  because  of  excess  overpressure. 

The  computational  procedures  used  in  the  calculation  of  the  overpressure  envelopes  are  listed 
in  Appendix  D. 

1.3.5  Triple-Point  Path.  The  vicinity  of  the  triple-point  path  still  remains  uninvestigated, 
insofar  as  the  effect  of  this  reinforced  shock  wave  on  aircraft  structure  is  concerned.  Posi¬ 
tioning  safety  criteria  for  full-scale  nuclear  tests  includes  a  requirement  that  test  aircraft  be 
positioned  no  closer  to  the  triple-point  path  than  15  percent  of  the  computed  horizontal  dis¬ 
tance  of  the  path  from  ground  zero.  Consequently,  positioning  envelopes  for  Project  5.3  in¬ 
cluded  overlays  of  the  triple-point  path  and  the  required  safety  margin.  The  path  for  each  yield 
in  the  operation  was  computed  by  means  of  the  method  described  in  Appendix  D. 
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Chapter  2 

PROCEDURE 


2.1  OPERATIONS 

North  American  Aviation,  Inc.,  (NAA),  Columbus  Division,  supplied  two  identically  in¬ 
strumented  FJ-4  airplanes  to  the  Naval  Air  Special  Weapons  Facility  (NASWF)  for  use  during 
Operation  Plumbbob.  In  addition  to  supplying  the  test  aircraft,  NAA  supplied  personnel  and 
materials  to  adequately  maintain  the  airplanes  and  Instrumentation  and  analytical  effort  to  as¬ 
sure  safe  positioning,  collection  of  significant  data,  data  reduction,  and  report  writing.  Project 
5.3  participated  during  Shots  Boltzmann,  Hood,  Diablo,  Kepler,  Shasta,  Doppler  (both  aircraft), 
and  Smoky  (both  aircraft).  During  Shot  Doppler,  the  oscillographs  for  FJ-4,  BUNO 139310  were 
turned  on  11  seconds  late,  which  invalidated  the  data  both  for  time  of  burst  and  time  of  shock 
arrival.  In  Shot  Smoky,  the  M-33  radar  van  broke  lock  on  FJ-4,  BUNO  139467,  forcing  the  ship 
to  abort  the  mission.  The  test  airplanes  were  positioned  in  a  level  flight  altitude,  tail-on  to  the 
blast  at  time  of  shock  arrival  for  all  shots.  Test-site  activity,  other  than  maintenance  of  the 
test  aircraft  and  their  instrumentation,  consisted  of  immediate  reduction  of  data  recorded  from 
each  shot.  Analysis  of  the  raw  data  from  each  shot  was  used  to  review  the  preshot  positioning 
criteria.  Thus,  more  desirable  positioning  for  succeeding  shots  was  accomplished  by  the  ap¬ 
plication  of  the  new  data.  Because  of  the  limited  data  reduction  facilities  at  the  test  site,  con¬ 
tinuous  reduction  and  correlation  of  data  was  accomplished  at  NAA,  Columbus  Division,  and 
continuous  liaison  maintained  with  the  field  participants. 

2.2  DATA  REQUIREMENTS  AND  DATA  REDUCTION 

A  description  of  the  required  data  and  methods  of  data  reduction  follows.  Exact  methods 
of  recording  data  are  fully  described  in  Appendix  E. 

2.2.1  Data  Required.  The  basic  parameters  required  were  divided  into  six  categories: 

(1)  thermal  effects  and  response;  (2)  dynamic  response,  including  flight  and  stability  parame¬ 
ters;  (3)  wing  overpressure  distributions;  (4)  nuclear  radiation;  (5)  free-stream  overpressure; 
and  (6)  engine  response. 

The  thermal  effects  and  response  parameters  needed  from  the  tests  were:  (1)  time  histo¬ 
ries  of  irradiance  and  radiant  exposure;  and  (2)  temperature-time  histories  of  the  critical  skin 
panels  on  the  wings,  flaps,  elevator,  aileron,  and  fuselage. 

The  dynamic  response  parameters  required  from  the  tests  were  time  histories  of  shear 
and  bending-moment  stresses  in  the  wings  and  horizontal  stabilizer.  The  flight  and  stability 
parameters  required  for  analysis  and  correlation  of  theoretical  results  with  test  data  were: 
airspeed;  altitude;  angle  and  rate  of  yaw,  pitch,  and  roll;  angle  of  attack;  normal  and  lateral 
acceleration;  control  surface  position  of  the  ailerons,  elevator,  and  horizontal  stabilizer;  and 
fuei  quantity  remaining  in  the  tanks. 
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For  the  wing  pressure  survey,  time  histories  of  the  pressure  differential  above  and  below 
static  air  load  at  various  chord-wise  positions  on  the  upper  and  lower  surfaces  of  the  wing 
were  required. 

For  determining  pilot  safety,  total  gamma  dose  was  recorded  by  the  use  of  film  badges  in 
the  cockpit,  ammunition  bay  right-drop  tank,  and  nose-wheel  well. 

The  data  requirements  for  a  survey  of  free-stream  overpressure  were  met  by  recording 
time  histories  of  differential  pressure  on  the  nose  boom,  fin  cap,  right  and  left  wing  tips,  and 
the  left  side  of  the  fuselage. 

The  engine  response  parameters  required  were:  compressor  inlet  total  temperature  and 
pressure;  compressor  discharge  total  pressure;  tailpipe  total  temperature  and  pressure;  fuel 
flow;  and  tachometer  reading. 

2.2.2  Data  Reduction.  Standard  data  reduction  procedures  were  used  in  the  compilation 
of  workable  data  from  the  oscillograph  traces  and  photo-recorder  films.  Reduction  of  data 
from  the  calorimeters  and  radiometers  was  accomplished  by  the  Naval  Radiological  Defense 
Laboratory  (NRDL).  All  other  data  were  reduced  by  NAA. 

2.3  DATA  CORRELATION 

The  data  correlation  procedures  used  in  the  final  analysis  of  weapon  effects  input  and  re¬ 
sponse  data  are  presented  in  this  section. 

2.3.1  Thermal  Effects  and  Response  Data  Correlation.  For  purposes  of  data  correlation 
and  analysis,  the  thermal  effects  and  response  data  were  separated  into  two  parts:  radiant- 
exposure  data  and  temperature-response  data.  Measured  initial  conditions  and  final  yields 
were  used  for  all  correlation  and  analysis. 

In  addition  to  comparing  the  calculated  and  measured  values  of  radiant  exposure,  peak  ir- 
radiance,  time  to  peak  irradiance,  and  the  generalized  thermal  pulse,  a  reliability  check  was 
made  on  the  radiant  exposure  data  by  comparing  the  measured  calorimeter  time  histories  to 
the  integrated  measured  radiometer  data. 

In  order  to  check  for  validity  of  the  temperature  response  prediction  methods,  as  sepa¬ 
rated  from  the  radiant  exposure  prediction  method,  the  measured  temperature-time  histories 
for  selected  thermocouples  were  compared  to  corresponding  temperature-time  histories  which 
were  computed  using  measured  radiant  exposure  and  irradiance  data.  The  measured  maximum 
temperature  rise  data  were  also  compared  to  calculated  data  derived  from  calculated  radiant 
exposure  and  irradiance  values. 

In  all  comparisons  of  die  temperature  response  data,  the  theoretical  turbulent  heat  trans¬ 
fer  coefficient  was  used  in  the  computation  of  calculated  values,  although  a  comparison  was 
made  of  this  parameter  with  corresponding  values  derived  from  measured  temperature-time 
histories. 

2.3.2  Dynamic  Response  Data  Correlation.  The  measured  stability  variables  which  de¬ 
scribe  the  rigid-body  response  were  compared  directly  with  the  corresponding  values  calcu¬ 
lated  by  the  methods  discussed  in  Appendix  B.  These  variables,  consisting  of  normal  load 
factor,  angle  of  attack,  pitch  angle,  and  pitch  rate,  were  instrumental  in  the  validation  of  stress 
levels  existing  throughout  the  wing  structure,  since  the  longitudinal  motions  of  the  airplane 
induce  significant  air  loads  on  the  wing. 

The  shear  and  bending-moment  stresses  recorded  by  the  wing  strain  gages  were  corre¬ 
lated  with  the  corresponding  calculated  values  from  the  dynamic  response  analysis.  These 
correlations,  in  particular  those  for  the  bending-moment  stresses  at  the  most  critical  station. 
Wing  Ration  17.5,  were  used  to  modify  and  verify  the  positioning  criteria  during  and  after  the 
test  operation. 

2.3.3  Wing  Pressure  Data  Correlation.  The  dynamic  response  analysis  referred  to  in 
Section  2.3.2  utilized  a  gust  representation  of  the  blast  wave  as  the  primary  excitation  function. 
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To  achieve  a  better  knowledge  of  the  actual  function,  a  static -overpressure  survey  was  made  at 
Wing  Station  175.75.  The  time  histories  of  overpressure  recorded  by  14  pressure  transducers 
were  cross  plotted  at  numerous  points  in  time  following  arrival  of  the  shock  front.  The  re¬ 
sulting  plots  demonstrated  the  chord-wise  load  distribution  which  was  very  useful  in  analyzing 
the  nature  of  the  total  load.  These  cross  plots  were  then  integrated  over  the  chord  to  deter¬ 
mine  section  lift  and  center-of -pressure  location  as  functions  of  time. 
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Chapter  3 

RESULTS 


Included  in  this  chapter  are  positions  of  the  test  aircraft  for  each  shot,  atmospheric  and 
flight  conditions,  reduced  weapon-effects  data  for  each  shot,  and  comparisons  of  calculated  in¬ 
puts  and  responses  with  observed  phenomena. 

Table  3.1  shows  the  positions  of  the  test  aircraft  at  the  time  of  burst  and  time  of  shock 
arrival.  Tables  3.2  through  3.4  are  summaries  of  shot  data  and  atmospheric  and  aircraft- 
flight  conditions.  Tables  3.5  through  3.15  present  the  input  and  response  data  used  for  final 
analyses. 


3.1  THERMAL  EFFECTS  AND  RESPONSE  DATA 

Reduced  thermal  effects  and  response  data  required  for  the  final  analyses  consisted  of 
time  histories  of  readings  from  calorimeters,  radiometers,  and  selected  thermocouples;  max¬ 
imum  readings  and  their  associated  times  from  these  calorimeters,  radiometers  and  selected 
thormocouples;  and  readings  at  the  time  of  shock  arrival  for  selected  thermocouples.  Repre¬ 
sentative  data  for  these  effects  and  responses  are  presented  in  this  section.  Figure  3.1  pre¬ 
sents  the  time  histories  of  radiant  exposure  from  the  90-  and  180-degree  calorimeters  and  the 
integrated  180-degree  radiometer  data  for  Shot  Hood.  Figure  3.2  presents  a  comparison  of  the 
calculated  generalized  thermal  pulse  and  the  measured  pulse  based  on  the  180-degree  radiome¬ 
ter  data  for  Shot  Hood.  Figure  3.3  shows  the  measured  temperature-time  histories  for  the 
fuselage  and  right -outer -wing  panel  thin-skins  for  Shot  Hood.  Figure  3.4  compares  the  calcu¬ 
lated  and  measured  time  histories  of  temperature  for  the  thin-skins  during  Shot  Hood.  Figure 
3.5  presents  the  measured  temperature-time  histories  of  the  lower  face  plates  of  die  honey¬ 
comb  control  surfaces  for  Shot  Hood.  Figures  3.6  through  3.8  compare  the  calculated  and 
measured  time  histories  of  temperature  in  the  honeycomb  lower  face  plates  in  Shot  Hood. 
Figures  3.9  through  3.16  present  the  same  type  of  data  as  Figures  3.1  through  3.8,  respec¬ 
tively,  for  Shot  Smoky.  Table  3.5  compares  the  calculated  and  measured  values  of  radiant  ex¬ 
posure,  peak  lrradlance,  and  time  to  peak  lrradlance  for  all  shots.  Tables  3.6,  3.7,  and  3.8 
present  the  calculated  and  measured  values  of  maximum  temperature  rise,  time -to- maximum 
temperature  rise,  and  temperature  at  time  of  shock  arrival,  respectively,  for  selected  ther¬ 
mocouple  locations  for  all  shots.  Table  3.9  presents  the  theoretical  values  of  the  turbulent 
heat  transfer  coefficient  and  the  values  for  this  parameter  as  calculated  from  measured 
temperature-time  histories.  It  should  be  noted  that  the  theoretical  values  of  the  turbulent 
heat-transfer  coefficient  were  used  in  the  calculation  of  the  temperature  data  shown  in  Tables 
3.8,  3.7,  and  3.8.  Time  histories  of  radiant  exposure,  generalized  lrradlance,  and  temperatures, 
both  calculated  and  measured,  are  presented  in  Appendix  F  for  all  shots  except  Shots  Hood  sad 
Smoky. 
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TABLE  3.1  POSITIONS  OF  TEST  AIRCRAFT  AT  TIME  OF 
BURST  AND  TIME  OF  SHOCK  ARRIVAL 


Shot  , 

Altitude 
Above  Burst 

Horizontal 

Distance* 

Slant  Range 

Angle  of  Incidence 
Above  Horizontal 

ft 

ft 

ft 

deg 

Time  of  Burst: 

Boltzmann 

9,545 

44,650 

10,617 

64.03 

Hood 

14,599 

-3,216 

14,949 

77.60 

Diablo 

9,905 

-1,692 

10,048 

80.30 

Kepler 

9,645 

-214 

9,647 

88.80 

Shasta 

7,092 

+2,719 

7,595 

69.03 

Doppler 

10,905 

-3,957 

11,601 

70.06 

Smoky 

9,992 

+7,518 

12,502 

53.05 

Time  of  Shock  Arrival: 

Boltzmann 

9,710 

+17,545 

20,055 

28.96 

Hood 

14,532 

+5,010 

15,372 

70.97 

Diablo 

9,904 

+4,191 

10,754 

67.07 

Kepler 

9,681 

+6,122 

11,454 

57.69 

Shasta 

7,244 

+9,197 

11,707 

38.23 

Doppler 

10,900 

+2,802 

11,254 

75.59 

Smoky 

10,024 

+20,436 

22,762 

26.13 

*-  Short  at  Ground  Zero;  +  beyond  Ground  Zero. 


TABLE  3.2  SUMMARY  OF  SHOT  DATA 


Shot 

Final 

Yield 

Burst  Height 
Above  Terrain 

Terrain 
Elevation 
Above  MSL 

kt 

ft 

ft 

Boltzmann 

11.5  ±  0.8 

500  Tower 

4,245 

Hood 

74.1  ±  5  pet 

1,500  Balloon 

4,230 

Diablo 

18.7  ±  1.5 

500  Tower 

4,469 

Kepler 

10.3  ±  0.5 

500  Tower 

4,309 

Shasta 

16.8  ±  5  pet 

500  Tower 

4,382 

Doppler 

10.7  ±  0.5 

1,500  Balloon 

4,186 

Smoky 

43.7  ±  5  pet 

700  Tower 

4,479 
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TABLE  3.3  ATMOSPHERIC  CONDITIONS  AT  TIME  OF  SHOT 


Visibility  was  unrestricted  for  all  shots.  Data  was  obtained  from  official  meteorological  reports. 


Shot 

Boltzmann 

Hood 

Diablo 

Kepler 

Shasta 

Doppler 

Smoky 

Date 

38  May 

5  July 

15  July 

24  July 

18  Aug. 

23  Aug. 

31  Aug. 

Pressure  at  Ground  Zero,  mb 

867 

875 

864 

865 

866 

877 

856 

Pressure  at  Flight  Altitude,  mb 

602 

483 

593 

601 

665 

559 

580 

Temperature  at  Ground  Zero,  °F 

65 

70 

74 

70 

80 

71 

57 

Temperature  at  Flight  Altitude,  °F 

27 

18 

38 

34 

51 

26 

20 

Density  at  Flight  Altitude,  slugs/ft* 

0.001521 

0.001225 

0.001410 

0.001488 

0.001581 

0.001401 

0.001469 

Speed  of  Sound  at  Flight  Altitude,  ft/sec 

1082 

1072 

1094 

1090 

1108 

1081 

1074 

Humidity  at  Ground  Zero,  pet 

40 

19 

20 

22 

33 

57 

31 

Humidity  at  Flight  Altitude,  pet 

61 

NAf 

13 

29 

39 

58 

NAt 

Low  Clouds,  Coverage  Type* 

None 

None 

None 

None 

None 

None 

None 

Medium  Clouds,  Coverage  Type* 

4/10  AC 

None 

None 

None 

None 

None 

None 

High  Clouds,  Coverage  Type* 

1/10  Ci 

None 

None 

None 

None 

None 

None 

•  Standard  Meteorological  abbreviations  are  used, 
t  NA:  Data  not  available. 


TABLE  3.4  FLIGHT  CONDITIONS  OF  TEST  AIRCRAFT  AT  TIME  OF 
BURST  AND  TIME  OF  SHOCK  ARRIVAL 


NA,  Data  not  available;  NR,  Data  not  recorded. 


Shot 

True 

Airspeed 

Mach  Number 

Angle  of 
Pitch* 

Gross  Weight 

Angle  of 
Attack* 

ft/sec 

deg 

lb 

deg 

Time  of  Burst: 

Boltzmann 

743 

0.6992 

0 

.... 

Hood 

684 

0.6611 

-0.3 

— 

Diablo 

692 

0.653 

+4.0 

— 

Kepler 

709 

0.669 

+3.1 

— 

— 

Shasta 

700 

0.635 

+2.2 

— 

— 

Doppler 

737 

0.681 

+1.6 

— 

— 

Smoky 

733 

0.679 

-1.5 

— 

— 

Time  of  Shock  Arrival: 

Boltzmann 

753 

0.713 

-- 

18,740 

NR 

Hood 

687 

0.663 

— 

18,834 

+2.7 

Diablo 

689 

0.648 

— 

18,752 

+2.5 

Kepler 

708 

0.663 

— 

18,459 

+1.7 

Shasta 

696 

0.632 

— 

18,059 

+2.0 

Doppler 

737 

0.682 

— 

18,584 

+1.5 

Smoky 

739 

0.682 

— 

18,609 

NA 

*  vindicates  nose  up;  -Indicates  nose  down. 
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lOMilniJ  tw  temperature-time  Matorlea  computed  from  measured  radiometer  data. 

ICUalaed  from  temperature -time  Matorlea  computed  from  differentiated  calorimeter  time  hlatorlea, 
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*  Computed  map  calculated  lrradlaace. 
t  Compded  ealap  meeaured  lrradlaace. 

t  Comfut*:;  oaiaf  value  of  it  attained  top  differentiating  colorimeter  time  hletorlee. 

I  Ottalaed  from  temperature-time  hletorlee  computed  from  meeeured  lrredlence  date. 


TABLE  3.8  CALCULATED  AND  MEASURED  VALUES  OF  TEMPERATURE  AT  TIME  OF  STOCK  ARRIVAL  FOR 
SELECTED  THERMOCOUPLES 


NA,  Data  not  available;  NR,  Data  not  recorded. 


Shot 

Measured  Time 
of 

Shock  Arrival 

Fuselage 

Station 

81.378 

Fuselage 

310.81 

Fuselage 

Station 

348.188 

Fuselage 

OtmHrm, 

333.638 

Fuselage 

phH«. 

388.138 

Wing 

Station 

170.0 

Wing 

Station 

333.138 

Wing 

Station 

333.78 

Bee 

“F 

°F 

°F 

•F 

°F 

°F 

"F 

”F 

Boltzmann 

Calculated* 

NA 

73.7 

68.6 

78.4 

88.4 

NA 

67.1 

64.3 

Calculated! 

NA 

68.0 

64.3 

78.8 

87.6 

NA 

64.1 

88.3 

Measured 

16.76 

NR 

67.8 

67.0 

78.6 

88.7 

NR 

64.2 

81.3 

Hood 

Calculated* 

NA 

100.8 

104.4 

118.1 

164.0 

NA 

83.8 

NA 

Calculated! 

NA 

60.8 

86.3 

108.8 

144.8 

NA 

67.6 

NA 

Measured 

11.06 

NR 

67.8 

81.6 

131.8 

148.6 

NR 

61.3 

NR 

Diablo 

Calculated* 

NA 

118.6 

118.3 

130.4 

180.4 

NA 

NA 

103.0 

Calculated! 

NA 

66.6 

63.8 

104.8 

116.8 

NA 

NA 

81.8 

Measured 

8.06 

NR 

66.0 

66.6 

107.8 

130.1 

NR 

NR 

87.4 

Kepler 

Calculated* 

NA 

63.6 

60.4 

104.8 

118.6 

131.0 

NA 

81.1 

Calculated!,} 

NA 

76.8 

78.4 

67.4 

100.0 

63.0 

NA 

68.0 

Measured 

8.66 

NA 

61.0 

77.8 

61.0 

101.4 

83.7 

NR 

73.7 

Shasta 

Calculated* 

138.7 

188.1 

140.0 

146.7 

178.0 

306.6 

NA 

118.8 

Calculated! 

66.6 

101.8 

104.0 

113.3 

138.4 

116.8 

NA 

86.8 

Measured 

8.88 

104.0 

106.4 

106.3 

118.3 

130.3 

117.0 

NR 

88.0 

Doppler 

Calculated* 

74.1 

71.7 

74.8 

64.8 

06.4 

66.1 

NA 

08.0 

Calculated! 

78.1 

68.3 

66.6 

64.8 

04.4 

77.7 

NA 

63.0 

Measured 

8.66 

78.3 

73.8 

78.4 

67.1 

117.1 

78.1 

NR 

66.0 

Smoky 

Calculated* 

83.3 

77.4 

83.4 

NA 

183.8 

100.7 

NA 

70.6 

Calculated! 

66.4 

71.3 

76.4 

NA 

136.3 

76.8 

NA 

60.7 

Measured 

17.83 

73.3 

78.7 

78.1 

NA 

111.0 

71.0 

NR 

60.3 

‘Computed  using  calculated  radiant  exposure  and  calculated  lrradt— re. 

^Computed  using  measured  radiant  exposure  and  measured  lr radiance. 

t  Computed  using  measured  radiant  eaposure  aad  irradUnce  oMaiaed  bp  dtBerei mating  calorimeter  data. 
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TABLE  3.9  CALCULATED  AND  THEORETICAL  VALUES  OF  THE  TURBULENT  HEAT 
TRANSFER  COEFFICIENT  FOR  SELECTED  THERMOCOUPLE  LOCATIONS 


Shot 

Fuselage  or 
Wing  Station* 

Calculated  From 
Skin  Temperature 
Data 

Theoretical 

Percent  Difference 
Based  on  Skin 
remp.  Calculations 

cal/cmVsec/°F 

cal/cmVsec/°F 

Boltzmann 

F, 248.155 

0.00297 

0.00260 

12.5 

F, 333.625 

0.00303 

0.00248 

18.2 

F,389.125 

0.00400 

0.00240 

40.0 

Hood 

F, 210.81 

0.00085 

0.00223 

162.4 

F, 333.625 

0.00141 

0.00203 

43.9 

F, 389.125 

0.00183 

0.00197 

7.65 

Diablo 

F, 248.155 

0.00231 

0.00254 

9.95 

F, 333.625 

0.00216 

0.00242 

12.0 

F, 389.125 

0.00239 

0.00234 

2.1 

Kepler 

F, 333.625 

0.00132 

0.00245 

93.2 

F, 389.125 

0.00444 

0.00237 

46.6 

W,223.75 

0.00213 

0.00397 

86.4 

Shasta 

F, 210.81 

0.00246 

0.00289 

17.4 

F, 248.155 

0.00302 

0.00276 

8.6 

F, 333.625 

0.00220 

0.00264 

20.0 

F.389.125 

0.00232 

0.00255 

9.91 

W,  170.0 

0.00298 

0.00426 

42.9 

Doppler 

F, 389. 125 

0.00183 

0.00228 

24.5 

W.170.0 

0.00325 

0.00382 

17.5 

Smoky 

F.210.81 

0.00258 

0.00267 

3.48 

F, 248. 155 

0.00255 

0.00255 

0.00 

F, 389.125 

0.00357 

0.00235 

34.2 

W,  223.75 

0.0069 

0.00394 

42.9 

W.170.0 

0.0036 

0.00394 

9.44 

F.81.375 

0.00341 

0.00323 

5.28 

*  F  designates  Fuselage  Station;  W  designates  Wing  Station. 


TABLE  3.10  CALCULATED  VALUES  OF  SHOCK  WAVE  PARAMETERS 
BASED  ON  FINAL  YIELDS 


Shot 

Peak  Material 
Velocity 

Shock  Front 
Propagation 
Velocity 

Peak  Static 
Overpressure 

Duration  of 
Positive  Phase 

ft/sec 

ft/sec 

psi 

sec 

Boltzmann 

31.80 

1101 

.37 

1.44 

Hood 

76.36 

1119 

.73 

1.81 

Diablo 

63.07 

1132 

.72 

1.12 

Kepler 

45.49 

1118 

.52 

1.01 

Shasta 

70.95 

1151 

.89 

1.26 

Doppler 

48.10 

1110 

.52 

1.04 

Smoky 

47.06 

1102 

.53 

2.09 
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TABLE  3.11  CALCULATED  AND  MEASURED  MAXIMUM  VALUES  OF  AIRPLANE  FLIGHT  AND 
STABILITY  VARIABLES  AFTER  SHOCK  ARRIVAL 


NR,  Data  not  recorded. 


Shot 

Normal 

Load 

Factor 

Perturbation 
of  Angie 
of  Attack 

Rate  of 
Change 
of  Angle 
of  Attack 

Perturbation 
of  Angle 
of  Pitch 

Rate  of 
Change 
of  Angle 
of  Pitch 

Acceleration 
of  Angle 
of  Pitch 

gravities 

degrees 

degrees/sec 

degrees 

degrees/sec 

degrees/sec* 

Boltzmann 

Calculated* 

1.37 

1.03 

-4.58 

-2.23 

-4.13 

-25.0 

Calculated! 

1.49 

1.22 

-3.09 

-1.27 

-2.60 

-11.6 

Measured 

1.54 

NR 

NR 

-1.0 

-3.8 

NR 

Hood 

Calculated* 

2.70 

5.29 

-18.59 

-10.08 

-17.61 

-97.1 

Calculated! 

2.82 

6.32 

-12.15 

-5.39 

-10.71 

—40.4 

Measured 

3.85 

5.9 

NR 

-5.1 

-12.1 

NR 

Diablo 

Calculated* 

2.50 

4.58 

-14.98 

-7.58 

-14.32 

-93.2 

Calculated! 

2.61 

5.00 

-9.74 

-3.87 

-7.97 

-35.0 

Measured 

3.85 

5.4 

NR 

-4.5 

-11.1 

NR 

Kepler 

Calculated* 

1.85 

3.41 

-11.60 

-5.88 

-10.70 

-84.1 

Calculated! 

2.17 

3.22 

-8.70 

-2.42 

-5.39 

-24.9 

Measured 

2.63 

2.8 

NR 

-2.8 

-7.9 

NR 

Shasta 

Calculated* 

2.02 

3.51 

-13.08 

-7.38 

-12.37 

-70.9 

Calculated! 

2.30 

3.92 

-8.48 

-3.42 

-6.76 

-29.8 

Measured 

2.63 

3.0 

NR 

-3.3 

-7.2 

NR 

Doppler 

Calculated* 

2.00 

3.87 

-13.17 

-6.65 

-12.37 

-72.7 

Calculated! 

2.38 

3.68 

-7.56 

-2.84 

-6.25 

-29.8 

Measured 

2.89 

3.4 

NR 

-2.6 

-8.3 

NR 

Smoky 

Calculated* 

1.47 

1.45 

-6.11 

-3.31 

-5.76 

-30.2 

Calculated! 

1.61 

1.70 

-4.18 

-1.86 

-3.58 

-15.2 

Measured 

1.77 

NR 

NR 

-2.0 

-4.9 

NR 

*  Calculated  from  original  equations, 
t  Calculated  from  modified  equations. 
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f Calculated  from  modified  equations. 


TABLE  3.13  MEASURED  MAXIMUM  VALUES  OF  HORIZONTAL 
STABILIZER  BENDING-MOMENT  STRESS  LEVELS 
AFTER  SHOCK  ARRIVAL 


NA,  Data  not  available. 


Shot 

Horizontal 
Stabilizer 
Station  33.75 
(Right) 

Horizontal 
Stabilizer 
Station  33.75 
(Left) 

Horizontal 
Stabilizer 
Station  60.5 
(Left) 

psi 

psi 

psi 

Initial  Stress 

1,720 

1,720 

1,730 

Limit  Allowable  Stress 

40,700 

40,700 

36,700 

Boltzmann 

2,680 

NA 

2,147 

Hood 

6,719 

6,715 

3,332 

Diablo 

6,770 

6,010 

3,234 

Kepler 

4,455 

4,137 

2,512 

Shasta 

4,484 

3,861 

2,531 

Doppler 

3,546 

2,598 

1,941 

Smoky 

3,178 

2,293 

NA 

TABLE  3.14  COMPARISON  OF  CALCULATED  AND  MEASURED  VALUES  OF  NUCLEAR  DOSE 


NA,  Data  Not  Available. 


Shot 

Cockpit 

* 

Measured  Gamma  Dose 

Cockpit  Right  Ammo 

t  Wing  Tank  Bay 

Wheel 

Well 

Postshot 
Calculated 
Gamma  Dose 

Postshot 
Calculated 
Total  Dose 

Shielding 
of  Shot 

r 

r 

r 

r 

r 

r 

rem 

Boltzmann 

1.02 

0.4  to  0.9 

1.15 

1.56 

0.76 

1.27 

1.48 

Light 

Hood 

1.10 

0.76 

1.32 

1.47 

1.04 

1.721 

1.75 

None 

Diablo 

0.70 

0.35 

0.83 

1.10 

0.59 

1.365,1 

1.441 

Heavy  on  one 

side 

Kepler 

0.58 

0.44 

0.48 

0.64 

0.32 

0.281 

0.351 

Very  heavy 

Shasta 

1.90 

NA 

NA 

3.90 

NA 

2.891 

4.531 

Heavy  on  one 

side 

Dopplert 

0.00 

0.90 

1.10 

1.32 

0.75 

1.751 

1.85 

None 

Doppler, 

0.90 

0.77 

1.40 

2.90 

2.90 

1.751 

1.85 

NOne 

Smoky 

0.71 

0.35 

1.30 

1.30 

1.002 

1.15 

1.25 

None 

•  Left  side  of  pilot's  Mat.  1  FJ-4  BUNO  13*467 

t  Aft  of  pilot's  right  shoulder.  3  FJ-4  BUNO  130310 

1  Externally  mounted  on  left  tank. 

I  One  roentgen  added  to  account  (or  aircraft  being  more  than  1000  feet  short  of  ground  sero  at 
time  sero. 

t  Values  shown  have  been  reduced  by  a  (actor  of  10  to  account  for  heavy  shielding. 
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Figure  3.1  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree  calo¬ 
rimeters  and  integrated  180  degree  radiometer,  Shot  Hood. 


Figure  3.2  Comparison  of  the  calculated  and  measured  generalized  thermal  pulse  taken  from 
180  degree  radiometer,  Shot  Hood. 
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Figure  3.3  Temperature  versus  time  data  of  the  thin  skins  on  the  fuselage 
and  right  outer  wing  panel,  Shot  Hood. 
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Figure  3.4  Comparison  of  the  calculated  and  measured  temperature -time 
histories  of  selected  thin  skins,  Shot  Hood. 


37 

CONFIDENTIAL 


o 


0  2  4  6  8  10  12 

Time,  seconds 


Figure  3.5  Temperature  versus  time  data  of  the  lower  faceplates  of  the 
honeycomb  control  surfaces,  Shot  Sood. 


38 

CONFIDENTIAL 


0  24  6  8  10  12 

Time,  seconds 

Figure  3.6  Comparison  of  the  calculated  and  measured  temperature-time 
history  of  the  elevator  honeycomb  faceplate,  Shot  Hood. 


I 
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Figure  3.7  Comparison  of  the  calculated  and  measured  temperature-time 
history  of  the  aileron  honeycomb  faceplate,  Shot  Hood. 
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Figure  3.8  Comparison  of  the  calculated  and  measured  temperature-time 
history  of  the  flap  honeycomb  faceplate,  Shot  Hood. 
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Figure  3.9  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree  calorim¬ 
eters  and  integrated  180  degree  radiometer,  Shot  Smoky. 


Figure  3.10  Comparison  of  the  calculated  and  measured  generalised 
thermal  pulse  taken  from  180  degree  radiometer.  Shot  Smoky. 
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Figure  3.11  Temperature  versus  time  data  of  the  thin  skins  on  the  fuselage 
and  right  outer  wing  panel,  Shot  Smoky. 
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Figure  3.12  Comparison  of  the  calculated  and  measured  temperature-time  his¬ 
tories  of  selected  thin  skins,  Shot  Smoky. 
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Figure  3.13  Temperature  versus  time  data  of  the  lower  faceplates  of  the 
honeycomb  control  surfaces,  Shot  Smoky. 
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Figure  3.14  Comparison  of  the  calculated  and  measured  temperature-time 
history  of  the  elevator  honeycomb  faceplate,  Shot  Smoky. 
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Figure  3. 15  Comparison  of  the  calculated  and  measured  temperature-time 
history  of  the  aileron  honeycomb  faceplate,  Shot  Smoky. 
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Figure  3.16  Comparison  of  the  calculated  and  measured  temperature-time 
history  of  the  flap  honeycomb  faceplate,  Shot  Smoky. 
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3.2  DYNAMIC  RESPONSE  DATA 


The  dynamic  response  data  measured  during  the  test  operations  consisted  of  basic  flight 
and  stability  parameters  and  representative  stresses.  In  addition,  the  blast  wave  static  over¬ 
pressure  was  measured  as  an  initial  means  of  determining  the  applied  air  loads  which  Induce 
the  response. 

Table  3.10  presents  the  values  of  peak  overpressure  behind  the  shock  front  which  were 
calculated  from  the  final  yields.  This  table  also  lists  computed  values  of  other  blast  wave  pa¬ 
rameters  which  were  utilized  in  the  dynamic  response  analysis. 

Table  3.11  itemizes  the  calculated  and  measured  maximum  values  of  airplane  flight  and 
stability  variables  following  impingement  of  the  blast  wave. 

Table  3.12  presents  the  calculated  and  measured  maximum  values  of  shear  and  bending- 
moment  stresses  in  the  wing  structure.  The  corresponding  limit-allowable  and  initial  stress 
levels  are  also  tabulated. 

Table  3.13  itemizes  the  measured  maximum  values  of  representative  bending-moment 
stresses  in  the  horizontal  stabilizer  structure.  The  corresponding  limit -allowable  and  initial- 
stress  levels  are  also  listed. 

Figures  3.17  through  3.21  illustrate  the  dynamic  response  of  the  airplane  as  a  rigid-body 
during  Shot  Shasta.  This  response  is  described  by  the  calculated  and  measured  stability  varia¬ 
bles  as  perturbations  about  their  corresponding  values  existing  prior  to  shock  arrival. 

Figures  3.22  through  3.33  illustrate  the  structural  response  of  the  wing  during  Shot  Shasta. 
This  response  is  described  by  the  calculated  and  measured  variations  of  shear  and  bending- 
moment  stresses  for  selected  points  in  the  wing  structure. 

Additional  dynamic  and  structural  response  data  is  presented  in  Appendix  G  for  Shots 
Boltzmann,  Hood,  Diablo,  Kepler,  Doppler,  and  Smoky. 


3.3  WING  PRESSURE  DATA 

The  wing  pressure  survey  conducted  during  the  test  operation  consisted  of  static  over¬ 
pressure  measured  at  14  points  on  the  upper  and  lower  surfaces  of  the  chord  at  Wing  Station 
175.75,  as  shown  in  Figure  3.34.  This  data  was  recorded  during  and  after  arrival  of  the 
shock  front. 

Figures  3.35  and  3.36  illustrate  the  variations  of  static  overpressure  as  measured  by  each 
transducer  during  Shot  Shasta. 

Figures  3.37  through  3.39  illustrate  the  cross-plotted  values  of  overpressure  for  numerous 
points  in  time  following  arrival  of  the  shock  front  during  Shot  Shasta.  These  values  are  pre¬ 
sented  as  the  chord-wise  distributions  of  positive  overpressure  along  the  upper  and  lower  wing 
surfaces. 

Figures  3.40  and  3.41  illustrate  the  variation  of  section  lift  and  center  of  pressure  follow¬ 
ing  arrival  of  the  blast  wave  for  Shot  Shasta.  The  center -of -pressure  location  of  the  perturbed 
air  load  is  described  by  units  of  percent  chord  aft  of  the  leading  edge.  Corresponding  section- 
lift  and  center-of-pressure  data  are  presented  in  Appendix  H  for  Shots  Hood,  Diablo,  Kepler, 
and  Doppler. 

Figures  3.42  and  3.43  illustrate  the  computer  simulation  of  measured  section  lift  and  the 
corresponding  effect  on  the  maximum  critical-stress  level  for  Shot  Hood. 


3.4  NUCLEAR  RADIATION  DATA 

The  nuclear  radiation  data  obtained  consisted  of  gamma  dose  measurements  from  film 
badges  located  at  four  positions  in  the  test  aircraft.  The  readings  of  these  film  badges  and  the 
postshot  calculated  values  of  gamma  dose  for  each  shot  are  listed  in  Table  3.14.  Postshot  cal¬ 
culated  values  for  total  nuclear  dose  are  also  included  to  show  the  additive  effects  of  neutron 
radiation. 
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Figure  3.17  Comparison  of  calculated  and  measured  gust  Induced  perturba¬ 
tion  of  normal  load  factor  versus  time.  Shot  Shasta. 
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Figure  3.18  Comparison  of  calculated  and  measured  gust  induced  perturba' 
tion  of  pitch  angle  versus  time,  Shot  Shasta. 
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Figure  3.19  Comparison  of  calculated  and  measured  pitch  rate  versus  time. 
Shot  Shasta. 


Figure  3.20  Measured  gust  induced  perturbation  of  angle  of  attack  versus 
time.  Shot  Shasta. 
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Figure  3.21  Measured  gust  induced  perturbation  of  horizontal  stabilizer 
position  versus  time,  Shot  Shasta. 
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Figure  3.22  Comparison  of  calculated  and  measured  gust  induced  incremental 
iring  bending  moment  stress  level  at  Right  Wing  Station  17.5  versus  time,  Shot 
Shasta. 
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Figure  3.23  Comparison  of  calculated  and  measured  gust  induced  incremental 
wing  bending  moment  stress  level  at  Left  Wing  Station  17.5  versus  time,  Shot 
Shasta. 


Figure  3.24  Comparison  of  calculated  and  measured  gust  induced  Incremental 
wing  bending  moment  stress  level  at  Wing  Station  36.5  versus  time.  Shot 
Shasta. 
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Figure  3.25  Comparison  of  calculated  and  measured  gust  Induced  incremental 
wing  bending  moment  stress  level  at  Hear  Spar  Station  173  versus  time.  Shot 
Shasta. 


Figure  3.26  Comparison  of  calculated  and  measured  gust  Induced  incremental 
wing  shear  stress  level  at  Wing  Station  48  versus  time,  Shot  Shasta. 
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Figure  3.27  Comparison  of  calculated  and  measured  gust  Induced  incre¬ 
mental  wing  shear  stress  level  at  Rear  Spar  Station  50.987  versus  time.  Shot 
Shasta. 
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Figure  3.28  Comparison  of  calculated  and  measured  gust  Induced  incre¬ 
mental  wing  shear  stress  level  at  Rear  Spar  Station  83.012  versus  time.  Shot 
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Figure  3.31  Comparison  of  calculated  and  measured  gust  induced  Incremental 
■ting  shear  stress  level  at  Rear  Spar  Station  86.2  versus  time.  Shot  Shasta. 
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Figure  3.32  Comparison  of  calculated  and  measured  gust  Induced  incremental 
wing  shear  street  level  at  Rear  Spar  Station  115.129  versus  time,  Shot  Sheets. 
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Figure  3.33  Comparison  of  calculated  and  measured  gust  induced  incremental 
wing  shear  stress  level  at  Rear  Spar  Station  143.25  versus  time.  Shot  Shasta. 


Figure  3.34  Chord -wise  positions  of  wing  pressure  transducers  at  Wing  Station  175.75. 
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Figure  3.35  Measured  static  overpressure  at  Wing  Station  175.75, 
upper  chord,  versus  time,  Shot  Shasta. 
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Figure  3.36  Measured  static  overpressure  at  Wing  Station  175.75, 
lover  chord,  versus  time,  Shot  Shasta. 
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Figure  3.37  Measured  chord-wise  distribution  of  static  overpressure  at 
Wing  Station  175.75  from  2  milliseconds  to  10  milliseconds  after  shock 
arrival,  Shot  Shasta. 
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Figure  3.38  Measured  chord-wise  distribution  of  static  overpressure  at  Wing 
Station  175.75  from  12  milliseconds  to  22  milliseconds  after  shock  arrival, 

Shot  Shasta. 
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Figure  3.39  Measured  chord-wise  distribution  of  static  overpressure  at 
Wing  Station  175.75  from  24  milliseconds  to  50  milliseconds  after  shock 
arrival,  Shot  Shasta. 
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Figure  3.40  Measured  section  lift  at  Wing  Station  175.75  versus  time,  Shot  Shasta. 
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Figure  3.41  Measured  section  center  of  pressure  at  Wing  Station  175.75 
versus  time,  Shot  Shasta. 
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Figure  3.42  Measured  section  lift  and  computer  simulation  of  measured  section  lift 
versus  time.  Shot  Hood. 
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Figure  3.43  Effect  of  peak  diffraction  load  on  bending-moment  stress  levels  at  Wing 
Station  17.5,  Shot  Hood. 

65 


CONFIDENTIAL 


3.5  STATIC  OVERPRESSURE  AND  TIME-Or-SHOCK- ARRIVAL  DATA 


In  addition  to  the  overpressure  data  obtained  In  the  wing  pressure  survey,  measurements 
were  also  msde  at  the  peak  static  overpressure  experienced  at  various  other  locations  on  the 
aircraft  Table  3.15  presents  the  maximum  values  of  measured  overpressure  tor  selected 
transducer  locations  and  a  comparison  of  the  calculated*  and  average-measured  peak-static 
overpressures  for  each  shot  The  average-measured  peak  overpressure  values  from  Table 
3.15  were  reduced  to  correspond  to  a  1  kt  burst  In  a  homogeneous  sea-level  atmosphere  by  use 
at  the  scaling  equations  outlined  In  Appendix  D.  These  scaled  values  are  shown  in  Figure  3.44 
and  are  compared  with  the  peak  overpressure  versus  slant  range  curve  f or  a  1  kt  burst  In  a 
homogeneous  sea-level  atmosphere  obtained  from  Reference  5. 

The  tlme-of-shock-front -arrival  data  presented  In  Table  3.8  were  reduced  to  correspond 
to  a  1  kt  burst  In  a  homogeneous  sea-level  atmosphere  by  use  of  the  scaling  equations  outlined 
In  Appendix  B.  These  scaled  values  are  shown  In  the  form  of  data  points  In  Figure  3.45.  The 
basic  curve  shown  In  this  figure  was  taken  from  Reference  1. 

3.6  ENGINE  RESPONSE  DATA 

A  review  of  all  of  the  engine  response  data  collected  Indicated  similar  responses  in  all 
shots.  The  data  from  Shot  Boltzmann  was  selected  for  presentation  in  this  section  as  being 
representative  at  power  plant  response.  Figures  3.46,  3.47,  and  3.48  present  the  time  histo¬ 
ries  of  compressor  inlet  total  pressure  and  temperature,  compressor  discharge  total  pres¬ 
sure,  tailpipe  total  temperature,  engine  speed,  fuel  flow  rate,  and  calibrated  airspeed  for  Shot 
Boltamaim .  Appendix  I  presents  similar  data  recorded  for  the  other  shots. 
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Figure  3.44  Correlation  of  static  overpressure  data  scaled  to  one  Uloton, 
homogeneous  sea  level  atmospheric  conditions. 
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Figure  3.47  Engine  compressor  inlet  total  temperature,  tailpipe  total  temperature,  and  com¬ 
pressor  discharge  total  pressure  time  histories.  Shot  Boltzmann. 
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Figure  3.48  Calibrated  airspeed,  engine  speed,  and  fuel  flow  time  histories,  Shot  Boltzmann. 
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Chapter  4 

DISCUSSION 


4.1  THERMAL  EFFECTS  AND  RESPONSE 

4.1.1  Measured  Response.  The  data  used  in  the  thermal  effects  and  response  analysis 
were  obtained  from  calorimeters,  radiometers,  and  thermocouples.  The  calorimeters  and 
radiometers  were  supplied,  calibrated,  and  the  data  from  them  reduced  by  NRDL.  The  thermo* 
couples  were  designed  specifically  for  recording  transient  temperatures.  Therefore,  the  in¬ 
stallation  technique  described  in  Appendix  E  was  used  to  increase  the  thermocouple  response 
to  such  a  degree  that  lag  was  considered  negligible.  The  installation  technique,  duplication  of 
thermocouples,  and  careful  preshot  inspections  were  considered  responsible  for  the  high  de¬ 
gree  of  thermocouple  reliability.  Only  one  thermocouple  malfunction  was  observed  during  the 
test.  This  was  at  Fuselage  Station  333.625  during  Shot  Smoky. 

The  only  damage  sustained  by  the  test  aircraft  was  the  blistering  of  the  neoprene  rain- 
erosion  coating  on  the  fiberglas  nose  radome  during  Shot  Hood,  as  shown  in  Figure  4.1.  The 
aircraft  was  positioned  short  of  ground  aero  at  the  time  of  burst  and  received  approximately 
10  cal/  cm1.  The  damaged  radome  was  sent  to  the  Naval  Air  Development  Center,  (NADC), 
Johnsville,  Pennsylvania,  for  evaluation  of  the  effect  of  the  damage  upon  the  electrical  per¬ 
formance.  As  a  result  of  this  evaluation  it  was  concluded  that  negligible  deterioration  in  the 
radome  electrical  performance  resulted  from  the  blistered  neoprene  rain-erosion  coating. 

4.1.2  Radiant  Exposure  and  Irradiance.  The  measured  and  calculated  values  of  radiant 
exposure,  second  peak  irradiance,  and  time  to  the  second  peak  irradiance  are  compared  in 
Table  3.5.  For  purposes  at  data  comparison  and  thermal  response  predictions,  the  values  of 
radiant  exposure  were  obtained  from  the  calorimeter  and  integrated  radiometer  curves  at  a 
value  at  ten  times  the  time  to  second  peak  irradiance  (10  tj)  and  multiplied  by  1.04.  This  pro¬ 
cedure  was  in  accordance  with  the  assumption  in  Reference  1  that  96  percent  of  the  total 
thermal  energy  is  emitted  by  the  time  10  q.  The  thermal  response  prediction  methods  are 
also  based  on  this  assumption.  The  calculated  values  of  radiant  exposure  were  computed  from 
the  methods  in  Appendix  A,  Section  A.1;  and  peak  irradiance  from  Appendix  A,  Section  A.2. 

The  calculations  were  made  in  accordance  with  the  orientation  of  the  calorimeters  and  radi¬ 
ometers  (Table  3.5).  The  values  of  the  more  pertinent  factors  used  in  Equations  A.1  and  A.2 
are  presented  in  Table  4.1.  A'^^.n /A' factors  which  convert  direct  radiant  exposure 
values  to  normal  values  for  use  in  the  temperature  response  analysis  are  also  presented  in 
Table  4.1. 

In  comparing  the  measured  and  calculated  radiant  exposure,  three  measured  values  were 
neglected  because  they  were  inconsistent  with  the  other  values  for  that  particular  event.  These 
were  the  90-degree  calorimeter  from  Shot  Diablo,  180-degree  calorimeter  from  Shot  Shasta, 
and  the  integrated  180-degree  radiometer  from  Shot  Smoky.  The  percentage  of  over  prediction 
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Figure  4.1  FJ-4  radome  showing  thermal  damage  incurred  to  the  neoprene  rain  erosion  coating  during 
Shot  Hood. 


TABLE  4.1  FACTORS  USED  IN  THE  CALCULATION  OF 
THE  RADIANT  EXPOSURE  AND  THE 
IRRADIANCE  VALUES  IN  TABLE  3.5 


Symbols  are  defined  In  Appendix  A.  NA,  not  applicable. 
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was  calculated  using  the  average  of  the  remaining  measured  values,  and  it  was  observed  that 
the  percentage  of  error  increased  with  the  degree  of  shielding  of  the  nuclear  device.  The  two 
shots  known  to  be  unshielded.  Shot  Hood  and  Shot  Smoky,  resulted  in  radiant  exposure  over 
predictions  of  37.5  percent  and  26.4  percent,  respectively,  while  Shot  Kepler,  which  was  heavily 
shielded,  resulted  in  an  over  prediction  of  281.4  percent.  Shot  Boltzmann  had  light  shielding, 
and  the  error  was  48.4  percent.  Shots  Diablo  and  Shasta  were  shielded  heavily  on  one  side, 
resulting  in  over  predictions  of  111.1  percent  and  199.6  percent,  respectively.  Shot  Doppler, 
which  is  over -predicted  by  50.9  percent,  is  considered  separately  since  the  effect  of  Its 
«hi«»iding  is  questionable.  The  comparison  of  calculated  to  measured  radiant  exposure  shows 
that  shioiding  does  lower  the  thermal  yield  and  that  the  values  for  the  unshielded  shots  fall 
within  the  expected  realm  of  accuracy  of  Equation  A.1.  Refinements  to  Equation  A.1  could  pos¬ 
sibly  be  made  by  reviewing  thermal  yield  as  a  function  of  radiochemical  yield,  considering 
atmospheric  attenuation  as  a  function  of  the  spectral  distribution  emitted  by  the  fireball,  in¬ 
vestigating  the  effect  an  ground  albedo  of  the  dust  layer  stirred  up  by  the  shock  wave,  revising 
the  generalized  thermal  pulse  to  the  time  10  ij,  and  reviewing  the  effects  of  fireball  distortion 
an  direct  and  reflected  radiant  energy. 

Comparisons  of  the  generalized  thermal  pulse  for  each  event  were  made  by  non- 
dumensionalizing  radiometer  data  and  calculated  ir radiance  from  Equation  A. 2.  Excluding 
Shot  Kepler,  for  which  no  radiometer  data  was  available,  it  is  observed  that  the  generalized 
thermal  pulses  for  the  shielded  shots  are  more  full  after  peak  irradiance  than  for  the  un¬ 
shielded  shots.  Also,  the  generalized  thermal  pulses  for  the  unshielded  shots  fall  within  the 
expected  realm  of  accuracy  of  ±15  percent.  A  discrepancy  in  the  decay  portion  of  the  meas¬ 
ured  thermal  pulse  from  Shot  Hood  is  noted.  This  decay  portion  is  considered  to  be  in  error 
because  the  comparison  of  measured  to  calculated  temperature -time  histories  for  selected 
thermocouples  in  Figures  3.4,  3.6,  3.7,  and  3.8  indicate  that  the  thermocouple  locations  did 
not  receive  the  additional  irradiance.  A  reasonable  explanation  to  account  for  this  is  that  the 
bottom  skin  of  the  wing  tank,  which  was  at  an  elevated  temperature,  was  in  the  field  of  view  of 
the  radiometer. 

4.1.3  Maximum  Temperature  Rise.  Three  independent  values  of  the  maximum  tempera¬ 
ture  rise  for  selected  thermocouple  locations  were  calculated  for  comparison  with  the  meas¬ 
ured  maximum  temperature  rise.  These  comparisons  are  presented  in  Table  3.6.  The  maxi¬ 
mum  temperature  rises  were  calculated  by  using  the  normalized  value  of  calculated  radiant 
exposure  (Table  3. 5 j  in  Equation  A.3,  using  the  normalized  value  of  average  measured  radiant 
exposure  Table  3.5;  in  Equation  A.3,  and  computing  the  temperature-time  history  from  Equa¬ 
tion  A.  12  using  measured  normalized  irradiance. 
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The  skin-absorption  coefficients  which  were  used  in  the  above  calculations  are  presented 
in  Table  4.2.  Since  a  small  variation  in  magnitude  of  the  absorptivity  can  result  in  an  appre¬ 
ciable  variation  in  temperature  rise  predictions,  a  white  90-degree  calorimeter  waa  installed 
in  the  wing  tank  to  measure  the  white  paint's  absorptivity;  however,  an  unreliable  value  of  0.65 
was  indicated.  Reference  6  shows  that  the  absorptivity  of  a  similar  white  paint  is  0.70  for 
monochromatic  incident  radiation  having  a  wave  length  of  0.4  microns,  0.14  at  1.0  microns, 
and  0.9  at  10.0  microns.  Since  the  absorptivity  of  white  paint  is  so  sensitive  to  spectral  dis¬ 
tribution.  the  reason  for  the  discrepancy  in  the  white  calorimeter  data  is  considered  to  be  due 
to  the  difference  in  spectrum  of  the  calibration  tests  and  thermal  radiation  from  the  nuclear 
fireball.  Since  the  surface  was  covered  with  an  oil  film  from  the  engine  oil  mist  vents  on  the 
underside  of  the  fuselage,  0.30  was  used  in  the  temperature  rise  calculations  for  Fuselage 
Station  389.125. 


TABLE  4.2  SKIN  ABSORPTION  COEFFICIENTS 


Type  of 

Skin 

Absorption  Coefficient 

White  (Clean) 

0.25 

White  (Oil  Film) 

0.30 

Grey  (Clean) 

0.60 

Insignia  Blue 

0.90 

Black 

0.95 

Table  4.3  presents  the  results  of  the  correlation  between  calculated  temperature  rise 
using  measured  radiant  exposure  and  the  measured  temperature  rise.  The  values  in  Table  4.3 
are  fractions  of  points  which  are  categorized  into  magnitude  of  temperature  rise  versus  per¬ 
centage  of  error.  It  shows  that  correlation  is  good  for  temperature  rises  greater  than  30F, 
and  it  is  slightly  improved  for  rises  greater  than  40  F.  The  purpose  of  this  type  of  correlation 
is  to  show  that  data  reading  and  instrumentation  errors  for  low  responses  result  in  errors 
which  do  not  represent  the  validity  of  the  temperature  rise  prediction  methods. 

The  same  method  of  correlation  is  used  to  present  in  Table  4.4  the  comparison  of  calcu¬ 
lated  temperature  rise  using  the  measured  ir radiance  in  a  time-dependent  solution  to  measured 
temperature  rise.  The  same  statement  can  be  made  concerning  good  correlation  for  tempera¬ 
ture  rises  greater  than  30  F.  It  should  be  noted  that  approximately  the  same  percentage  of  data 
paints  fall  within  the  ±  20  percent  error  region  in  both  Tables  4.3  and  4.4.  However,  an  appre¬ 
ciable  increase  in  the  percentage  of  data  points  in  the  ±10  percent  error  region  is  noticed  in 
Table  4.4.  This  concurs  with  the  expectation  that  time -dependent  solutions  result  in  better 
accuracy  than  time -independent  solutions.  The  summaries  presented  In  Tables  4.3  and  4.4 
show  that  both  methods  of  temperature  rise  prediction  for  thin  skins  and  honeycomb  panels 
are  satisfactory  for  low-yield  weapons. 

4.1.4  Time  to  Maximum  Temperature  Rise.  Three  values  of  the  time  to  maximum  tem¬ 
perature  rise  for  selected  thermocouple  locations  were  calculated  for  comparison  with  the 
measured  time  to  maximum  temperature  rise.  These  comparisons  are  presented  in  Table  3.7. 
The  times  to  maximum  temperature  rise  were  calculated  by  using  calculated  time  to  second 
peak  irradiance  (y)  in  the  method  presented  in  A.3.2,  using  measured  y  (Table  3.5)  in  the 
method  presented  in  A.3.2,  and  computing  the  temperature  time  history  from  equation  A.12 
using  measured  irradiance. 

Table  3.7  shows  that  the  largest  discrepancies  exist  between  the  calculated  time  to  maxi¬ 
mum  temperature  rise  using  the  calculated  y  and  the  measured  values.  The  best  correlation 
is  found  between  the  calculated  temperature-time  histories  and  the  measured  data.  This  indi¬ 
cates  that  the  radiometer  lag  is  negligible.  The  comparison  of  calculated  values  using  meas- 
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TABLE  4.3  CORRELATION  BETWEEN  CALCULATED  TEMPERATURE 
RISE  USING  MEASURED  RADIANT  EXPOSURE  AND 
MEASURED  TEMPERATURE  RISE 


Values  are  fractions  of  points  categorised  into  magnitude  of  temperature 
rise  versus  percentage  of  error. 


Percent 

Error 

Region 

Measured 

AT  >  10®F 

Measured 

AT  >  20°F 

Measured 

AT  >  30® F 

Measured 
AT  >  40® F 

±50 

59/63 

49/49 

35/35 

27/27 

(0.94) 

(1.00) 

(1.00) 

(1.00) 

±25 

47/  63 

39/49 

30/35 

22/27 

(0.75) 

(0.80) 

(0.86) 

(0.81) 

±20 

38/63 

32/49 

26/35 

20/27 

(0.60) 

(0.65) 

(0.74) 

(0.74) 

±15 

29/63 

25/49 

21/35 

17/27 

(0.46) 

(0.51) 

(0.60) 

(0.63) 

±10 

16/63 

14/49 

11/35 

9/27 

(0.25) 

(0.29) 

(0.31) 

(0.33) 

TABLE  4.4  CORRELATION  BETWEEN  CALCULATED  TEMPERATURE 
RISE  USING  THE  MEASURED  IRRADIANCE  IN  A  TIME- 
DEPENDENT  SOLUTION  AND  MEASURED  TEMPERATURE 
RISE 

Values  are  fractions  of  points  categorized  into  magnitude  of  temperature 
rise  versus  percentage  of  error. 


Percent 


Error 

Region 

Measured 

AT  >  10®F 

Measured 

AT  >  20®F 

Measured 

AT  >  S0°F 

Measured 
AT  >  40®F 

±50 

48/53 

43/45 

34/34 

21/21 

(0.91) 

(0.96) 

(1.00) 

(1.00) 

±25 

36/53 

32/45 

28/34 

22/27 

(0.68) 

(0.71) 

(0.82) 

(0.81) 

±20 

30/53 

28/45 

25/34 

21/27 

(0.57) 

(0.62) 

(0.74) 

(0.78) 

±15 

22/53 

22/45 

20/34 

18/27 

(0.42) 

(0.49) 

(0.59) 

(0.67) 

±10 

19/53 

19/45 

17/34 

15/27 

(0.36) 

(0.42) 

(0.50) 

(0.56) 

74 


CONFIDENTIAL 


ured  n  and  the  calculated  values  from  temperature-time  histories  shows  good  agreement  for 
unshielded  shots,  while  the  correlation  is  considerably  worse  for  shielded  shots. 

4.1.5  Temperature  at  Shock  Arrival.  One  major  objective  of  the  thermal  analysis  was 
the  development  of  a  time-independent  solution  for  the  prediction  of  temperatures  at  the  time 
of  shock  arrival  for  thin  skins.  Two  values  of  temperature  at  the  time  of  shock  arrival  were 
computed,  using  the  method  presented  in  Appendix  A,  Section  A.3.3  for  comparison  with  the 
measured  values  in  Table  3.8.  One  value  was  calculated  using  calculated  maximum  tempera¬ 
ture,  calculated  time  to  maximum  temperature,  calculated  turbulent  convection  heat  transfer 
coefficient,  and  measured  time-of-shock  arrival.  The  other  value  was  calculated  using  meas¬ 
ured  maximum  temperature,  measured  time  to  maximum  temperature,  calculated  turbulent 
convection  heat  transfer  coefficient,  and  measured  time-to-shock  arrival.  Table  3.8  shows 
that  the  comparison  of  the  latter  calculated  temperature  at  shock  arrival  to  measured  temper¬ 
ature  at  shock  arrival  is  good.  The  prediction  method  is  considered  satisfactory;  however, 
higher  yield  devices  with  thermal  pulses  of  long  duration  will  result  in  large  error.  Therefore, 
Equation  A.  7  should  be  restricted  to  analysis  involving  yields  less  than  100  kt. 

4.1.6  Turbulent  Heat  Transfer  Coefficient.  Using  the  method  presented  in  Appendix  A, 
Section  A.3.4,  Equation  A.9,  values  of  the  flat-plate-turbulent -heat-transfer  coefficient  were 
calculated  from  the  decay  portion  of  the  measured  temperature-time  histories  from  selected 
thermocouples.  The  theoretical  flat  -plate -turbulent -heat -transfer  coefficients  were  also  com¬ 
puted  for  these  thermocouple  locations  by  the  method  presented  in  Appendix  A,  Section  A.3.4, 
Equation  A.8.  These  values  are  compared  for  each  shot  in  Table  3.9.  With  the  exception  of  a 
few  paints,  the  correlation  shows  reasonable  agreement  between  the  two  values.  The  average 
difference  for  all  points  is  30.9  percent,  while  the  mean  value  for  all  points  is  12.5  percent. 

The  theoretical  value  of  the  flat -plate-turbulent-heat-transfer  coefficient  was  therefore  con¬ 
sidered  to  be  a  reasonable  approximation;  therefore,  theoretical  values  were  used  in  the  com¬ 
putation  of  temperature-time  histories,  maximum  temperature  rises,  time  to  maximum  tem¬ 
perature  rises,  and  temperatures  at  the  time  of  shock  arrival. 


4.2  DYNAMIC  RESPONSE 

4.2.1  Measured  Response.  The  measured  data  of  Operation  Plumbbob  which  described 
the  dynamic  response  of  the  airplane  to  the  blast  wave  was  reliable  and  consistent.  The  re¬ 
liability  of  this  data  may  best  be  demonstrated  by  the  fact  that,  of  the  29  recording  channels 
required  for  adequate  analysis  of  the  longitudinal  response  on  seven  flights,  there  were  only 
four  traces  that  were  considered  to  be  unusable.  These  figures  imply  a  reliability  factor  of  98 
percent,  which  is  considered  excellent.  The  consistency  of  both  the  rigid-body  and  structural 
motions,  as  recorded  for  each  shot,  appeared  to  be  good.  The  similarity  of  the  recorded  vari¬ 
ables  as  functions  of  time  after  shock  arrival  may  be  seen  in  Appendix  G.  Because  of  this  con¬ 
sistency  of  measured  data  the  following  discussions  will  deal  primarily  with  the  analysis  of 
Shot  Shasta,  as  typical. 

4.2.2  Calculated  Response.  The  calculated  dynamic -response  data  presented  in  Chapter  3 
were  determined  by  the  methods  discussed  in  Appendix  B.  The  atmospheric  and  flight  condi¬ 
tions  that  determined  the  parametric  values  for  the  analysis  were  those  measured  at  the  time 
of  each  shot  and  presented  in  Tables  3.1,  3.3,  and  3.4.  These  measured  conditions  were  used  In 
an  attempt  to  simulate  as  accurately  as  possible  the  actual  response  characteristics  of  the 
airplane  during  each  event.  The  blast-wave  parameters  that  were  used  in  the  response  analysis 
are  listed  in  Table  3.10.  These  parameters  were  computed  from  the  final  measured  yields  of 
Table  3.2  by  the  methods  described  In  Section  B.l.  It  was  felt  that  more  reliable  data  would  be 
obtained  from  an  analysis  based  on  the  measured  yield  rather  than  one  which  utilised  the  peak 
overpressure  behind  the  shock  front  as  measured  by  the  test  vehicle,  because  of  the  relative 
inconsistencies  in  the  peak  values  recorded  by  the  transducers  for  each  shot. 
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As  Indicated  in  Chapter  S.  two  methods  of  analysis  were  used  to  calculate  the  dynamic  re¬ 
sponse.  These  methods  are  referred  to  as  the  original  equations  and  the  modified  equations. 

A  thorough  discussion  of  the  development  of  each  method  is  presented  In  Section  B.2.  it  Is 
sufficient  to  say  that  the  original  equations  were  developed  before  the  test  operations  and  that 
the  modified  equations  were  formulated  after  the  testing.  The  modified  equations  Incorporate 
both  simplifications  and  refinements  that  facilitate  a  more  accurate  determination  of  the  air¬ 
plane’s  response  characteristics. 

4.1.3  Correlation  of  Maximum  Responses.  The  calculated  and  measured  maximum  values 
of  the  airplane  stability  variables  after  shock  arrival  are  presented  in  Table  3.11  for  all  shots. 
The  consistency  of  relative  magnitudes  between  the  calculated  and  measured  values  for  each 
shot  is  apparent.  Perturbation  of  angle  of  attack,  which  is  a  direct  indication  of  the  blast-wave 
gust  effect  as  well  as  the  incremental  normal  load  factor  and  perturbation  of  angle  of  pitch, 
which  indicate  the  amplitude  of  rigid-body  response,  should  be  associated  directly  with  the 
blast  incidence  angle  and  calculated  peak  overpressure  of  Tables  3.1  and  3.10,  respectively. 
This  comparison  Indicates  that  the  blast -wave  gust  effect  and  rigid-body  response  are  both 
very  nearly  proportional  to  peak  overpressure  and  blast  angle.  It  should  be  noted  that  the  peak 
overpressures  considered  are  relatively  small  and  that  larger  values  of  overpressure  for  low 
blast  incidence  angles  would  Introduce  a  significant  nonlinear  effect.  This  nonlinearity  is  due 
to  an  alleviation  of  wing  load  by  the  reduction  of  local  dynamic  pressure.  It  Is  the  horizontal 
component  of  material  velocity  behind  the  shock  front  which  reduces  the  resultant  airspeed  and 
consequently  lowers  the  dynamic  pressure.  This  load  alleviation  at  low  incidence  angles  is  of 
extreme  Importance  to  the  establishment  of  the  airplane’s  delivery  capabilities. 

Mavimiun  shear  and  bending-moment  stress  levels  as  measured  In  the  wing  during  the 
test  operation  are  listed  in  Table  3.12  with  the  corresponding  values  from  the  theoretical 
analysis.  Although  the  consistency  of  relative  magnitudes  between  the  calculated  and  measured 
values  for  each  shot  is  apparent,  there  are  large  discrepancies  between  the  computed  and 
measured  absolute  magnitudes  of  stress,  with  but  one  exception.  The  exceptional  case  is  that 
of  the  most  critical  structural  component  at  Wing  Station  17.5.  The  stress  level  at  this  loca¬ 
tion  is  due  to  bending-moment  in  the  wing  carry-through  structure,  which  is  Initially  a  direct 
result  of  the  total  wing  loading.  Since  the  total  wing  loading  directly  determines  this  stress  by 
Its  integrated  value  of  bending-  moment  and  since  the  straightforward  design  of  the  carry- 
through  structure  lends  Itself  easily  to  theoretical  stress  analysis,  the  stress  at  Wing  Station 
17.5  may  be  computed  with  a  high  degree  of  accuracy.  The  peak-stress  level  at  this  location 
occurs  at  approximately  0.1  second  after  shock  arrival,  which  is  coincident  with  maximum 
normal  load  factor. 

The  bending-moment  stress  was  measured  at  Wing  Station  17.5  on  both  the  right  and  left 
wings.  This  duplication  of  gages  was  made  to  insure  the  accuracy  of  this  most  critical  stress. 
From  Table  3.12  it  Is  seen  that  the  average  variation  of  measured  maximum  percent  of  limit- 
allowable  stress  between  the  two  wings  was  only  1.26  percent  for  all  shots,  and  that  the  maxi¬ 
mum  variation  was  2.12  percent  for  Shot  Hood.  This  excellent  correlation  proves  the  validity 
of  measured  stress  levels  at  this  structural  point. 

Comparison  of  the  stresses  calculated  by  the  modified  equations  to  the  respective  meas¬ 
ured  values  at  Wing  Station  17.5  demonstrates  the  high  degree  of  capability  of  the  theoretical 
analysis  to  predict  this  most  critical  stress  in  the  FJ-4  airplane  during  blast  Impingement.  It 
is  semi  that  the  average  variation  of  the  maximum  percent  of  limit-allowable  stress  between 
the  calculated  and  measured  values  was  only  2.55  percent,  based  on  the  results  of  all  shots. 
The  largest  corresponding  error  was  6.57  percent  of  limit -allowable  stress  for  the  left  wing 
in  Shot  Diablo. 

The  calculated  values  of  both  shear  and  bending-moment  stresses  outboard  of  Wing  Station 
17.5  show  poor  correlation  with  their  respective  measured  values.  These  Inconsistencies  have 
been  tentatively  attributed  to  the  extreme  difficulties  encountered  in  the  determination  of  the 
stress  distribution  from  the  modal  deflection  shapes  for  such  a  complex  structure.  These 
stress  distributions  are  necessary  to  establish  the  modal  stresses  which  are  used  in  conjunc¬ 
tion  with  the  mode-displacement  method  discussed  In  Appendix  B.  The  modal  stresses  were 
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used  In  a  direct  determination  of  stress  response  from  the  mode  deflections  of  the  theoretical 
response  analysis.  It  is  now  clear  that  experimental  evaluation  of  modal  stresses  for  such 
complex  structures  as  the  FJ-4  wing  should  be  made  during  the  ground -vibration  tests  of  the 
airframe. 

Maximum  bending-moment  stress  levels,  as  measured  In  the  horizontal  stabiliser  during 
the  test  operation,  are  presented  In  Table  3.13  with  the  corresponding  values  of  limit -allowable 
stress.  It  is  noted  that  the  maximum  stabiliser  stresses  recorded  during  the  tests  are  much 
less  critical  than  those  of  the  wing.  This  characteristic  may  be  attributed  to  two  separate  fac¬ 
tors.  The  first  factor  is  that  pitching  motion  of  the  airplane  has  a  large  aerodynamic  alleviat¬ 
ing  effect  on  the  tall  surfaces.  As  seen  by  the  values  of  pitch  rate  In  Table  3.11,  pitching 
motion  is  evident  characteristic  of  the  airplane  response,  and  as  a  consequence  this  alleviat¬ 
ing  effect  is  significant.  The  second  factor  is  the  difference  of  the  design  requirements  for  the 
two  surfaces.  The  FJ-4  airplane  has  a  horizontal  stabilizer  which  was  designed  on  the  basis  of 
combined  strength  and  flutter  requirements,  whereas  the  wing  was  designed  primarily  on  the 
basis  of  strength  requirements  alone.  As  a  result,  the  tall  possesses  basically  higher  strength 
than  the  wing:  accordingly,  the  percent  of  allowable  stress  levels  under  similar  loadings  tends 
to  be  lower. 

4.3.4  Rigid  Airplane  Response.  As  the  wave  front  rapidly  envelopes  the  airplane,  the  lift¬ 
ing  surfaces  experience  an  increased  aerodynamic  loading  because  of  the  material  velocity 
behind  the  shock  front.  These  transient  aerodynamic  loads  result  in  a  combined  heaving  and 
pitching  motion  of  the  airplane.  These  motions  may  be  characterized  by  an  increased  normal 
load  factor  and  a  build-up  of  negative  pitch  rate.  After  these  parameters  have  reached  their 
maximum  values  in  responding  to  the  gust-like  blast  wave,  the  airplane  experiences  a  highly- 
damped  short -period  oscillation  about  the  trimmed  flight  conditions  existing  prior  to  shock 
arrival. 

The  variation  of  normal  load  factor  as  illustrated  in  Figure  3.17  for  Shot  Shasta  was  de¬ 
termined  from  normal  acceleration  recorded  by  a  linear  accelerometer  located  in  the  fuselage 
structure  near  the  airplane  center  of  gravity.  In  this  location,  the  recorded  acceleration  signal 
contained  accelerations  due  to  structural  vibrations  as  well  as  those  due  to  the  total  airplane 
normal  loading.  This  effect  is  evident  in  the  measured  normal  load  factor  of  Figure  3.17. 
During  the  first  0.1  second  following  arrival  of  the  shock  front,  local  high-frequency  vibrations 
at  the  fuselage  structure  may  be  seen.  After  this  effect  has  been  sufficiently  damped,  excita¬ 
tion  at  the  first  vibration  mode  of  the  wing  becomes  apparent.  As  a  result  of  these  structural 
effects,  it  was  difficult  to  correlate  the  maximum  measured  normal  acceleration  with  the  cor¬ 
responding  value  obtained  from  the  dynamic -response  analysis.  The  superposition  of  the 
structural  vibrations  on  the  normal-acceleration  signal  may  be  reduced  by  installing  the  ac¬ 
celerometer  at  the  intersection  of  the  first  two  wing -mode  node  lines.  The  negative  portion  of 
the  rigid-airplane  response,  described  by  the  normal  load  factor,  compares  well  with  the  cor¬ 
responding  values  calculated  from  the  modified  equations.  This  noted  improvement  over  the 
results  of  the  original  equations  is  primarily  because  of  the  misrepresentation  of  the  airplane’s 
pitching  characteristics  in  the  original  analysis. 

The  pitching  motion  of  the  airplane  because  of  the  induced  loads  of  the  blast  wave  is  Il¬ 
lustrated  in  Figures  3.18  and  3.19.  Predictions  of  the  maximum  negative  values  of  both  the 
pitch  rate  and  perturbation  of  angle  of  pitch  by  the  modified  equations  are  in  good  agreement 
with  the  respective  measured  amplitudes.  Although  these  maximum  values  are  in  agreement, 
the  basic  frequency  of  the  response  appears  to  be  in  error.  This  pitching  frequency  during  the 
early  portion  of  the  response  is  directly  dependent  upon  the  nature  of  the  material  velocity 
behind  the  shock  front  and  the  response  characteristics  of  the  airplane.  The  representation  of 
the  airplane’s  basic  pitching-response  characteristics  have  been  verified  by  correlation  with 
flight  test  data  of  prior  studies.  It  has  therefore  been  concluded  that  the  error  In  the  frequency 
of  response  is  due  to  the  misrepresentation  of  the  material  velocity  as  a  function  of  time. 

The  theoretical  analysis  utilized  a  material -velocity  function  similar  to  the  exponential 
function  of  static  overpressure  behind  the  shock  front.  B  Is  presently  felt  that  the  actual 
material-velocity  profile  has  a  longer  positive  phase  and  a  much  faster  initial  rate  of  decay. 
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If  tills  is  the  esse,  then  the  characteristic  frequency  of  the  calculated  pitching  response  would 
be  reduced  and  correlation  with  the  measured  data  would  be  greatly  Improved. 

The  relatively  poor  similarity  between  the  pitching  response  of  the  original  analysis  and 
the  measured  data  Is  largely  attributed  to  two  misrepresentations  in  the  aerodynamic  loading. 
The  first  discrepancy  was  found  to  be  an  error  in  the  analog  mechanisation  which  reduced  the 
air  loads  on  the  horlaontal  stabiliser  and  consequently  lowered  the  short-period  pitching  fre¬ 
quency  of  the  airplane.  The  second  discrepancy  was  the  effective  center-of -pressure  location 
for  the  blast  loading  on  the  wing.  This  center  of  pressure  was  assumed  to  be  too  far  aft,  and 
as  a  result  the  blast -induced  negative  pitching  moment  was  greatly  amplified.  These  two  mis¬ 
representations  In  the  original  equations  made  the  calculated  pitching  motions  too  violent  and 
of  too  low  a  frequency  as  seen  in  Figures  3.18  and  3.19. 

The  airplane  angle  of  attack  as  Illustrated  In  Figure  3.20  was  measured  by  free -stream 
vanes  mounted  on  the  end  of  a  relatively  long  nose  boom.  The  recorded  signal  consisted  of 
both  the  airplane  angle  of  attack  and  induced  angle  of  attack  of  the  forward  end  of  the  nose 
boom.  This  induced  angle  of  attack  was  generated  by  structural  vibrations  of  the  boom  which 
were  excited  by  the  sharp-edged  blast  wave.  Airplane  angle  of  attack  was  measured  as  a 
means  of  determining  the  normal-gust  effect  of  the  blast  wave.  A  mean  initial  value  of  the 
oscillatory  response  shown  in  Figure  3.20  was  used  to  indicate  this  effect  and  is  compared  in 
Table  3.11  to  the  corresponding  values  utilized  in  the  two  theoretical  analyses. 

The  positions  of  all  the  longitudinal  control  surfaces  were  recorded  on  each  flight  to  detect 
airy  extraneous  inputs  due  to  pilot  reaction.  The  only  measurable  control-surface  displace¬ 
ments  were  found  to  be  of  the  horizontal  stabilizer.  As  illustrated  in  Figure  3.21  for  Shot 
Shasta,  the  stabilizer  motions  were  quite  small.  Although  the  aerodynamic  loads  because  of 
these  small  displacements  of  the  horizontal  stabilizer  were  nearly  negligible,  they  were  in¬ 
cluded  in  the  dynamic  response  analysis  in  an  attempt  to  simulate  as  accurately  as  possible 
the  response  characteristics  of  the  airplane  during  each  test. 

4.2.5  Wing  Structural  Response.  As  the  blast  wave  sweeps  across  the  wing,  a  substantial 
increase  in  static  and  dynamic  pressure  Is  experienced  by  the  underside  of  the  lifting  surface. 
This  abrupt  increased  loading  tends  to  excite  the  normal  vibration  modes  of  the  wing.  These 
normal  modes  are  seen  as  complex  deformations  of  the  structure,  vibrating  at  relatively  high 
frequencies.  These  motions  served  as  the  basic  means  of  evaluating  the  stress  levels  existing 
throughout  the  structure  fallowing  blast  impingement. 

The  most  critical  stress  experienced  by  the  wing  was  at  Wing  Station  17.5  of  the  carry- 
through  structure.  Failure  at  this  location  would  occur  as  compression  in  the  upper  member 
of  the  box-like  structure  due  to  the  accumulated  bending  moment  from  the  normal  dynamic 
wing  loading.  The  variation  of  bending-moment  stress  during  Shot  Shasta  for  both  the  right  and 
left  wings  at  Wing  Station  17.5  is  illustrated  in  Figures  3.22  and  3.23,  respectively.  Also  shown 
are  the  corresponding  perturbation  stress  levels  as  calculated  by  both  the  original  and  modified 
equations  of  motion. 

The  bending-moment  stress  at  this  most  critical  structural  point  is  largely  due  to  excita¬ 
tion  of  the  first  normal  mode  of  vibration  of  the  wing,  as  seen  in  the  Illustration  by  the  charac¬ 
teristic  frequency  of  approximately  6  cps.  It  should  be  noted  that  this  relatively  high-frequency 
stress  variation  is  superimposed  on  the  perturbation  stress  level  induced  by  the  rigid-body 
motions  of  the  airplane  and  that  the  maximum  combination  of  stress  occurs  at  approximately 
0.1  second  following  arrival  of  the  shock  front.  At  this  time,  both  contributions  to  the  stress 
level  reach  their  peak  values,  which  results  In  an  overstress  or  dynamic -magnification  factor 
of  approximately  1.8.  This  relatively  large  dynamic  -magnification  factor  is  due  to  the  discon¬ 
tinuous  nature  of  the  material  velocity  at  the  shock  front,  which  Is  Instrumental  In  exciting  the 
structural  modes  of  vibration. 

The  corresponding  structural  response  determined  by  the  modified  equations  agrees 
closely  with  the  measured  data.  The  only  deviation  of  the  two  response  curves  is  due  to  the 
stresses  induced  by  the  airplanes  pitching  motions,  and  this  error  is  not  evident  until  approxi¬ 
mately  0.5  second  after  maximum  stress  occurs.  The  response  calculated  from  the  original 
equations  also  predicts  a  valid  maximum  stress  level,  although  the  deviation  from  the  meas- 
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ured  curve  which  occurs  after  this  level  has  been  reached  is  quite  sltable.  The  apparent  inac¬ 
curacies  of  both  calculated  responses  are  primarily  due  to  the  misrepresentations  of  the  air¬ 
plane  motions  in  pitch,  which  have  previously  been  discussed. 

Other  wing-bending  moment  stresses  that  were  recorded  during  Shot  Shasta  and  presented 
in  Chapter  3  were  at  Wing  Station  36.5  in  the  forward  spar  and  at  Wing  Station  178  in  the  root  of 
the  outer  wing  panel.  Figure  3.24  illustrates  fairly  good  agreement  between  the  measured  and 
calculated  structural  responses  of  Wing  Station  36.5.  This  bending-moment  stress  Is  primarily 
due  to  the  differential  bending  moment  carried  by  the  forward  torque  box  In  reaction  to  the 
aerodynamic  torque  acting  on  the  lifting  surface.  The  calculated  bending-moment  stress  at 
Wing  Station  173  of  Figure  3.25  demonstrates  fairly  poor  agreement  with  the  corresponding 
value  of  measured  stress.  This  disagreement  has  been  attributed  to  the  misrepresentation  of 
the  airplane  pitching  motions  and  the  excess  contribution  of  the  second  vibration  mode.  Both 
these  effects  have  been  previously  discussed  as  erroneous  airplane  pitching  recovery  and  Im¬ 
proper  evaluation  of  the  modal  stresses. 

The  wing  shear  stresses  as  illustrated  in  Figures  3.26  to  3.33  for  Shot  Shasta  indicate 
clearly  the  erroneous  contributions  of  the  higher  vibration  modes  to  the  calculated  shear  dis¬ 
tribution  in  the  wing.  As  discussed  in  Sect! cm  4.2.3,  this  poor  correlation  of  the  calculated 
shear  stresses  with  the  measured  data  has  been  attributed  to  the  improper  evaluation  of  the 
modal  stresses.  It  is  suggested  that  the  mode -displacement  method,  which  was  used  to  calcu¬ 
late  these  values,  be  discarded.  This  eliminates  the  requirement  of  having  to  perform  a  de¬ 
tailed  calculation  of  the  inertial-load  distribution  and  the  subsequent  static-stress  analysis  of 
this  complex  structure  for  each  normal  mode  utilized  by  the  analysis.  Instead,  it  is  proposed 
that  a  more  direct  approach  be  taken  to  the  evaluation  of  the  modal  stresses.  During  the  con¬ 
ventional  ground-vibration  tests  of  the  aircraft,  all  strain-gage  outputs  can  be  recorded  and 
correlated  with  high-speed  movies  of  the  wing-tip  deflections  as  each  normal  mode  is  excited. 
This  correlation  will  determine  the  modal  stresses  directly  with  reduced  cost,  time  and  man¬ 
power.  The  accuracy  of  the  modal  stresses  will  then  be  within  the  tolerances  of  the  recording 
equipment,  which  is  sufficient  to  fully  evaluate  the  methods  of  dynamic-response  analysis  used 
in  this  operation. 

4.2.6  Factors  Affecting  the  Delivery  Capability.  For  an  airplane  to  deliver  its  weapon 
with  the  maximum  of  accuracy  over  a  target,  it  is  a  generally  accepted  fact  that  the  carrier 
must  be  as  close  over  the  point  of  impact  as  possible  at  the  time  of  release.  This  introduces 
the  problem  of  the  airplane  having  to  experience  the  violent  shock  wave  generated  by  the 
weapon  during  its  escape  maneuver.  Mentioned  in  Section  4.2.3  were  two  important  character¬ 
istics  of  the  dynamic  response  which  have  a  direct  bearing  on  the  choice  of  the  optimum  mode 
of  delivery.  The  first  of  these  factors  is  that  the  magnitude  of  response  decreases  with  the 
blast  incidence  angle  as  a  result  of  the  reduced  normal  component  of  material  velocity  behind 
the  shock  front.  And  secondly,  there  is  a  definite  alleviation  of  aerodynamic  loading  because 
of  the  increased  horizontal  component  of  material  velocity  at  low  blast  angles.  Both  these  ef¬ 
fects  indicate  that  a  low-angle-escape  maneuver  will  enable  the  airplane  to  be  as  close  to  the 
point  of  impact  as  is  practicable  at  the  time  of  release. 

To  fully  evaluate  the  delivery  capability  of  the  FJ-4/4B  airplane,  the  validity  of  die 
dynamic -response  analysis  should  be  proven  at  conditions  of  low  blast  angle  and  high  over¬ 
pressure.  These  conditions  were  not  achieved  during  Operation  Plumbbob  because  of  the  com¬ 
bined  restrictions  of  the  low  yields  detonated  and  geographical  features  of  the  test  site.  It  is 
strongly  recommended  that  further  flight  tests  be  conducted  under  these  optimum  escape  con¬ 
ditions  to  establish  the  validity  of  the  dynamic -response  analysis  and  subsequently  facilitate 
the  establishment  of  the  maximum  delivery  capabilities  of  both  the  FJ-4/4B  and  future  de¬ 
livery  systems. 

4.3  WING  PRESSURE  SURVEY 

A  wing  static -overpressure  survey  of  the  FJ-4  was  made  on  five  flights  of  the  operation. 
The  objective  of  this  investigation  was  to  obtain  a  better  knowledge  of  the  nature  of  the  Induced 
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aerodynamic  loading  on  a  lifting  surface  during  impingement  of  a  high -Intensity  blast  wave. 
Static  overpressure  was  measured  at  14  points  on  the  upper  and  lower  surfaces  of  the  chord 
at  Wing  Station  175.75.  These  measurements  were  made  by  pressure  transducers  installed  at 
points  along  the  surfaces  in  a  manner  that  would  facilitate  an  adequate  description  of  the  upper 
and  lower  chord-wise  pressure  distributions.  The  recorded  fluctuation  of  static  pressure  by 
each  of  these  gages  is  illustrated  in  Figures  3.35  and  3.36  for  Shot  Shasta. 

These  time  histories  of  static  overpressure  were  cross  plotted  with  percent  chord  at  con¬ 
secutive  intervals  in  time  following  arrival  of  the  shock  front.  These  chord-wise  distributions 
are  typified  by  those  of  Shot  Shasta  which  are  presented  in  Figures  3.3?  to  3.39.  Since  the  air¬ 
plane  was  positioned  at  a  relatively  low  blast-incidence  angle  for  Shot  Shasta,  the  chord-wise 
traversal  of  the  shock  front  is  indicated  by  the  progressively  increased  pressure  along  the 
lower  surface  of  the  wing  during  the  first  9  msec  of  blast-wave  impingement.  The  pressure- 
distribution  recovery  from  this  transient  effect  of  the  shock  front  is  noted  as  an  Increased 
loading  along  the  aft  portion  of  the  upper  surface.  A  quasi-steady  condition,  which  is  similar 
to  that  of  an  increased  angle  erf  attack  at  an  elevated  ambient  pressure,  is  then  assumed  by  the 
pressure  distribution.  This  condition  is  a  natural  consequence  of  the  quasi-steady  nature  of 
the  local  overpressure  and  material  velocity  within  the  blast  wave. 

The  net  values  of  these  chord-wise  distributions  of  static  overpressure  were  integrated  to 
determine  the  nature  of  the  section  lift  and  center  of  pressure  corresponding  to  the  blast- 
induced  loading.  The  integrated  values  for  Shot  Shasta  are  illustrated  in  Figures  3.40  and  3.41. 
A  notable  effect  is  the  peak  in  section  lift  occurring  at  approximately  8  msec  after  arrival  of 
the  shock  front  at  the  trailing  edge.  This  peak  load  is  a  direct  result  of  the  progressively  In¬ 
creasing  pressure  on  the  underside  of  the  wing,  as  indicated  in  the  preceding  paragraph.  The 
effect  has  been  tentatively  attributed  to  the  lagging  of  the  shock  wave  on  the  upper  surface  be¬ 
hind  its  lower  counterpart  because  of  its  diffraction  and  subsequent  expansion  about  the  trailing 
edge.  This  phenomenon  will  be  referred  to  as  the  diffraction  pulse.  Associated  with  this  peak 
loading  is  a  transient  negative  moment  due  to  the  travel  of  the  center  of  pressure  from  the 
trailing  edge  forward.  This  moment  becomes  positive  for  approximately  2  msec  as  the  in¬ 
duced  loading  attains  a  minimum  value  before  recovery  to  its  quasi-steady  condition  referred 
to  above. 

An  attempt  was  made  to  determine  the  effect  of  the  diffraction  pulse  on  the  structural  re¬ 
sponse  of  the  wing.  A  pulse  generator  was  incorporated  in  the  analog  mechanization  of  the 
dynamic-response  equations  to  simulate  the  measured  loading.  The  generated  pulse  is  shown 
as  an  addition  to  the  lift -lag  function  of  the  basic  analysis  in  Figure  3.42  for  Shot  Hood.  For 
comparative  purposes  the  measured  section  lift  is  also  presented.  The  diffraction  pulse  was 
essentially  described  by  two  basic  parameters.  These  parameters  are  pictorially  defined  in 
Figure  3.42  as  amplitude  ratio  and  time-to-peak  diffraction  load.  The  calculated  diffraction 
effect  on  the  maximum  value  of  percent  limit -allowable  stress  at  Wing  Station  17.5  is  illus¬ 
trated  in  Figure  3.43  as  a  function  of  these  parameters.  It  can  be  seen  that  the  diffraction 
pulse  has  a  negligible  effect  on  this  critically  stressed  member,  even  though  the  airplane  was 
positioned  at  a  relatively  high  angle  where  the  effect  may  be  expected  to  be  the  greatest.  Vari¬ 
ation  of  the  descriptive  parameters  demonstrated  that  the  diffraction  pulse  can  have  significant 
effects  on  the  structural  response.  As  a  consequence,  there  is  a  definite  requirement  for  a 
more  thorough  fundamental  knowledge  of  the  diffraction  phenomenon,  since  it  may  well  be  In¬ 
strumental  in  the  establishment  of  the  delivery  capabilities  of  future  aircraft. 

4.4  NUCLEAR  RADIATION 

There  was  considerable  variation  in  the  total  nuclear  dose  recorded  at  the  various  posi¬ 
tions  in  the  aircraft  (Table  3.14).  Maximum  values  were  consistently  obtained  at  the  ammuni¬ 
tion  bay  location  because  of  the  lesser  amount  of  shielding  from  the  aircraft  structures  com¬ 
pared  to  other  positions.  Although  the  limiting  positioning  criterion  for  nuclear  dose  was  five 
rem  for  each  shot,  the  inherent  safety  factors  in  the  positioning  procedures  precluded  obtaining 
values  this  large.  The  post  shot -calculated  values  of  total  nuclear  dose  shown  in  Table  3.14 
were  computed  by  the  methods  described  in  Appendix  C,  using  final  yields. 
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For  the  heavily  shielded  shots  the  positioning  criteria  included  the  reduction  of  the  total 
nuclear  dose  by  a  factor  of  1,000.  In  general,  the  correlation  between  postshot-calculated 
gamma  dose  and  the  test  measurements  was  fair;  however,  the  data  Indicates  that  for  heavily 
shielded  shots  the  reduction  of  the  calculated  dose  by  a  factor  of  1,000  would  reault  in  a  con¬ 
siderable  underprediction.  Hence,  a  reduction  by  a  factor  of  10  applied  to  the  heavily  shielded 
shots  was  made  and  can  be  noted  in  Table  3.14  as  giving  fair  correlation. 

For  positioning  purposes,  1  r  was  added  to  the  gamma  radiation  dose  for  an  aircraft  being 
more  than  1,000  feet  short  of  and  approaching  ground  aero  at  time  aero.  It  was  found  that  using 
this  procedure  in  the  final  analysis  improved  the  data  correlation  shown  in  Table  3.14. 

4.5  FREE  STREAM  OVERPRESSURE 

Peak  static  overpressure  data  were  obtained  from  transducers  located  in  the  nose  boom, 
fin  top,  and  right  and  left  wing  tips.  Because  of  erratic  behavior  of  the  left  wing  tip  transducer 
it  was  relocated  in  the  left  side  of  the  fuselage  after  Shot  Diablo.  It  was  also  noted,  after  com¬ 
pletion  of  field  participation,  that  the  oscillograph  traces  from  the  fin  top  transducer  for  all 
shots  and  the  traces  from  the  right  wing  tip  transducer  for  Shots  Doppler  and  Smoky  were  so 
erratic  that  even  the  peak  readings  were  considered  to  be  unreliable.  The  tabulated  data  in 
Table  3.15  therefore  does  not  include  the  readings  for  the  fin  top  transducer,  and  the  average 
measured  overpressures  reported  therein  do  not  include  these  and  other  erratic  data  In  their 
calculation. 

Overpressures  were  calculated  using  the  data  listed  in  Tables  3.1  through  3.3  and  the 
method  outlined  in  Appendix  A.  A  comparison  of  these  calculated  values  with  the  average 
measured  peak  static  overpressures  shows  excellent  correlation,  with  a  maximum  variation 
of  0.08  psi  and  an  average  variation  of  0.036  psi.  There  was  also  no  significant  bias  toward 
either  overprediction  or  underprediction.  A  graphical  comparison  (Figure  3.44)  of  the  meas¬ 
ured  data,  scaled  to  1  kt,  homogeneous  sea-level  atmospheric  conditions,  with  the  basic 
curve  obtained  from  Reference  5  further  emphasizes  the  excellent  correlation  in  the  low 
overpressure  region  for  which  data  was  obtained. 

4.6  TIME  OF  SHOCK  ARRIVAL 

Time  erf  shock  arrival  values  were  obtained  from  the  various  pressure  data  from  all  four 
oscillographs.  Correlation  of  these  sets  of  data  was  excellent,  with  a  maximum  deviation  of 
0.04  second  for  all  shots.  The  average  deviation  was  approximately  0.015  second.  The  meas¬ 
ured  values  are  listed  in  Table  3.8.  Multiple  shocks  were  recorded  for  Shots  Hood,  Diablo, 
Kepler,  and  Doppler.  In  particular,  three  shocks  were  recorded  for  Shots  Hood  and  Doppler  and 
two  shocks  were  observed  during  Shots  Diablo  and  Kepler.  The  responses  resulting  from  the 
second  and  third  shock  waves  were  observed  to  be  insignificant. 

In  Figure  3.45,  the  measured  values  of  the  time  of  arrival  of  the  first  shock,  scaled  ac¬ 
cording  to  the  methods  described  in  Appendix  B,  are  compared  with  the  1  kt,  homogeneous, 
sea-level  atmosphere  curve  from  Reference  1.  The  excellent  correlation  demonstrates  the 
accuracy  of  existing  methods  in  predicting  times  of  shock  arrival  for  the  low  yield  range  ex¬ 
perienced  in  Operation  Plumbbob. 


4.7  ENGINE  RESPONSE 

Engine  response  parameters  were  recorded  during  all  of  the  seven  shots  in  which  the 
FJ-4  aircraft  participated.  These  parameters  consisted  of  compressor  inlet  temperature  and 
pressure,  compressor  discharge  pressure,  tailpipe  total  pressure  and  temperature,  fuel  flow, 
and  engine  speed.  Time  histories  of  these  parameters  were  presented  in  Chapter  3  for  Shot 
Boltzmann.  Additional  time  histories  are  presented  in  Appendix  I  for  the  other  six  shots  in 
which  FJ-4  aircraft  participated.  The  data  from  Shot  Boltzmann  are  representative  and  suffi¬ 
cient  for  a  study  of  engine  response. 
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In  analysing  the  engine  operating  parameters  during  the  Boltzmann  participation  the  engine 
eras  noted  to  operate  at  essentially  steady  state  conditions  between  time  zero  and  time  of  shock 
arrival.  The  thermal  pulse  produced  a  slight  Increase  in  compressor  Inlet  temperatures  and 
exhaust  gas  temperatures,  while  engine  pressure  parameters  showed  minor  disturbances. 

Fuel  flow  variations  were  In  the  order  of  300  pounds  per  hour;  however,  this  amount  of  varia¬ 
tion  would  be  considered  normal  over  the  time  span  noted.  This  Is  based  on  the  fact  that  fuel 
control  bias  response  ratios  are  In  the  order  of  5  to  8  seconds  for  the  Inlet  temperature  sys¬ 
tem,  and  1  to  3  seconds  for  the  pressure  bias  system.  Therefore,  it  Is  believed  that  fuel  flow 
variations  noted  were  mainly  due  to  fuel  control  action  required  to  maintain  constant  speed. 

At  the  time  of  shock  arrival  sharp  increases  In  all  engine  parameters  were  noted.  In 
analyzing  the  compressor  inlet  pressure  data,  it  was  noted  that  the  net  effect  of  the  shock  wave 
was  a  general  increase  in  ambient  pressure.  While  some  variation  was  noted,  the  pressure  did 
not  change  in  an  oscillatory  manner.  Compressor  inlet  temperature  Increased  2F  while  the 
engine  exhaust  gas  temperature  increased  7F.  The  latter  was  attributed  to  the  2F  Inlet  tem¬ 
perature  increase  plus  an  increase  in  engine  speed  of  33  rpm.  Engine  fuel  flow  fluctuated 
from  plus  to  minus  450  pounds  per  hour  at  1  cps.  However,  this  flow  oscillation  rapidly  con¬ 
verged  and  was  completely  damped  in  4.5  seconds.  The  compressor  discharge  pressure  showed 
the  same  characteristics  as  the  inlet  pressure  and  was  in  phase,  indicating  that  compressor 
pressure  ratio  was  essentially  constant. 

In  comparing  the  engine  speed  and  fuel  flow  data  from  Shot  Boltzmann  to  the  remaining  six 
shots,  the  same  general  characteristics  were  noted.  However,  data  obtained  during  Shots  Hood 
and  Doppler  from  time  zero  to  time  of  shock  arrival,  and  during  Shot  Diablo  after  time  of 
shock  arrival,  show  abnormal  variations  in  engine  speed  and  fuel  flow.  It  has  been  confirmed 
that  during  the  participation  in  Shots  Hood  and  Doppler  the  variations  noted  were  due  to  pilot 
maneuver  of  the  power  lever.  This  was  also  believed  to  be  the  case  during  the  Shot  Diablo 
participation  when  a  fuel  flow  fluctuation  of  4,000  pounds  per  hour  was  noted  after  time  of 
shock  arrival. 

In  summarizing  the  data,  it  was  evident  that  the  engine  was  sensitive  mainly  to  the  shock 
wave  and  the  associated  increase  in  ambient  air  pressure;  variations  noted  were  not  sufficient 
to  cause  engine  damage.  The  fuel  flow  fluctuations  following  time  of  shock  arrival  were 
abnormal  for  steady  state  operation  and  if  continued  at  ±490  pounds  per  hour  about  the  steady 
state  flow  at  a  higher  frequency  in  an  undamped  maneuver  would  probably  cause  flame  out. 
However,  dampening  and  slow  variation  were  noted  in  the  data  for  all  flights.  It  is  significant 
that  pressure  variations  took  the  effect  of  a  general  increase  and  were  not  oscillatory.  It  is 
evident  from  this  behavior  that  the  pressure  variations  noted  will  not  cause  divergent  engine 
instability.  Power  plant  structural  damage  will  not  be  caused  by  these  pressure  variations 
since  maximum  operating  case  differentials  in  the  order  at  100  psi  are  normally  experienced 
over  the  operating  range  of  the  engine. 
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Chapter  5 

CONCLUSIONS  and  RECOMMENDATIONS 

5.1  CONCLUSIONS 

Tlie  following  conclusions  can  be  drawn  from  the  final  analyses  of  the  test  data. 

1.  The  radiant  exposure  data  correlation  indicates  that  the  thermal  output  of  the  nuclear 
device  is  affected  by  the  degree  of  shielding. 

2.  The  maximum  temperature  rise  prediction  methods  for  both  thin  skins  and  honeycomb 
panels  are  satisfactory. 

3.  The  method  of  predicting  temperatures  at  shock  arrival  for  thin  skins  exposed  to  nu¬ 
clear  detonations  of  less  than  100  kt  is  generally  acceptable. 

4.  The  most  critical  stress  level  in  the  airframe  structure  during  blast  wave  impingement 
occurs  as  bending-moment  stress  at  Wing  Station  17.5  of  the  carry-through  structure. 

5.  Maximum  stresses  of  the  FJ-4  horizontal  stabilizer  were  much  less  critical  than  those 
of  file  wing  following  arrival  of  the  shock  front. 

6.  The  maximum  amplitude  of  airplane  rigid-body  pitching  motion  was  predicted  with  good 
accuracy. 

7.  For  the  relatively  high  blast  angle  positions  of  the  airplane  In  this  operation,  the  blast 
effects  and  rigid-body  response  were  both  nearly  proportional  to  the  free  stream  overpressure 
and  blast  incidence  angle. 

8.  The  theoretical  analysis  indicates  a  nonlinear  alleviation  of  aerodynamic  loading  at  low 
blast  incidence  angles  which  is  of  importance  in  establishing  the  capabilities  of  the  airplane  to 
withstand  the  effects  of  intense  blasts. 

9.  There  was  a  discernible  diffraction  pulse  of  force  on  the  FJ-4  wing  during  blast-wave 
impingement;  however,  this  pulse  bad  a  negligible  effect  on  the  critical  stresses  of  the  FJ-4 
wing. 

10.  Correlation  of  calculated  postshot  gamma  doses  with  the  measured  values  were  good 
when  the  effects  of  shielding  were  considered. 

11.  Correlation  of  both  the  free  stream  overpressures  and  times  of  shock  arrival  with 
their  measured  values  was  excellent,  thus  confirming  the  accuracy  of  the  prediction  methods. 

12.  The  data  analyzed  showed  engine  performance  variations,  but  no  detrimental  effects 
were  experienced.  The  engine  environmental  pressure  did  not  vary  in  an  oscillatory  pattern, 
and  compressor  case  differentials  were  well  within  the  normal  values  experienced  over  gen¬ 
eral  operation  of  the  engine.  B  is  therefore  concluded  that  no  power  plant  structural  damage  or 
adverse  operational  effects  will  occur  in  the  vicinity  of  low  yield  nuclear  detonations  for  air¬ 
craft  positions  limited  by  other  criteria. 
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5.2  RECOMMENDATIONS 


It  is  recommended  that: 

1.  Additional  study  be  given  the  problem  of  predicting  radiant  exposure  and  irradiance. 

S.  Time-independent  solutions  for  the  prediction  of  the  maximum  temperature  rise  in 
honeycomb  panels  be  developed. 

3.  A  general  review  of  skin  absorptivitles  be  made  taking  into  account  the  wave  length  en¬ 
ergy  distribution  as  a  function  of  time  for  nuclear  pulses. 

4.  Further  flight  testing  of  the  FJ-4  should  be  performed  in  the  vicinity  of  nuclear  detona¬ 
tions  corresponding  to  low  blast  angle  and  high  overpressure  conditions. 

5.  Evaluation  of  modal  stresses  for  complex  structures  should  be  made  experimentally  by 
incorporating  strain-gage  instrumentation  in  conjunction  with  conventional  ground-vibration 
test  procedures. 

6.  A  comprehensive  investigation  should  be  performed  to  establish  a  more  precise  knowl¬ 
edge  of  the  time  history  of  material  velocity  behind  the  shock  front. 

7.  A  comprehensive  investigation  should  be  performed  to  establish  a  more  thorough 
knowledge  of  the  diffraction  phenomenon. 

8.  Additional  engine  response  data  be  collected  during  future  high  yield  tests,  since  the 
conclusions  made  are  based  upon  data  from  low  yield  tests  only. 
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Appendix  A 

THERMAL  EFFECTS  THEORY 


A.1  COMPUTATIONAL  PROCEDURE  FOR  RADIANT  EXPOSURE 


The  following  equation  was  developed  for  the  calculation  of  radiant  exposure  normal  to  the 
critical  surface: 


3.77  x  10T  W*M  A'  (1-F) 
=  R*eoR. 


(A.1) 


where  QTn 
W 
R 

R. 

a 

(1-F) 

7 


K» 

Vr 

i 

*ac 

A' 


ieff 


B 


Hb 

D 


radiant  exposure,  cal/cm* 

radiochemical  yield  of  the  device  considered,  Ulotons 
slant  range  from  burst  to  test  aircraft,  feet 
R,  1000,  feet 

atmospheric  attenuation  factor,  per  10"*  feet.  Figure  A.1 

flyaway  factor 

f(rj  x  Kj),  Reference  8 

time  to  peak  irradiance 

0.032  W\  seconds.  Reference  1 


radial  component  of  velocity,  10~*  ft/sec 
x  cos  i  x  10~*  Reference  8 

angle  of  incidence  of  irradiance  to  the  plane  of  the  critical  surface,  degrees 
true  airspeed,  ft,  sec 
non-dimensional  normalising  factor 

K.sinieff  +  ~B  ■(— ) 

i  where  i  a  20  degrees 
20  degrees  where  i  <  20  degrees 
ground  reflectance  (albedo)  =  0.44 

reflected  normal  coefficient  for  an  air  burst 

f  (Bk,,  D,  Hb),  Reference  11 

height  of  test  aircraft  above  ground,  feet 

height  of  burst,  feet 

horizontal  distance  of  test  aircraft  from  burst,  feet 
incident  correction  factor 
2 

fj  +  f,  •  fj  •  cosy  (90-i),  Reference  12 
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Figure  A.1  Atmospheric  attenuation  for  an  exceptionally  clear 
day,  50  miles  visibility,  versus  aircraft  altitude.  (From  unpub¬ 
lished  data  supplied  by  W.  E.  Schoor,  Sandia  Corporation). 


f.  =  6.21  x  10*  -  2.0  x  10-*],  for  177  s  369 

lHb’  J  W* 

Bk 

f,  =  1,  for  <  177,  surface  burst 
Hh 

f,  =  0,  for  ^5^  >  369,  air  burst 
*t  =  l-f> 

f,  =  3.41  x  10-’  (|b.)  -  1,  for  1*.  <  183 


f,  =  1.63,  for  —g  2  183 
i  ws 


J  = 


0.5  D* 


(1  +  sin  i)  (Dg  -  sin  i)1 

Dl* 


Kj  =  reflection  interpolation  factor 

Since  J  is  only  valid  for  a  hemispherical  fireball,  the  factor  Kj  is  Introduced  to  provide  for 
linear  interpolation  between  Hb  =  0  and  Hj,  =  R^,  where  Rs  is  the  radius  of  the  hemispherical 
fireball,  given  by 


R,  =  230  W4*,  feet,  Reference  11. 
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Kj  must  be  defined  so  that 

K,  =  1  when  Hb  =  0 
and 

Kj  =  J  when  Hb  =  Rs 
The  equation 


satisfies  these  conditions,  and  therefore 


Kj  = 


Ra  J 


*  Hb  (1  -  J)  +  Rs  J 


,  or 


£  _  Hb  (1  -  J)  +  Rs  J 


for  Hb  <  Rs 


-  =  lfor  Hb  2:  Rs 


Equation  A.1  was  used  to  calculate  radiant  exposure  for  comparison  with  measured  radi¬ 
ant  exposure. 


A.2  COMPUTATIONAL  PROCEDURE  FOR  IRRAD LANCE 


Equation  A.2  was  used  to  calculate  the  lrradiance-time  histories. 


,(t)  =  4x[R(t)]1  (929)  e«*iW 


(A.2) 


where  q(t)  =  irradiance  at  time  t,  cal/cml-sec 
qm  =  8.20  x  W*'44  x  10u,  cal/sec 
W  =  radiochemical  yield  of  the  device  considered,  kt 
A'  =  nondimensional  normalising  factor  (see  Equation  A.1) 
a  =  atmospheric  attenuation  coefficient,  per  10“*  feet 
R(t)  =  slant  range  from  the  center  of  burst  to  the  test  aircraft  at  time  t,  feet 

=  1  (H«-Hb)*  +  (D+VM  -'--v)1 

H|g  =  altitude  of  test  aircraft,  feet 
Hb  =  altitude  of  burst,  feet 

D  =  horizontal  distance  of  test  aircraft  from  burst  at  time  zero,  feet 
V  =  true  airspeed,  ft/sec 


t  =  time,  seconds 
tj  =  time  to  peak  irradiance 
=  0.032  x  w'* ,  seconds 


It  was  assumed  that  the  variation  in  all  factors  other  than  slant  range  and  q/qn  were  negligible 
during  the  duration  of  the  thermal  pulse.  Therefore,  these  factors  were  evaluated  for  the  test 
aircraft’s  position  at  time  zero  and  were  assumed  constant.  Computational  procedure  involved 
arbitrarily  selecting  values  of  t,  and  calculating  q(t)  versus  t. 


87 

CONFIDENTIAL 


An  IBM  704  Fortran  program  was  compiled  to  calculate  q(t)  versus  t  using  Equation  A.2. 
The  number  of  calculated  values  of  q(t)  versus  t  is  dependent  upon  the  number  of  points  from 
the  (q/qm)  versus  (t/rj)  curve  which  is  loaded  into  the  program  as  a  table.  The  program  also 
performs  a  stepwise  trapezoidal  integration  to  provide  a  calculated  radiant  exposure  time 
history. 


A.  3  THERMAL  RESPONSE 


A.3.1  Thin  Skin  Maximum  Temperature  Rise.  Since  the  conventional  skin-stringer  struc¬ 
ture  at  the  selected  thermocouple  locations  in  Table  3.6  satisfied  the  thin  skin  criteria  of 
Reference  9.  the  maximum  temperature  rise  was  calculated  by  Equation  A.3. 


AT 


m 


Qtn  '  >'  *  (1-FH) 
G  (1-F) 


,  degrees  F 


(A.3) 


where  ATm 

Qt 

> 

G 

Ps 

Cp 

bs 

(1-F) 

(1-FH) 


h. 


maximum  temperature  rise,  degrees  F 

radiant  exposure  from  Equation  A.1,  cal/cm* 

skin  absorption  coefficient 

PgCpbs,  cal,  cm1  degrees  F 

skin  density,  gm/cm* 

specific  heat,  cal/gm  degrees  F 

skin  thickness,  cm 

flyaway  factor  for  horizontal  flight  (See  Equation  A.1) 

combined  convective  cooling  and  flyaway  factor  for  horizontal  flight 

F^qKj,  Reference  8 

turbulent -heat -transfer-coefficient,  cal/cm2-sec -degrees  F,  (see  Equation  A.  8) 
and  K]  is  defined  under  Equation  A.1 


The  ratio  (1-FH)/  (1-F)  is  used  since  it  was  determined  that  the  flyaway  and  convective  codling 
effects  are  inter -dependent,  and  (1-F)  is  divided  into  QTN  to  cancel  the  (1-F)  incorporated  in 
Equation  A.1. 

A.  3.2  Time  to  Maximum  Temperature  Rise  in  Thin  Skin.  For  structure  satisfying  the 
thin  skin  criteria,  the  time  to  maximum  temperature  can  be  calculated  from  a  plot  of  (t^/q) 
versus  (hgjj/G)  presented  in  Reference  9.  hgif/G  is  the  same  term  as  is  calculated  in  obtaining 
(1-FH)  in  Equation  A.3.  tm  (time  to  ATm»)  is  necessary  in  the  prediction  of  temperature  at 
shock  arrival  which  is  presented  in  A.3.3. 


A.3.3  Temperature  at  Shock  Arrival  in  Thin  Skin.  The  calculation  of  the  temperature  at 
shock  arrival  is  based  on  three  primary  assumptions: 

(1)  The  structure  satisfies  the  thin  skin  criteria, 

(2)  Heat  flow  along  the  skin  is  negligible,  and 

(3)  Heat  loss  by  radiation  is  negligible. 


B  is  then  possible  to  write  a  simple  heat -balance  equation  around  such  a  skin  of  the  form: 
Heat  Out  =  Heat  In  +  Loss  of  Heat  Stored,  or: 


h8  (T(t)  -  TBL)  =  q i  -  (PgCpbg) 
The  internal  source  of  heat  must  be 


(A.4) 


qt=  hg(T,  -  Tbl) 


(A.5) 


where  T,  =  temperature  of  the  skin  at  time  zero,  degrees  F 
T  BL  =  boundary  layer  temperature,  degrees  F 
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Substituting  Equation  (A.4)  into  Equation  (A. 5),  and  rearranging 


dT  ^  h.T  _  h,,(T,  -  Tbl)  t  h„TBL 

®  PsCpb8  pBCpbg  PsCpbg 


or 


dT  ,  h„  T  _  hR  T, 

dt  pgCpbg  PgCpbg 


(A.6) 


where  t  =  time,  seconds 

T  =  temperature  of  skin  at  time  t,  degrees  F 

Substituting  the  .boundary  conditions  of  T.„  and  t^,  the  solution  to  the  differential  equa¬ 
tion  (Equation  A.6)  for  temperature  at  shock  arrival  is: 


AT 


=  exp 


f  bs 

IpsCpbg 


(tea 


(A.7) 


where  TM  =  temperature  at  time  of  shock  arrival,  degrees  F 
tgg  =  time  of  shock  arrival,  seconds 


ATmaz  and  t„,,»  are  found  in  Sections  A. 3.1  and  A3.2  respectively. 


A.3.4  Turbulent-Heat-Transfer-Coefficient.  The  turbulent -heat -transfer-coefficient  at 
each  selected  thermocouple  location  can  be  calculated  with  a  reasonable  degree  of  accuracy  by 
assuming  that  the  thermocouples  were  located  on  a  flat  plate  having  the  same  length  as  the 
actual  skin.  The  equation  presented  below  is  used  to  calculate  the  theoretical  h8. 


he  =  .00803  ~  Re'*  Pr'*1  cal/cn^-sec -degree  F 

Li 


where 


Re  =  P  Yes  —  =  Reynolds  Number,  dimensionless 

p  =  density,  slugs/ft* 

Vac  =  true  airspeed,  ft/sec 
L  =  characteristic  length,  feet 
p  -  viscosity,  slugs/ft -sec 

Pr  =  ^  =  Prandtl  Number,  dimensionless 

Cp  =  specific  heat,  BTU/  Slug -degree  F 
K  =  thermal  conductivity,  BTU -ft /ft 2 -sec -degree  F 


(A.8) 


Since  all  of  the  above  properties  are  based  on  the  temperature  of  the  boundary  layer,  the  fol¬ 
lowing  equation  is  used  to  calculate  this  reference  temperature,  T*. 


T*  =  0.5  (Ta  +  TW)  +0.162 


where  Ta  =  free  stream  static  temperature,  degrees  Ra 
Tw  =  turbulent  boundary  layer  temperature 

=  T,  +  0.746  ,  degrees  R. 

The  turbulent -heat-transfer-coefficient  can  also  be  calculated  by  an  equation  obtained 
from  rearranging  Equation  A.7. 
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(A.9) 


Computational  procedure  involved  selection  of  an  arbitrary  value  for  T  (or  t)  from  the  decay 
portion  of  the  temperature-time  history,  finding  the  corresponding  value  for  t  (or  T),  and 
solving  Equation  A.9  for  hs. 

A.3.5  Temperature-Time  ffistory  of  a  Thin  Skin.  Based  on  the  assumptions  that  the  skin 
is  a  thin  plate  with  an  adiabatic  wall  at  the  inner  surface  and  heat  flow  along  the  skin  is  negli¬ 
gible,  a  simple  heat -balance  equation  around  the  skin  can  be  written  in  the  form: 

Heat  In  =  Heat  loss  by  convection  +  heat  loss  by  radiation  +  gain  in  storage 

or 

«l(t)  +  qj  =  hs(T  -  Tbl)  +  otT4  +  PsCpb8^  (A.10) 

from  Equation  A.  5 

«li  =  hs(T,  -Tbl) 

Equation  A.10  then  becomes 

q(t)  =  ha(T  -  T,)  +  oe T4  +  p8Cpb8  (A.11) 

Since  no  mathematical  relationship  exists  for  q(t).  Equation  A.11  is  solved  for  T  versus  t  by 
the  method  of  finite  differences  in  the  following  manner: 

assume 

dT  _  Tn  —  T,,-!  _  T„  -  Tn_i 
dt  ~  tn  —  tn_j  "  At 

2 

T4=T4n-l 

where  Tn  =  temperature  at  tn 

Tn_i  =  temperature  at  t„  - 1  =  tn_j 

The  above  assumptions  linearize  Equation  A.11.  Substituting  these  relations  into  Equation  A.11, 
results  in  equation  A.  12. 


where  T„  *  temperature  of  skin  at  time  degrees  Rankine 

qn  =  ir radiance  at  time  tn,  cal/cm2-sec 
qn-i  =  irradiance  at  time  t„_i ,  cal/cm*-sec 
>  =  skin  absorptivity 

T,  =  initial  temperature  of  the  skin,  degrees  Rankine 
hs  =  turbulent  heat  transfer  coefficient,  cal/cm*-sec 
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7.52  x  10-1  =  factor  for  conversion  of  the  radiation  loss  term  from  btu/ft*-hr  to  cal/ 
cm*-sec 

a  =  Stephan  Boltzmann  coefficient 

G  =  PaCptg 

e  =  emissivity  for  re-radiation  in  the  infrared  wave  length  region 

An  IBM  704  Fortran  program  was  written  to  solve  the  finite  difference.  Equation  A.12. 

Computational  procedures  involve  the  loading  of  the  physical  and  thermal  properties  of 
the  skin  considered  and  selected  values  of  q(t)  from  the  thermal  pulse  as  input  data.  The 
program  then  computes  the  temperature  time  history  of  the  skin. 

A.3.6  Honeycomb  Thermal  Response.  Because  of  the  complexity  of  honeycomb  panel  con¬ 
struction,  temperature -time  history  solutions  are  required  to  calculate  the  thermal  response 
to  a  desirable  degree  of  accuracy.  The  temperature-time  history  is  solved  through  a  one¬ 
dimensional  heat  flow  IBM  704  program.  A  heat  balance  is  written  around  each  of  the  follow¬ 
ing  elements:  (1)  face  plate,  (2)  face  plate  adhesive,  (3)  core,  (4)  back  plate  adhesive,  and 
(5)  back  plate.  This  results  in  five  simultaneous  equations  which  are  solved  by  the  IBM  704  by 
matrix  mathematics.  Loading  the  physical  and  thermal  properties  of  the  Individual  components 
and  the  thermal  pulse  as  input  data,  the  program  computes  the  temperature-time  history  of 
each  element 

Many  of  the  variables  involved  in  the  honeycomb  thermal  response  could  be  eliminated  for 
the  specialised  conditions  in  this  operation  as:  (1)  only  one  type  of  honeycomb,  having  a  */u-inch 
thick  core,  an  0.016-inch  thick  exposed  face  plate,  and  an  0.008-to-0.10-inch  thick  back  plate, 
is  used  on  the  FJ-4  airplane.  (2)  The  range  of  Initial  temperatures,  boundary  layer  tempera¬ 
tures,  and  turbulent  heat  transfer  coefficients  were  relatively  constant.  (3)  For  the  test  yields 
of  less  than  80  kt,  the  thermal  pulse  duration  was  short,  with  lOij  being  less  than  3  seconds. 

It  was  found  that  the  maximum  temperature  rise  could  be  expressed  as  a  function  of  total 
heat  absorbed  with  a  fair  degree  of  accuracy.  Variations  in  Items  2  and  3  above  had  less  than 
a  5  percent  effect  on  the  maximum  temperature  rise.  Thus,  the  following  equation  was  derived 
for  subsonic,  ICAO  standard-day,  low-altitude  conditions  with  yields  of  less  than  80  kt. 

ATffi  =  32.0  r  QtN  (A.18) 

where  ATm  =  maximum  temperature  rise,  degrees  F 
32.0  =  constant,  °F  —  cm*/cal 

QTN  =  radiant  exposure  from  Equation  A.1,  cal/cm2 
y  -  skin  absorption  coefficient 

It  should  be  emphasized  that  Equation  A.  13  is  not  for  general  application  to  the  problem  of 
temperature  rise  in  honeycomb  but  is  valid  for  the  specific  honeycomb  structure  configuration 
of  the  FJ-4B. 
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Appendix  B 


DYNAMIC  RESPONSE  THEORY 


The  purpose  of  the  following  theoretical  development  was  to  determine,  by  a  comprehen¬ 
sive  analysis,  the  detailed  motions  of  the  FJ-4  airplane  In  response  to  a  nuclear  blast  wave. 
This  included  not  only  the  motions  of  the  airplane  in  space  but  also  the  transient  vibratory  mo¬ 
tions  of  the  structure.  These  latter  motions  served  as  the  basic  means  of  evaluating  the  stress 
levels  existing  throughout  the  air  frame  following  blast  impingement. 


B.l  DEVICE  EFFECTS  PARAMETERS 


In  order  to  analyze  the  dynamic  response  of  aircraft  structures  to  the  blast  wave 
from  a  nuclear  detonation,  it  is  necessary  to  predict  the  variation  of  the  significant  blast  wave 
parameters  with  distance  and  time  behind  the  shock  front. 


B.1.1  Particle  Velocity  Behind  the  Shock  Front.  The  peak  velocity  of  the  air  behind  the 
shock  front  was  determined  from  Reference  1  using  the  equation: 


where  u  =  particle  velocity  behind  the  shock  front,  ft/sec 

Cac  =  sound  velocity  in  front  of  the  shock  at  receiver  altitude,  ft/sec 
Ap  =  peak  overpressure  behind  shock  front,  psi 
pM  =  ambient  pressure  at  receiver  altitude,  psi 

Values  of  Ap  were  determined  by  the  equations  described  in  Appendix  D.  Values  of  pac  and 
Cac  were  based  on  the  ICAO  standard  atmosphere. 


B.1.2  Peak  Density  Behind  the  Shock  Front.  The  peak  density  behind  the  shock  front  was 
determined  by  the  following  relationship: 


where  p'  =  peak  density  behind  shock  front,  slugs/ft1 

p^  =  ambient  density  at  receiver  altitude,  slugs/ ft* 
Ap  =  overpressure  behind  shock  front,  psi 
pM=  ambient  pressure  at  receiver  altitude,  psi 


(B.2) 
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B.1.3  Peak  Dynamic  Pressure  Behind  the  Shock  Front.  The  maximum  value  of  the  dy¬ 
namic  pressure,  q,  due  to  the  particle  velocity  behind  the  shock  front,  was  determined  by  the 
relationship: 

q  =  l/i  P'(u)1,  psf  (B.3) 


B.1.4  Time  History  of  the  Shock  Front.  Values  to  t0,  free  air  duration  of  positive  phase, 
were  determined  from  the  following  equations  (Reference  10): 


(tf 

(B.4) 

-s  a)®' (&r 

(B.5) 

t^  =f(Ri> 

(B.6) 

where 


R  =  slant  range  to  the  test  aircraft,  feet 

R,  =  slant  range  for  a  1  kt  burst  in  a  homogeneous  sea-level  atmosphere,  feet 
W  =  radiochemical  yield  of  the  device  under  consideration,  kt,  or  twice  the  radiochemi¬ 
cal  yield,  kt,  as  dictated  by  the  triple-point  path  criterion. 

Wt  =  1  kt 

p  #  =  ambient  pressure  at  sea  level,  psi 
pK  =  ambient  pressure  at  the  altitude  of  the  test  aircraft,  psi 
C,  =  velocity  of  sound  at  sea  level,  ft/sec 
C*c  =  velocity  of  sound  at  the  altitude  of  the  test  aircraft,  ft/sec 
t(j  =  positive  phase  duration  for  a  1  kt  burst  in  a  homogeneous  sea-level  atmosphere, 
seconds 

t,  =  positive  phase  duration  of  a  burst  of  W  yield  at  the  altitude  of  the  test  aircraft, 
seconds 


The  functional  relationship  t0j  =  f(R,)  is  presented  in  graphical  form  in  Reference  4. 

Raving  found  the  values  of  t0,  the  time  history  of  the  shock  front  was  found  by  use  of  the 
following  equation: 


Ap{t)  =  Ap  1 


(B.7) 


where  p(t)  =  free  stream  overpressure  at  time  t,  psi 
t  =  time,  sec,  0  <  t  t0 

The  values  for  Ap  were  found  by  the  equations  described  in  Appendix  D. 


B.1.5  Time  of  Arrival  of  the  Shock  Front.  Shock  time  of  arrival  can  be  determined  from 
the  following  equations  which  were  obtained  from  Reference  10. 


where 


(B.8) 

(a#) 

(B.10) 

time  of  arrival  of  shock  front,  seconds 

time  of  arrival  of  shock  front  of  a  1  kt  burst  in  a  homogeneous  sea  level  atmos¬ 
phere,  seconds 
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(R,  Rj,  W,  Wl(  p*  pK,  C„  and  C,c  have  been  previously  defined  in  paragraph  B.1.4.)  The  func 
ttonal  relationship  t8,  =  f(R,)  is  found  in  graphical  form  in  Reference  1. 

The  method  outlined  in  Appendix  D  requires  that  twice  the  actual  yield  be  used  if  a  ground 
burst  were  considered  or  if  the  test  aircraft  were  outside  the  triple-point -path.  Since  Project 
5.3  did  not  participate  in  any  ground  burst  shots,  only  the  triple -point -path  criterion  was  con¬ 
sidered. 

B.1.6  Shock  Front  Propagation  Velocity.  Having  previously  determined  the  particle 
velocity  of  air  behind  the  shock  front  (Section  B.1.2),  the  shock  front  velocity  was  determined 
from  the  following  relationship. 


shock  front  propagation  velocity,  ft/sec 
particle  velocity  behind  the  shock  front,  ft/sec 
peak  overpressure  behind  the  shock  front,  psi 
ambient  pressure  at  the  altitude  of  the  test  aircraft,  psi 


where  U  = 
u  = 
Ap  = 
P«C  = 


B.2  DYNAMIC-RESPONSE  EQUATIONS 

B.2.1  Discussion.  The  purpose  of  this  section  is  to  describe  the  basic  considerations  in¬ 
volved  in  the  establishment  of  the  equations  that  were  used  to  analyze  an  elastic  airplane  in 
response  to  a  blast  wave.  As  the  wave  front  rapidly  envelops  the  airplane,  the  sudden  pressure 
rise  and  associated  material  velocity  cause  dynamic  loads  to  be  developed  of  extremely  short 
rise  time  and  of  sufficient  intensity  to  cause  a  significant  excitation  of  structural  vibrations  in 
combination  with  deviations  in  flight  path  and  altitude.  The  problem  is  therefore  in  the  form  of 
a  special  type  of  gust  loading  condition  in  combination  with  secondary  effects  associated  with 
the  pressure  rise.  The  method  of  analysis  would  be  equally  applicable  to  cases  in  which  the 
pressure-rise  effects  were  more  dominant  than  the  gust-loading  effects. 

The  analytical  work  for  this  investigation  was  performed  in  two  stages.  Prior  to  the  test 
series,  the  first  analysis  was  performed  and  the  equations  that  were  derived  at  that  time  are 
referred  to  as  the  original  equations.  Subsequent  to  the  test  series,  a  second  analysis  was 
performed,  leading  to  the  development  of  what  are  referred  to  as  the  modified  equations.  The 
modified  equations  incorporate  a  number  of  improvements  to  eliminate  deficiencies  in  the 
original  analysis  that  were  revealed  by  the  test  results.  In  addition,  these  dynamic -response 
equations  were  reformulated  in  a  manner  to  take  advantage  of  improved  procedures  that  have 
evolved  in  the  course  of  development  concurrently  being  performed  by  the  Dynamic  Science 
Section  of  NAA  in  the  field  of  steady-state  aeroelasticlty. 

The  basic  formulations  of  both  the  original  and  modified  equations  were  similar  in  that 
two  rigid-body  degrees  of  freedom,  heaving  and  pitching,  and  three  flexural  degrees  of  free¬ 
dom,  the  first  three  symmetrical  wing  vibration  modes,  were  used  to  represent  the  dynamic 
system.  The  heaving  and  pitching  modes  were  represented  in  the  form  of  conventional 
longitudinal-stability  equations  in  order  that  these  rigid-body  motions  may  properly  influence 
the  structural  equations  which  were  derived  by  application  of  Lagrange’s  Equations  to  the  first 
three  orthogonal  or  normal  wing  vibration  modes.  This  description  of  the  degrees  of  freedom 
for  the  dynamic  system  constitutes  the  left-hand  side  of  the  equations,  in  which  the  modified 
equations  differed  from  the  original  set  only  by  a  new  evaluation  of  the  generalized  masses  and 
frequencies  which  reflected  the  latest  available  data  as  pertained  to  the  specific  wing  mass 
configuration  of  the  test  operation. 

The  representation  of  the  forces  acting  on  the  system;  however,  were  quite  different  in 
the  two  analyses.  These  forces  constitute  all  the  terms  on  the  right-hand  side  of  the  equations. 
The  basic  concept  of  the  original  analysis  was  to  divide  the  airplane  into  seven  lifting-surface 
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components,  six  streamwise  wing  sections  and  the  horizontal  tall,  and  to  calculate  the  total 
applied  force  acting  on  each  section  at  any  instant  due  to  the  combination  of  the  blast  wave  and 
the  resulting  airplane  and  structural  motions.  This  approach  necessitated  the  utilization  of 
aerodynamic  strip  theory  with  empirical  three-dimensional  spanwlse  corrections  for  the  six 
wing  sections.  The  response  data  recorded  during  the  test  operation  indicated  that  this  sec¬ 
tional  approach  to  the  representation  of  the  applied  air  forces  inadequately  described  the  actual 
rigid-body  and  structural  response  characteristics  of  the  FJ-4.  On  the  basis  of  these  findings, 
another  approach  was  adopted  for  the  modified  analysis.  This  new  method  separately  treated 
each  applied  force  arising  from  several  independent  causes.  Each  of  these  forces  was  then  al¬ 
lowed  to  act  on  the  entire  structure  in  each  of  its  five  degrees  of  freedom  rather  than  inde¬ 
pendently  affecting  the  seven  components  of  the  previous  analysis.  These  causes  can  be  cate¬ 
gorized  as  follows:  initial  blast  inputs,  induced  aircraft  rigid  motions,  and  induced  structural 
deformations  and  velocities.  The  manner  in  which  the  forces  were  computed  due  to  each  of 
these  causes  will  be  discussed  individually  in  the  following  section. 

B.2.2  Summary.  The  dynamic -response  equations  are  illustrated  below  in  a  very  general 
form.  Each  form  of  these  equations  will  be  technically  described  in  the  following  paragraphs. 
Heaving  Equation: 

MV (9  -a)  =  LQa(t)  +  lJ (t)  +  L*(t  -  tw) 

Pitching  Equation: 

13  =  Paar(t)  +  P^a(t)  +  P^9(t)  +  P66(t)  +  Pj(t)  +  Pj(t  -t^ 

First  Modal  Equation: 

G,(«i  +  g^,  +  «?*,)  =  S,(t)  +  S1/?(t  -tw)  +  Ey(t  -tw) 

Second  Modal  Equation: 

G,(I,  +  +  ulti)  =  s*(t)  +  EJjg(t  —  tw)  +  —  tw) 

Third  Modal  Equation: 

G,(I,  +  gjb-'slj  +  =  **(t)  +  Hjp(t  —  tw)  +  E}£>(t  -tw) 

The  left-hand  terms  of  the  heaving  and  pitching  equations,  respectively,  represented  the 
total  airplane  inertial  forces  acting  normal  to  the  flight  path  and  inertial  moments  acting  about 
the  center  of  gravity.  These  terms  depended  directly  upon  the  airplane  mass  characteristics 
and  forward  speed  as  functions  of  the  airplane  rigid-body  accelerations. 

The  first  term  on  the  right-hand  side  of  the  heaving  equation  described  the  total  lift  force 
acting  on  the  airplane  due  to  perturbations  of  the  airplane  angle  of  attack.  This  function  was 
based  on  the  experimentally-determined  airplane  lift  coefficient  which  incorporated  the  effect 
of  load  alleviation  due  to  quasi -static  structural  deformations  of  the  air  frame. 

The  first  three  terms  on  the  right-hand  side  of  the  pitching  equation  represented  the  total 
pitching  moment  acting  on  the  airplane  about  its  center  of  gravity  due  to  perturbations  of  the 
airplane  angle  of  attack,  angle -of -attack  rate  and  angle-of-pitch  rate,  respectively.  All  these 
functions  were  based  on  the  corresponding  experimentally-determined  longitudinal-stability 
derivatives  which  also  incorporated  the  effect  of  load  alleviation  due  to  quasi-static  structural 
deformations  of  the  airframe. 

The  fourth  term  on  the  right-hand  side  of  the  pitching  equation  described  the  additional 
airplane  pitching  moment  due  to  rigid  motions  of  the  horizontal  stabilizer.  This  term  also  in¬ 
corporated  the  quasi-static  structural  alleviation  of  the  air  loads. 

The  second  and  third  terms  on  the  right-hand  side  of  the  heaving  equation,  plus  the  fifth 
and  sixth  terms  on  the  right-hand  side  of  the  pitching  equation  respectively,  represented  the 

95 


CONFIDENTIAL 


total  blast-induced  forces  and  moments  acting  on  the  airplane.  In  each  of  these  equations,  the 
blast -induced  load  was  described  by  two  separate  terms  as  Indicated.  The  first  terms  denoted 
the  blast -induced  loading  on  the  horizontal  stabilizer  while  the  second  terms  depicted  a  cor¬ 
responding  load  acting  on  the  wing-fuselage  combination.  These  latter  terms  incorporated  a 
finite  time  delay  to  account  for  the  time  required  by  the  shock  front  to  travel  from  the  hori¬ 
zontal  stabilizer  to  the  effective  center  of  pressure  of  the  wing-fuselage  combination.  This 
time  delay  was  dependent  upon  the  propagation  speed  and  orientation  of  the  shock  front,  and 
the  true  speed  of  the  airplane.  These  blast-induced  forces  and  moments  were  calculated  by 
converting  blast-induced  angles  of  attack  into  circulatory  lift  forces  using  aerodynamic  coef¬ 
ficient  data  applicable  to  the  separate  lifting  surfaces.  These  coefficients  incorporated  the 
quasi-static  structural  alleviation  of  air  loads  previously  mentioned.  The  blast-wave  gust  ef¬ 
fect  consisted  of  a  simultaneous  variation  of  local  angle  of  attack  and  dynamic  pressure  due  to 
the  material  velocity  and  over  density  behind  the  shock  front  as  prescribed  in  Reference  16. 

The  blast  wave  was  considered  to  instantaneously  envelop  the  entire  stabiliser  surface  and,  at 
a  later  instant  of  time,  to  Instantaneously  envelop  the  wing-fuselage  combination.  All  these 
blast -induced  forces  and  moments  Incorporated  the  time  lag  associated  with  the  build-up  of 
circulatory  lift  in  response  to  an  abrupt  variation  in  local  angle  of  attack  and  airspeed.  This 
lift  lag  was  incorporated  by  use  of  Duhamel’s  integral  in  conjunction  with  a  modified  Wagner 
Function  of  Reference  IS. 

All  three  modal  equations  were  identical  in  form  and  hereafter  will  be  described  as  a 
single  typical  equation  associated  with  its  corresponding  symmetrical  normal-vibration  mode. 
The  fundamental  approach  which  was  utilized  in  the  formulation  of  the  structural  equations  of 
motion  is  outlined  in  Reference  13  as  a  representation  of  an  unrestrained  elastic  airplane  in 
terms  of  normal  coordinates  by  means  of  Lagrange’s  Equation. 

All  the  terms  on  the  left-hand  side  of  the  modal  equation  represented  the  Inertial,  damping, 
and  elastic  properties  of  the  wing  structure  associated  with  the  normal  mode  being  considered. 
The  first  term  described  the  inertial  effect  of  the  generalized  mass  and  was  based  on  the  total 
kinetic  energy  of  the  modal  system.  The  second  term  represented  an  approximation  to  a  dis¬ 
sipative  force  as  a  means  of  removing  energy  from  the  system  due  to  structural  damping. 
Treatment  of  this  term  is  thoroughly  discussed  in  References  13  and  14.  The  last  term  on  the 
left-hand  side  of  the  modal  equation  represented  the  effective  elastic  restoring  force  of  the 
modal  system  and  was  based  on  the  total  strain  energy  of  the  structure  associated  with  the 
normal  mode  being  considered. 

The  first  term  on  the  right-hand  side  of  the  modal  equation  represented  the  summation  of 
all  the  generalized  aerodynamic -coupling  forces  associated  with  the  mode.  These  forces  arose 
from  the  induced  angle  of  attack  at  the  wing  due  to  both  the  rigid-body  displacements  and  rates, 
and  the  structural  deformations  and  rates  of  deformation,  described  as  Causes  2  and  3  in  Sec¬ 
tion  B.2.1.  It  was  through  this  term  that  the  rigid-body  motions  were  permitted  to  influence 
the  structural  responses  properly.  A  span-wise  variation  in  angle  of  attack  was  computed  for 
unit  amplitude  of  each  mode  of  structural  deformation  and  deformation  rate,  in  addition  to  the 
constant  angle  of  attack  associated  with  rigid-body  motions.  From  these  span-wise  distribu¬ 
tions  of  induced  angle  of  attack,  a  span-wise  variation  of  aerodynamic  load  coefficient  was 
computed.  A  conventional  19-point  method  was  incorporated  by  this  computation  which  was 
based  on  the  Weissinger  theory  for  subsonic  aerodynamic  load  distributions.  These  load  coef¬ 
ficient  functions  were  then  converted  to  their  generalized-force  contributions  to  each  struc¬ 
tural  mode.  The  summation  of  all  these  generalized  aerodynamic-coupling  forces  associated 
with  the  mode  was  entered  in  the  equations  in  a  Duhamel  integral  in  conjunction  with  the  modi¬ 
fied  Wagner  Function  previously  mentioned  to  account  for  the  time  lag  associated  with  the  cir¬ 
culatory  lift  build-up. 

The  second  term  on  the  right-hand  side  of  the  modal  equation  represented  the  generalized 
force  acting  on  the  wing  due  to  the  blast-wave  gust  effect.  This  generalized  force  was  funda¬ 
mentally  based  on  the  same  considerations  discussed  for  the  corresponding  terms  of  the 
heaving  and  pitching  equations.  The  aerodynamic -lift  coefficient  utilized  in  calculating  this 
load  was  for  the  exposed  wing  area  and  assumed  an  elliptical  span-wise  lift  distribution  without 
the  quasi -static  structural  alleviation  of  loading.  This  blast -induced  generalized  force  was 
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utilized  in  analytical  form  by  special  application  of  the  Duhamel  integral  to  incorporate  the  ef¬ 
fect  of  the  finite  period  of  time  for  the  blast  wave  to  sweep  across  the  wing  from  tip  to  root  as 
a  continuous  function,  rather  than  as  sequenced  inputs  to  six  wing  sections,  as  it  was  handled 
in  the  original  analysis.  The  comprehensive  treatment  of  this  continuous  sweeping  effect  of  the 
blast  wave  was  considered  to  be  the  most  original  innovation  of  all  the  changes  introduced  in 
the  modified  dynamic -response  analysis. 

The  last  term  on  the  right-hand  side  of  the  modal  equation  described  an  exponential  ap¬ 
proximation  to  the  diffraction  pulBe  measured  during  the  test  operation.  This  diffraction  load 
was  generalized  and  incorporated  in  the  sweeping  function  of  the  blast-wave  gust-effect  term 
as  an  additional  structural  loading.  This  effect  was  included  in  a  purely  empirical  manner, 
based  on  the  recorded  data  from  the  test  series,  enabling  an  evaluation  of  the  effect  of  the  dif¬ 
fraction  phenomenon  on  the  FJ-4  structural  response. 

In  summary,  it  can  be  seen  that  the  modified  dynamic  analysis  of  the  FJ-4  treated  the 
rigid  system  and  the  elastic  system  as  essentially  two  independent  dynamic  systems.  The  ef¬ 
fects  of  structural  deformation  on  rigid-body  forces  was  taken  into  account  by  the  use  of  flexi¬ 
ble  aerodynamic  derivatives.  This  takes  the  quasi-steady  mode  deformations  into  account, 
using  previously  established  data  confirmed  by  flight  test,  and  neglects  the  small  effects  on 
rigid  motions  caused  by  dynamic  excitation.  Thus,  it  is  unnecessary  to  include  structural 
coupling  terms  in  the  rigid-body  equations,  i.e.,  these  equations  contain  no  terms  including 
the  variables  |2,  |Jt  or  their  derivatives.  Outputs  from  the  rigid-body  equations  do  enter 
the  structural  equations,  however,  as  well  as  the  effects  of  each  mode  on  each  other.  Conse¬ 
quently,  rigid  aerodynamic  data  are  used  to  convert  angles  of  attack  to  forces  in  the  structural 
equations,  since  these  equations  automatically  calculate  their  own  flexibility  corrections.  This 
manner  of  separating  all  basic  input  forces  and  responses  has  been  found  to  provide  good  ac¬ 
curacy  and  efficiency  in  performing  solutions  on  an  analog  computer  and  analysis  of  the  sig¬ 
nificance  of  each  part  of  the  problem. 
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Appendix  C 

NUCLEAR  RADIATION  THEORY 


C.l  TOTAL  NUCLEAR  DOSE 


The  methods  tor  calculating  the  total  nuclear  radiation  dose  based  on  final  yields  and  ob¬ 
served  atmospheric  conditions  are  as  follows: 


Dt  =  Dy+  WDn 


where 


Dt  =  total  nuclear  dose,  rem 
Dy=  gamma  dose,  r 

=  f(Rac,  <rr,  W),  Reference  3 
Dn=  neutron  dose  for  1  kt,  rem 
=  f(Rac>  °r)>  Reference  2 
Rac  =  slant  range  from  aircraft  to  burst,  yards 
aT=  relative  air  density  between  aircraft  and  burst 
W  =  weapon  yield,  kt 


C.2  RELATIVE  AIR  DENSITY 


The  relative  air  density  used  in  this  analysis  is  defined  as  follows  (References  2  and  3): 
Given  air  pressures  at  burst  and  receiver  altitudes,  then 


25.8 P-  Pac 
hac  —  hb 


where  p^,  =  atmospheric  pressure  at  burst  altitude,  mb 

pac  =  atmospheric  pressure  at  receiver  altitude,  mb 
hac  ~  altitude  of  receiver,  feet 
hb  =  altitude  of  burst,  feet 

The  atmospheric  pressures  at  burst  altitudes  were  determined  by  applying  a  standard 
lapse  rate  of  33  mb  per  1,000  feet  altitude  to  the  known  values  at  ground  zero. 

When  an  aircraft  was  more  than  1,000  feet  short  of  and  approaching  ground  zero  at  time 
zero,  1  r  was  added  to  the  gamma  dose  to  account  for  the  aircraft’s  flying  toward  the  rising 
radioactive  cloud.  The  effects  of  flyaway  have  been  disregarded  as  a  conservatism. 
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Appendix  D 

OVERPRESSURE  THEORY 


The  method  used  for  computing  static  overpressure  was  based  on  weapon-effect  data  ob¬ 
tained  from  References  5  and  17. 

The  method  obtained  from  Reference  17  required  that  twice  the  actual  yield  be  used  for 
overpressure  prediction  if  a  ground  burst  were  considered  or  if  the  test  aircraft  were  outside 
the  triple-point  path.  Since  Project  5.3  did  not  participate  in  any  ground  burst  shots,  only  the 
triple -point -path  criterion  was  considered.  The  method  for  calculating  the  triple-point  path  for 
each  shot  is  defined  by  the  following  relationships: 


KD  _  /K^H,,  KYm\ 

w*  “  £VW’  w*  / 


(D.1) 


where  O  =  horizontal  distance  from  the  burst  of  a  point  on  the  path,  feet 
W  =  radiochemical  yield  of  the  device  under  consideration,  kt 


ph  =  ambient  pressure  at  the  altitude  of  a  point  on  the  path,  psi 
pb  =  ambient  pressure  at  the  burst  altitude,  psi 
p0  =  ambient  pressure  at  sea  level,  psi 
H,  =  height  of  burst  above  ground  zero,  feet 
Ym  =  H-Hg 

=  Mach  stem  height,  feet 
H  =  altitude  of  a  point  on  the  path,  feet 
Hg  =  altitude  of  the  ground,  feet 

The  functional  relationship  described  above  is  presented  in  graphical  form  in  Figure  D.l. 
By  selecting  a  series  of  H  values  greater  than  Hg,  corresponding  D  values  were  computed.  The 
coordinates  (Ym,  O)  then  defined  the  triple-point  path  desired. 

Static  overpressure  was  then  computed  for  positions  both  inside  and  outside  the  triple¬ 
point-path  by  the  use  of  the  following  equations  obtained  from  Reference  10. 


(D.3) 


AP!  =  f(Rj) 
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Scaled  Bunt  Haight,  K0Hb/W'$  ,  103  faat 


0  1  2  3^4  5  6  7  8 

KD/W  3  ,  io3  foot 


Figure  D.1  Triple-point  paths. 


Ap  = 


ApiPac 

P« 


(D.4) 


where  Ap,  =  peak  overpressure  for  a  1-kt  burst  In  a  homogeneous  sea -level  atmosphere,  psi 
Ap  =  peak  overpressure  at  altitude  of  the  test  aircraft 
p(  =  ambient  pressure  at  sea  level,  psi 
Pac  -  ambient  pressure  at  the  altitude  of  the  test  aircraft,  psi 
R,  =  slant  range  for  a  1-kt  burst  in  a  homogeneous  sea  level  atmosphere,  feet 
R  =  slant  range  to  the  test  aircraft,  feet 
W,  =  1  kt 


W  =  radiochemical  yield  of  the  device  under  consideration,  kt,  or  twice  the  radiochemi¬ 
cal  yield,  kt,  as  dictated  by  the  triple-point-path  criterion 


The  functional  relationship,  Ap,  =  f(R,)  is  presented  in  graphical  form  in  Figure  3.44. 
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Appendix  E 

INSTRUMENTATION 


The  instrumentation  specified  in  this  appendix  was  installed  by  the  contractor  identically 
in  two  FJ-4  airplanes,  (BUNO  139310,  NAA  No.  30,  and  BUNO  139467,  NAA  No.  87).  Primary 
positioning  equipment  consisted  of  a  Bendix-Pacific  X-Band  radar  beacon,  type  RBX-2,  and 
power  supply,  type  RPX-1,  used  to  insure  positive  lock-on  by  the  modified  M-33  gun  laying 
radars.  ARN-21  (TACAN)  was  also  installed  as  a  back  up  positioning  system. 

FJ-4  Aircraft,  BUNO  139467,  participated  in  Shots  Boltzmann,  Hood,  Diablo,  Kepler, 

Shasta,  and  Doppler.  FJ-4  Aircraft  BUNO  139310,  participated  in  Shots  Doppler  and  Smoky. 

The  flight  configurations  listed  in  this  appendix  are  applicable  to  all  shots  except  where  noted. 

E.l  PRIMARY  INSTRUMENTATION 

The  nose  compartment,  right  and  left  gun  bays,  and  all  ammo  bays  were  utilized  for  in¬ 
strumentation  equipment  installation.  A  Type  H,  200-gallon  wing  tank  was  utilized  for  additional 
oscillograph  and  related  equipment  Installation.  Two  CEC  26-channel  oscillographs  were 
mounted  in  the  nose,  two  similar  oscillographs  in  the  right  wing  tank,  and  a  12-hole  photo  re¬ 
corder  in  the  right  gun  bay. 

For  all  equipment,  except  the  structural  strain  gages,  calorimeters  and  radiometer,  the 
calibration  curves  were  published  in  Reference  18.  Calibration  curves  for  all  structural  strain 
gages  can  be  obtained  from  Reference  19.  Calibration  of  all  the  calorimeters  and  radiometers 
was  accomplished  by  NRDL.  Tables  E.5,  E.6,  E.7,  and  E.8  present  the  instrumentation  used  for 
the  collection  of  thermal,  overpressure,  structural  and  supporting  data,  respectively.  Nuclear 
radiation  data  was  acquired  from  film  badges  located  in  the  cockpit,  ammunition  bay,  right 
drop  tank,  and  nose  wheel  well.  In  Tables  E.5  through  E.8,  the  designations  of  “A,”  “B,”  "C,” 
and  “D”  for  the  oscillographs  refer  to  the  left  and  right  nose  compartment  installations  and  the 
rear  and  forward  wing  tank  installations,  respectively.  The  CF  numbers  shown  are  North 
American  Aviation  inventory  numbers  which  are  necessary  if  it  is  desirable  to  refer  to  the 
calibrations  listed  in  Reference  18. 

Figures  E.l  through  E.3  and  their  accompanying  legends,  and  Tables  E.l,  E.2,  E.3,  and 
E.4  show  the  installation  of  most  of  the  instrumentation. 

E.2  COMMENTS  AND  MISCELLANEOUS  INSTRUMENTATION 

The  following  is  a  list  of  items  of  miscellaneous  instrumentation  and  comments  concerning 
the  instrumentation  in  general. 

1.  Two  VDR-5  (16  mm)  cameras  were  installed  in  the  dorsal  fairing  to  photograph  left 
wing  and  empennage  flexure.  Painted  stations  and  flood  lighting  on  the  wing  and  marker  lights 
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TABLE  E.l  LEGEND  FOR  FIGURE  E.l  FOR  THERMOCOUPLE  INSTALLATIONS 


Thermo¬ 

couple 

Number 

Installation  Location 

NAA  Installation 
Drawing 
Number 

1 

Nose  Wheel  Well  Door,  Station  81.375 

209-75743 

2 

Aft  Engine  Access  Door,  Station  210.81 

209-75743 

3 

Skin,  Bottom  Center  Line,  Station  248.155 

209-75743 

4 

Inside  ARA  Antenna  Cover,  Station  313.625 

209-75743 

5 

Skin,  Bottom  Center  Line,  Station  333.625 

209-75743 

6 

Tail  Bumper  Bungee  Access  Door,  Station  389.125 

209-75743 

7 

Rudder  Splitter  Plate,  Inside  Left  Skin,  Water 

Plane  Station  69.375 

209-75742 

*8 

Rudder  Splitter  Plate,  Inside  Left  Skin,  Water 

Plane  Station  69.375 

209-75742 

0 

Right  Elevator,  Inside  Upper  Skin,  Station  45.875 

209-75742 

10 

Right  Elevator,  Inside  Lower  Skin,  Station  45.875 

209-75742 

11 

Right  Flap,  Inside  Lower  Skin,  Flap  Station 

59.233  Canted 

209-75744 

*12 

Right  Flap,  Inside  Lower  Skin,  Flap  Station 

59.233  Canted 

209-75744 

13 

Right  Wing  Rear  Spar,  Lower  Cap,  Rear  Spar 

Station  86.5 

209-75745 

14 

Right  Wing,  Inside  Lower  Skin,  Rear  Spar  Station 
138.454 

209-75745 

15 

Right  Wing,  Inside  Lower  Skin  Stiffener,  Rear 

Spar  Station  138.454 

209-75745 

16 

Right  Outer  Wing  Panel,  Outside  Lower  Skin,  Wing 
Station  222.125 

209-75745 

17 

Right  Outer  Wing  Panel,  Inside  Lower  Skin,  Wing 
Station  223.75 

209-75745 

18 

Right  Aileron,  Inside  Lower  Skin 

EOH-114948 

19 

Right  Outer  Wing  Panel  Inside  Lower  Skin  on  Blue 
Door,  Wing  Station  170 

20 

Nose,  Radome 

EOH-128454 

*  Backup  installation;  not  used. 
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TABLE  E.2  LEGEND  FOR  FIGURE  E.l  FOR  STRAIN  GAGE  INSTALLATIONS 


Strain  Gage 
Number 

Installation  Location 

NAA  Installation 
Drawing  Number 

1 

Fuselage  Axial  Stress,  Upper  Left  Longeron,  Out¬ 
board  Flange,  Station  250 

209-75736 

2 

Fuselage  Axial  Stress,  Upper  Left  Longeron,  In¬ 
board  Flange,  Station  250 

209-75736 

3 

Fuselage  Axial  Stress,  Upper  Left  Longeron,  Low  Web, 
Station  250 

209-75736 

4 

Fuselage  Axial  Stress,  Upper  Right  Longeron,  Out¬ 
board  Flange,  Station  250 

209-75736 

5 

Fuselage  Axial  Stress,  Upper  Right  Longeron,  In¬ 
board  Flange,  Station  250 

209-75736 

6 

Fuselage  Axial  Stress,  Upper  Right  Longeron,  Low 

Web,  Station  250 

209-75736 

7 

Fuselage  Torsion,  Left  Side,  Station  334.5 

209-75738 

8 

Fuselage  Torsion,  Right  Side,  Station  334.5 

209-75736 

9 

Fin  Bending,  Lower  Main  Beam,  Station  14.5 

209-75737 

10 

Fin  Torsion,  Forward  Spar,  Station  46.025 

209-75737 

11 

Fin  Torsion,  Forward  Spar,  Station  75.000 

209-75737 

12 

Fin  Torsion,  Rear  Spar,  Station  47.125 

209-75737 

13 

Fin  Torsion,  Rear  Spar,  Station  74.188 

209-75737 

14 

Fin  Bending,  Upper  Main  Beam,  Station  89.188 

209-75737 

15 

Fin  Cap  Bending,  Station  94.5 

209-75737 

*16 

Fin  Cap  Bending,  Station  99.228 

209-75739 

*17 

Fin  Cap  Bending,  Station  102.608 

209-75739 

*18 

Fin  Cap  Bending,  Station  105.328 

209-75739 

IS 

Rudder  Hinge  Moment 

-- 

20 

Horizontal  Stabilizer  Torsion,  Left,  Station  30.5 

209-75738 

21 

Horizontal  Stabilizer  Bending,  Left,  Station  33.75 

209-75738 

22 

Horizontal  Stabilizer  Torsion,  Left,  Station  58.5 

209-75738 

23 

Horizontal  Stabilizer  Bending,  Left,  Station  60.5 

209-75738 

24 

Horizontal  Stabilizer  Torsion,  Right,  Station  30.5 

209-75738 

25 

Horizontal  Stabilizer  Bending,  Right,  Station  33.75 

209-75738 

26 

Stabilizer  Actuator  Hinge  Moment 

— 

27 

Elevator  Hinge  Moment 

— 

28 

Horizontal  Stabilizer  Pivot  Load,  Left 

— 

29 

Wing  Bending,  Left  Outer  Wing  Panel,  Rear  Spar 

Station  173 

209-75804 

30 

Left  Aileron  Hinge  Moment 

— 

31 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  143.25 

209-75803 

32 

Wing  Shear,  Left  Rear  Spar,  Rear  Spar  Station  115.125 

209-75803 

33 

Wing  Bending,  Left  Rear  Spar,  Rear  Spar  Station  87.26 

209-75803 

34 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  86.2 

209-75803 

35 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  50.987 

209-75803 

36 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  63.612 

209-75803 

37 

Wing  Shear,  Left  Rear  Spar,  Rear  Spar  Station  55.051 

209-75803 

38 

Wing  Shear,  Left  Rear  Spar,  Rear  Spar  Station  67.801 

209-75803 

39 

Wing  Shear,  Left  Front  Spar,  Wing  Station  48 

209-75803 

40 

Wing  Shear,  Left  Landing  Gear  Front  Spar,  Wing 

Station  47.125 

209-75803 

41 

Wing  Shear,  Left  Landing  Gear  Rear  Spar,  Wing 

Station  39.84 

209-75803 

42 

Wing  Bending,  Left  Root,  Wing  Station  17.5 

209-75734 

43 

Wing  Bending,  Right  Root,  Wing  Station  17.5 

209-75734 

*44 

Left  Flap  Hinge  Moment 

— 

45 

Right  Flap  Hinge  Moment 

-- 

46 

Wing  Bending,  Right  Root,  Wing  Station  36.5 

— 

47 

Wing  Shear,  Right  Front  Spar,  Rear  Spar  Station  50.987 

209-75803 

48 

Wing  Shear,  Right  Front  Spar,  Rear  Spar  Station  63.612 

209-75803 

49 

Wing  Shear,  Right  Front  Spar,  Rear  Spar  Station  86.2 

209-75803 

*50 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  87.26 

209-75803 

51 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  55.051 

209-75803 

52 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  67.801 

209-75803 

53 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  115.125 

209-75803 

54 

Right  Aileron  Hinge  Moment 

— 

55 

Wing  Bending,  Right  Outer  Wing  Panel,  Rear  Spar 

Station  173 

209-75804 

56 

Horizontal  Stabilizer  Pivot  Load,  Right 

— 

*  Installation  not  used 
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TABLE  E.3  LEGEND  FOR  FIGURE  E.2 


Item 

Number 

Description  of  Item 

NAA  Installation 
Drawing  Number 

1 

Nose  Boom  and  Mount  Assembly 

194-75204-101 

2 

Dipole  Antenna  Radome 

209-931115 

3 

Calorimeter  and  Radiometer 

209-75721 

4 

Nose  Compartment  Instrumentation 

209-75048 

5 

Magnetic  Power  Supply 

209-75105 

6 

Vertical  Gyro 

209-75102 

7 

Cooling  Inverter 

209-75732 

8 

Right  Gun  Bay  Inverter  and  Electronic 

Disconnect  Panel 

209-75718 

9 

Doelcam  Gyro 

209-75118 

10 

Right  Gun  Bay  Instrumentation 

209-75103 

11 

Power  Distribution  Unit 

209-75110 

12 

Forward  Fuel  Cell  Electrical  Installations 

209-75507 

13 

External  Tank  Assembly 

Modified  Type  II  Instrumented  Tank  Assembly 
External  Tank  Cover  Assembly 

External  Tank  Cover  Latch  Assembly 

209-75726 

209-75727 

209-75728 

209-75729 

14 

Wing  Flexure  Lights 

209-75723 

15 

Wing  and  Empennage  Flexure  Camera 

Wing  and  Empennage  Flexure  Camera  Fairing 

209-75724 

209-75719 

16 

X-Band  Beacon  Equipment 

209-75127 

17 

Vertical  Stabilizer  Overpressure  Probe 

209-75741 

18 

Fin  Cap  X-Band  Antenna 

209-75126 

19 

Vertical  and  Horizontal  Stabilizer  Target 

Lights 

209-75722 

20 

Rudder  Position  Transmitter 

209-75611 

21 

Elevator  and  Rudder  Cone  Radar  Fairing 

L-Band  Recessed  Stub  APX-6  Antenna 

209-25007-21 

209-71622 

22 

Inboard  Elevator  Twist  Transmitter 

208-75616 

23 

Horizontal  Stabilizer  Position  Transmitter 

208-75604 

24 

Speed  Brake  Wells  Ballast 

209-75903 

25 

Accelerometer  at  Center  of  Gravity 

209-75117 

26 

Left  and  Right  Aileron  Position  Transmitter 

208-75627 

27 

Wing  Pressure  Survey  Transducers 

209-75758 

28 

External  Fuel  Tank  Float  Switch 

209-75716 

29 

Wing  Fuel  Gaging  Relay  Panel 

209-75714 

30 

Forward  Fuselage  Fuel  Cell  Capacitor 

209-75713 

*31 

TACAN  Computer -2  Units 

— 

32 

Ammunition  Can  Ballast 

208-75906 

33 

Left  Gun  Bay  Instrumentation 

209-75067 

34 

Ammunition  Bay  Ballast 

209-75902 

35 

Rudder  Force  Dynamometer  Ring  Assembly 

209-75701 

36 

Gun  Stabilization  Mounts  Ballast 

209-75905 

37 

Nose  Compartment  Plate  Assembly 

209-75111 

38 

Fuselage  Overpressure  Probe 

— 

*  Not  installed  in  final  configuration. 
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TABLE  E.4  LEGEND  FOR  FIGURE  E.3 


Item 

Number 

Description  of  Item 

NAA  Installation 
Drawing  Number 

1 

Tailpipe  Total  Temperature  and  Pressure 
Instrumentation 

209-75334-201 

2 

Compressor  Discharge  Pressure  Probe 

XFT-5433 

3 

Revere  Fuel  Flow  Transmitter 

209-75510 

4 

Power  Plant  Instrumentation 

209-75733 

5 

Compressor  Inlet  Pressure  and  Temperature 
Probes 

XFT-5432 

6 

MRP-12-4N  Magnetic  Power  Supply 

— 

7 

Miniature  18  Channel  Balance  Box 

XFT-5200 

8 

Miniature  18  Channel  Balance  Box 

XFT-5200 

9 

Calorimeter  and  Radiometer  Plugs 

— 

10 

RJ-22  Arnoux  Junction  Box 

— 

11 

Haydon  Timer 

— 

12 

Oscillograph  Disconnect  Panels 

XFT-5304 

13 

MRE-75S-G  and  MRE-75P-G  Quick  Disconnect 

Pull -Away  Plugs 

"** 

14 

USNRDL  MK-6F  Calorimeter 

209-75731 

USNRDL  MK-6F  Radiometer 

209-75731 

VDR-5  Camera 

209-75731 

15 

APX-6B  Transponder 

— 

16 

Gamma  Ray  Transducer 

— 

17 

APX-6  Antenna 

— 

.18 

5-114-P3  Modified  Oscillograph 

— 

19 

Miniature  18  Channel  Balance  Box 

XFT-5200 

20 

5-114-P3  Modified  Oscillograph 

— 

i 


i 
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TABLE  E.5  CHANNELIZATION  OF  THERMAL  INSTRUMENTATION 


Instrument  or  Parameter 

Units 

FJ-4  BuNo.  139467 

Recording 

Instrument  Channel 

FJ-4  BuNo.  139310 

Recording 

Instrument  Channel 

Calorimeter,  90°,  Forward  and  55°  Up 

cal/cm* 

Oscillograph  "A” 

25 

Oscillograph  “A” 

25 

Calorimeter,  180°,  Tank,  Painted  White 

cal/cmJ 

Oscillograph  “B" 

m 

Oscillograph  “B” 

19 

Calorimeter,  90°,  Tank,  Hot  Reference 

cal/cm* 

Oscillograph  "C” 

l 

Oscillograph  “C" 

1 

Calorimeter,  160°,  Tank,  Internal  Reference 

cal/cm* 

Oscillograph  “D" 

1 

Oscillograph  “D” 

1 

Radiometer,  180°,  Forward  and  Facing  Right 

cal/cm*/sec 

Oscillograph  “A” 

2 

Oscillograph  "A” 

2 

Radiometer,  90°,  Forward  and  55°  Up 

cal/cm*/sec 

Oscillograph  "A” 

26 

Oscillograph  “A" 

26 

Radiometer,  180°,  Tank 

cal/cm*/sec 

Oscillograph  “D” 

2 

Oscillograph  “DT 

2 

Outside  Air  Temperature,  CF72-45/72-47 

degrees  C 

Photo  Recorder 

- 

Photo  Recorder 

- 

Temperature,  Inside  ARA  Antenna  Cover,  Station  313.625 

degrees  F 

Oscillograph  "A" 

6* 

Oscillograph  “A” 

6 

Temperature,  Skin,  Aft  Engine  Access  Door,  Station  210.81 

degrees  F 

Oscillograph  "B" 

2 

Oscillograph  “B" 

2 

Temperature,  Skin,  Bottom  Center  Line,  Station  248.155 

degrees  F 

Oscillograph  “A* 

4 

Oscillograph  “A” 

4 

Temperature,  Skin,  Bottom  Center  Line,  Station  333.625 

degrees  F 

Oscillograph  **B” 

3 

Oscillograph  “B" 

3 

Temperature,  Tail  Bumper  Bungee  Door,  Station  389.125 

degrees  F 

Oscillograph  "B" 

4 

Oscillograph  “B” 

4 

Temperature,  Nose  Wheel  Well  Door,  Station  81.375 
Temperature,  Right  Wing,  Rear  Spar,  Lower  Cap,  Rear  Spar 

degrees  F 

Oscillograph  "B” 

199 

Oscillograph  “B” 

24 

Station  86.5 

Temperature,  Right  Wing,  Inside  Lower  Skin  Stiffener,  Rear 

degrees  F 

Oscillograph  "B" 

94* 

None 

None 

Spar  Station  138.454 

Temperature,  Right  Wing,  Inside  Lower  Skin,  Rear  Spar 

degrees  F 

Oscillograph  “B" 

26* 

None 

None 

Station  138.454 

Temperature,  Right  Outer  Wing  Panel,  Inside  Lower  Skin.  Wing 

degrees  F 

Oscillograph  “A” 

5 

Oscillograph  “A" 

5 

Station  223.75 

Temperature,  Right  Outer  Wing  Panel,  Inside  Lower  Skin,  Wing 

degrees  F 

Oscillograph  “A" 

12* 

Oscillograph  “ A ” 

12 

Station  223.75 

Temperature,  Right  Outer  Wing  Panel,  Outside  Lower  Skin, 

degrees  F 

Oscillograph  “B" 

249 

*“ 

— 

Wing  Station  222.125 

Temperature,  Right  Outer  Wing  Panel,  Inside  Lower  Skin  on 

degrees  F 

Oscillograph  “B" 

18t 

None 

None 

Blue  Door,  Wing  Station  170 

Temperature,  Right  Flap,  Inside  Lower  Skin,  Flap  Station 

degrees  F 

Oscillograph 

189 

Oscillograph  “B" 

18 

59.233  Canted 

degrees  F 

Oscillograph  “B** 

12 

Oscillograph  “B" 

12 

Temperature,  Right  Aileron,  Inside  Lower  Skin 

Temperature,  Rudder  Splitter  Plate,  Inside  Left  Skin,  Water 

degrees  F 

Oscillograph  “B" 

6 

Oscillograph  “B” 

6 

Plane  Station  69,375 

degrees  F 

Oscillograph  “B” 

5 

Oscillograph  “B” 

5 

Temperature,  Right  Elevator,  Inside  Lower  Skin,  Station  45.875 

degrees  F 

Oscillograph  “A” 

3 

Oscillograph  “A" 

3 

Temperature,  Right  Elevator,  Inside  Upper  Skin,  Station  45.875 

degrees  F 

Oscillograph  “A” 

18 

Oscillograph  “A” 

18 

Temperature,  Nose,  Radome 

degrees  F 

None 

None 

Oscillograph  “B" 

17** 

*  Applicable  on  Shot  Boltzmann  only, 
t  Applicable  on  all  shots  except  Boltzmann, 
t  Applicable  on  Shots  Boltzmann,  Hood,  and  Diablo  only. 
(  Applicable  on  Shots  Kepler,  Doppler,  and  Shasta  only. 
T  Applicable  on  Shot  Doppler  only. 

••  Applicable  on  Shot  Smoky  only. 
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TABLE  E.7  CHANNELIZATION  OF  STRUCTURAL  INSTRUMENTATION 


Instrument  or  Parameter 

Acceleration,  Lateral,  C  FI -166/34-291 
Acceleration,  Normal,  CF1-168/1-167 
Left  Aileron  Hinge  Moment 
Right  Aileron  Hinge  Moment 
Flap  Hinge  Moment 
Elevator  Hinge  Moment 
Stabilizer  Actuator  Hinge  Moment 
Rudder  Hinge  Moment 

Fuselage  Axial  Stress,  Upper  Left  Longeron,  Outboard  Flange, 
Station  250 

Fuselage  Axial  Stress,  Upper  Left  Longeron,  Inboard  Flange, 
Station  250 

Fuselage  Axial  Stress,  Upper  Left  Longeron,  Low  Web, 

Station  250 

Fuselage  Axial  Stress,  Upper  Right  Longeron,  Outboard  Flange, 
Station  250 

Fuselage  Axial  Stress,  Upper  Right  Longeron,  Inboard  Flange, 
Station  250 

Fuselage  Axial  Stress,  Upper  Right  Longeron,  Low  Web, 

Station  250 

Riselage  Torsion,  Left  Side,  Station  334.5 
Fuselage  Torsion,  Right  Side,  Station  334.5 
Wing  Bending,  Left  Root,  Wing  Station  17.5 
Wing  Bending,  Right  Root,  Wing  Station  17.5 
Wing  Bending,  Right  Root,  Wing  Station  36.5 
Wing  Bending,  Left  Rear  Spar,  Rear  Spar  Station  87.26 
Wing  Bending,  Left  Outer  Wing  Panel  Root,  Rear  Spar 
Station  173 

Wing  Bending,  Right  Outer  Wing  Panel  Root,  Rear  Spar 
Station  173 

Wing  Shear,  Left  Front  Spar,  Wing  Station  48 

Wing  Shear,  Left  Landing  Gear  Front  Spar,  Wing  Station  47.125 

Wing  Shear,  Left  Landing  Gear  Rear  Spar,  Wing  Station  39.84 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  50.987 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  63.612 

Wing  Shear,  Right  Front  Spar,  Rear  Spar  Station  50.987 

Wing  Shear,  Right  Front  Spar,  Rear  Spar  Station  63.612 

Wing  Shear,  Left  Rear  Spar,  Rear  Spar  Station  55.051 

Wing  Shear,  Lett  Rear  Spar,  Rear  Spar  Station  67.801 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  55.051 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  67.801 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  86.2 

Wing  Shear,  Right  Front  Spar,  Rear  Spar  Station  86.2 

Wing  Shear,  Left  Rear  Spar,  Rear  Spar  Station  115.125 

Wing  Shear,  Right  Rear  Spar,  Rear  Spar  Station  115.125 

Wing  Shear,  Left  Front  Spar,  Rear  Spar  Station  143.25 

Horizontal  Stabilizer  Pivot  Load,  Left 

Horizontal  Stabilizer  Pivot  Load,  Right 

Horizontal  Stabilizer  Bending,  Left,  Station  33.75 

Horizontal  Stabilizer  Bending,  Right,  Station  33.75 

Horizontal  Stabilizer  Bending,  Left,  Station  60.5 

Horizontal  Stabilizer  Torsion,  Left,  Station  30.5 

Horizontal  Stabilizer  Torsion,  Right,  Station  30.5 

Horizontal  Stabilizer  Torsion,  Left,  Station  58.5 

Fin  Bending,  Lower  Main  Beam,  Station  14.5 

Fin  Bending,  Upper  Main  Beam,  Station  89.188 

Fin  Cap  Bending,  Station  94.5 

Fin  Torsion,  Forward  Spar,  Station  46.625 

Fin  Torsion,  Forward  Spar,  Station  75.00 

Fin  Torsion,  Rear  Spar,  Station  47.125 

Fin  Torsion,  Rear  Spar,  Station  74.188 


*  Applicable  on  Shot  Boltzmann  only, 
t  Applicable  on  all  shots  except  Boltzmann, 
t  Applicable  on  Shots  Boltzmann,  Hood,  and  Diablo  only. 


FJ-4  BuNo.  139467  FJ-4  BuNo.  139310 

Recording  Recording 

Units  Instrument  Channel  Instrument  Channel 


g’s 

Oscillograph  “A” 

21 

Oscillograph  “A” 

21 

g’s 

Oscillograph  “A” 

11 

Oscillograph  “A” 

11 

lbs 

Oscillograph  “A” 

20* 

Oscillograph  “A” 

20 

lbs 

Oscillograph  “B” 

23 

Oscillograph  "B" 

23 

inch-lbs 

Oscillograph  “B” 

7* 

Oscillograph  “B" 

7 

lbs 

Oscillograph  “B” 

10 

Oscillograph  “B" 

10 

lbs 

Oscillograph  “B" 

9 

Oscillograph  “B” 

9 

inch-lbs 

Oscillograph  "B" 

11 

Oscillograph  “B" 

11 

psi 

Oscillograph  “C” 

19 

Oscillograph  “C" 

19 

psl 

Oscillograph  “D” 

13 

Oscillograph  “D” 

13 

psi 

Oscillograph  “D” 

12 

Oscillograph  “D” 

12 

psi 

Oscillograph  “D” 

10 

Oscillograph  “D” 

10 

psi 

Oscillograph  “D” 

11 

Oscillograph  “D” 

11 

psi 

Oscillograph  “C” 

21 

Oscillograph  “C” 

21 

psi 

Oscillograph  *C” 

20 

Oscillograph  “C” 

20 

psi 

Oscillograph  “D” 

5 

Oscillograph  “D” 

14 

psi 

Oscillograph  “C” 

3 

Oscillograph  “C" 

3 

psi 

Oscillograph  “D” 

4 

Oscillograph  “D” 

4 

psi 

Oscillograph  “D” 

26 

Oscillograph  “D” 

26 

psi 

Oscillograph  “C” 

4 

Oscillograph  “C" 

4 

psi 

Oscillograph  “C” 

5 

Oscillograph  “C” 

5 

psi  ^ 

Oscillograph  “D” 

6 

Oscillograph  “D” 

6 

psl 

Oscillograph  “C” 

24 

Oscillograph  “C” 

24 

psi 

Oscillograph  “C” 

7 

Oscillograph  “C” 

7 

psi 

Oscillograph  “D” 

8 

Oscillograph  “D” 

8 

psi 

Oscillograph  “C” 

6 

Oscillograph  “C” 

6 

psi 

Oscillograph  “C” 

22 

Oscillograph  “C” 

22 

psi 

Oscillograph  “D” 

7 

Oscillograph  “D” 

7 

psi 

Oscillograph  “D” 

25 

Oscillograph  “D” 

25 

psi 

Oscillograph  “C” 

9 

Oscillograph  “C” 

9 

psi 

Oscillograph  “C” 

8 

Oscillograph  “C” 

8 

psi 

Oscillograph  “D” 

22 

Oscillograph  “D” 

22 

psi 

Oscillograph  “D” 

21 

Oscillograph  “D" 

21 

psi 

Oscillograph  “C” 

10 

Oscillograph  “C” 

10 

psi 

Oscillograph  “D” 

23 

Oscillograph  “D” 

23 

psi 

Oscillograph  “C” 

25 

Oscillograph  “C” 

25 

psi 

Oscillograph  “D” 

24 

Oscillograph  “D” 

24 

psi 

Oscillograph  “C~ 

23 

Oscillograph  “C” 

23 

psi 

Oscillograph  “C” 

12 

Oscillograph  “C” 

12 

psl 

Oscillograph  “C” 

15 

Oscillograph  “D” 

15 

psi 

Oscillograph  “C” 

15 

Oscillograph  "C” 

15 

psi 

Oscillograph  “D” 

16 

Oscillograph  “D” 

16 

psi 

Oscillograph  “C” 

16 

Oscillograph  “C” 

16 

psi 

Oscillograph  “C” 

13 

Oscillograph  “C” 

13 

psi 

Oscillograph  “D” 

17 

Oscillograph  “D” 

17 

psi 

Oscillograph  “C” 

14 

Oscillograph  “C” 

14 

psi 

Oscillograph  “C” 

11 

Oscillograph  “C” 

11 

psi 

Oscillograph  "D” 

20 

Oscillograph  “D" 

20 

psl 

Oscillograph  “D” 

9 

Oscillograph  “D" 

9 

psi 

Oscillograph  “C” 

17 

Oscillograph  “C” 

17 

psi 

Oscillograph  “C” 

18 

Oscillograph  “C" 

18 

psi 

Oscillograph  “D" 

18 

Oscillograph  “D " 

18 

psi 

Oscillograph  **D” 

19 

Oscillograph  “D” 

19 

i  Applicable  on  Shots  Kepler,  Doppler,  and  Shasta  only. 

1  Applicable  on  Shot  Doppler  only. 

*•  Applicable  on  Shot  Smoky  only. 
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TABLE  E.8  CHANNELIZATION  OF  SUPPORTING  INSTRUMENTATION 


Instrument  or  Parameter 

Units 

FJ-4  BuNo.  139467 

Recording 

Instrument  Channel 

FJ-4  BuNo.  130310 

Recording 

Instrument  Chennai 

Airspeed,  CFS4-293/S4-294 

inches  of  Hg 

Oscillograph  “A” 

16 

Oscillograph  “A” 

16 

Knotmeter,  CP2-26/2-99 

knots 

Photo  Recorder 

- 

Photo  Recorder 

. 

Altitude,  CF34-S08/34-328 

inches  of  Hg 

Oscillograph  “A" 

17 

Oscillograph  "A" 

17 

AlUmeter,  CF3-1J/3-79 

feet 

Photo  Recorder 

- 

Photo  Racordar 

- 

Magnetic  Compass,  CF84-2/84-1 

degrees 

Photo  Recorder 

- 

Photo  Recorder 

• 

Angle  of  Pitch,  CF39-5S/39-8/81-S2 

degrees 

Oscillograph  "A" 

22 

Oscillograph  "A" 

22 

Angle  of  Beak,  CF39-55/39-6/81-14 

degrees 

Oscillograph  MAM 

13 

Oscillograph  "A" 

13 

Angle  of  Yaw 

degrees 

Oscillograph  “A” 

10 

Oscillograph  “A" 

10 

Bate  of  Pitch,  CF39-44/39-54 

degrees/sec 

Oscillograph  MB" 

15 

Oscillograph  “V 

15 

Rate  of  Roll,  C 939-41/39-53 

degrees/sec 

Oscillograph  “B r 

14 

Oscillograph  "B“ 

14 

Rate  of  Yaw,  CFS9-37/39-57 

degrees/sec 

Oscillograph  "8” 

10 

Oscillograph  *B” 

16 

Angle  of  Attack 

degrees 

Oscillograph  “B” 

»6t 

Oscillograph  “A" 

28 

Left  Aileron  Position 

degrees 

Oscillograph  “A” 

15* 

Oscillograph  “A" 

IS 

Right  Aileron  Poeitlon 

degrees 

Oscillograph  "B" 

13 

Oscillograph  “B" 

13 

Elevator  Position 

degrees 

Oscillograph  “A” 

13 

Oscillograph  MAT 

33 

Horisontal  Stabiliser  Position 

degrees 

Oscillograph  “A" 

7 

Oscillograph  “A” 

7 

Rudder  Position 

degrees 

Oscillograph  “A" 

14 

Oscillograph  “A” 

14 

Rudder  Pedal  Force 

lbs 

Oscillograph  "BT 

8* 

Oscillograph  “B” 

8 

Clock,  CF9-1 7/9-41 

Photo  Recorder 

- 

Photo  Recorder 

• 

Voltage  Monitor  Photocell 

volts 

Oscillograph  “A" 

24 

Oscillograph  “A” 

24 

Voltage  Monitor  Photocell 

volts 

Oscillograph  *C" 

28 

Oscillograph  "C* 

28 

Tailpipe  Temperature 

degrees  F 

Oscillograph  "A” 

1* 

Oscillograph  "A" 

1 

Tailpipe  Total  Pressure,  CF34-329/34-330 

P«U 

Oscillograph  "B” 

«• 

Oscillograph  “B” 

21 

Compressor  Inlet  Temperature 

degrees  F 

Oscillograph  "B" 

!• 

Oscillograph  “B" 

1 

Compressor  Inlet  Pressure  *1,  CF34-301/34-333 

psla 

Oscillograph  “B" 

17* 

Oscillograph  "A" 

9 

Compressor  Inlet  Pressure  *2,  CF34-302 

psia 

Oscillograph  “B" 

25* 

None 

None 

Compressor  Inlet  Pressure  #4,  CF34-304/34-336 

psla 

Oscillograph  “A" 

8* 

Oscillograph  "A” 

* 

Compressor  Inlet  Pressure  #5,  CF34-305 

psls 

Oscillograph  “A" 

«• 

Nona 

Rons 

Compressor  Discharge  Pressure,  CF34-306/34-307 

psla 

Oscillograph  "B" 

20* 

Oscillograph  “B” 

20 

Tachometer,  CF23-44/23-45 

rpm 

Photo  Recorder 

. 

Photo  Recorder 

Kiel  Flow,  CF53-S1/53-SJ 

gallons/hr 

Photo  Recorder 

Photo  Recorder 

Fuel  Temperature  at  Flow  Transmitter,  CF72 -42/72-43 

degrees  C 

Photo  Recorder 

Photo  Recorder 

* 

Fuel  Inlet  Pressure,  CF5 -34/5- 126 

peig 

Photo  Recorder 

Photo  Racordar 

Frame  Fuel  Counter,  CF12-136/12-145 

— 

Photo  Recorder 

Photo  Recorder 

Left  Fuel  Quantity,  CF53-29 

lbs 

Photo  Recorder 

Photo  Recorder 

„ 

Right  Fuel  Quantity,  CFS3-27/53-33 

lbs 

Photo  Recorder 

Photo  Recorder 

•  Applicable  on  Shot  Boltzmann  only, 
t  Applicable  on  all  shots  except  Boltzmann. 

I  Applicable  on  Shots  Boltzmann,  Hood,  and  Diablo  only. 
I  Applicable  on  Shots  Kepler,  Doppler,  and  Shasta  only, 
f  Applicable  on  Shot  Doppler  only. 

*•  Applicable  on  Shot  Smoky  only. 
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on  the  leading  edges  of  the  vertical  and  horizontal  stabilizers  provided  reference  indications 
for  observing  the  motions  of  these  structures.  These  VDR-5  cameras  were  set  for  66  frames 
per  second  for  Shot  Boltzmann  and  200  frames  per  second  for  all  other  shots. 

2.  One  VDR-5  camera  was  installed  in  the  tail  section  of  the  right  wing  tank  to  photograph 
the  fireball.  This  camera  was  set  for  32  frames  per  second  for  Shot  Boltzmann  and  200  frames 
per  second  for  all  other  shots. 

3.  One  N-9  (16  mm)  gun  camera  was  installed  in  the  cockpit  to  photograph  the  Instrument 
panel.  Flood  lighting,  controllable  from  a  switch  separate  from  the  standard  lighting  was  pro¬ 
vided.  This  camera  speed  was  16  frames  per  second. 

4.  Strain  gages,  thermocouples  and  associated  wiring  were  duplicated  wherever  possible. 
Wherever  possible,  the  strain  gages  were  temperature  compensated. 

5.  The  oscillograph  paper  speed  was  12.8  in./sec  and  the  photorecorder  speed  was  8 
frames  per  second. 

6.  A  photocell  of  proper  intensity  was  installed  in  the  fuselage  to  multiplex  with  the  galva¬ 
nometer  in  the  Voltage  Monitor  No.  1  Transducer  Power  Supply. 

7.  A  single  record  switch  was  installed  to  start  all  recording  devices,  to  close  the  over¬ 
pressure  valves,  and  to  switch  the  fuel  quantity  indication  from  the  cockpit  to  the  photore¬ 
corder. 

8.  Correlation  for  the  photorecorder  and  pilot’s  blipper  marker  was  provided  by  the  right 
reference  galvanometers  in  the  nose  compartment  and  wing  tank.  The  VDR-5  camera  in  the 
wing  tank  was  correlated  with  the  oscillographs  through  the  left  reference  galvanometer  in  the 
nose  compartment.  The  VDR-5  flexure  cameras  were  correlated  through  the  left  reference 
galvanometer  in  the  wing  tank. 

9.  Ballast  was  kept  to  a  minimum  consistent  with  the  safe  center  of  gravity  range.  All 
flights  were  made  with  the  left  wing  tank  mounted  and  usable.  A  float  switch  was  provided  in 
the  left  wing  tank  to  automatically  stop  any  fuel  transfer  and  light  an  indicator  when  the  weight 
symmetry  level  was  reached.  A  switch  mounted  in  the  left-hand  shroud  permitted  the  pilot  to 
override  this  automatic  cutoff. 

10.  The  left  wing  tank  only  was  releasable  by  the  normal  electrical  system;  however,  the 
manual  emergency  release  could  drop  both  the  left  fuel  tank  and  the  right  instrumentation  tank. 

11.  In  order  to  maintain  uniform  or  flat  frequency  response  (±2  percent)  from  zero  to  67 
percent  of  the  galvanometer  natural  frequency  wherever  possible  without  sacrificing  sensitiv¬ 
ity,  damping  networks  were  provided  to  obtain  64  percent  of  critical  damping. 

12.  All  thermocouple  wires  were  capacitance  welded  to  the  aircraft  skins  Instead  of  being 
installed  with  conventional  washers  riveted  to  the  skin.  This  method  eliminated  the  inaccura¬ 
cies  inherent  in  the  older  method  by  not  increasing  the  effective  mass  of  the  Instrumented  point 
as  well  as  substantially  reducing  thermocouple  lag. 

13.  The  frequency  response  of  the  galvanometers  employed  for  the  wing  pressure  survey 
was  135  cps. 

14.  The  frequency  response  of  the  galvanometers  employed  for  measurements  of  free 
stream  overpressure  was  60  cps. 
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Appendix  F 

THERMAL  EFFECTS  DATA  AND  CORRELATION 
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Figure  F.l  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree  calorimeters 
and  integrated  180  degree  radiometer,  Shot  Boltzmann. 


Figure  F.2  Comparison  of  the  calculated  and  measured  generalized  thermal 
pulse  taken  from  180  degree  radiometer,  Shot  Boltzmann. 
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Figure  F.4  Comparison  of  the  calculated  and  measured  temperature-time  his¬ 
tories  of  selected  thin  skins,  Shot  Boltzmann. 
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Temperature 


Figure  F.5  Temperature  versus  time  data  of  the  lower  faceplate  of  the  honey¬ 
comb  control  surfaces,  Shot  Boltzmann. 
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Jo  'aminiaduiej. 


Figure  F.6  Comparison  of  the  calculated  and  measured  temperature-time  his¬ 
tories  of  selected  honeycomb  faceplates,  Shot  Boltzmann. 
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Figure  F.7  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree  calorime¬ 
ters  and  integrated  180  degree  radiometer,  Shot  Diablo. 
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Figure  F.9  Temperature  versus  time  data  of  the  thin  skins  on  the  fuselage  and 
right  outer  wing  panel,  Shot  Diablo. 
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Tamparatura,  °F 


Tima,  seconds 

Figure  F.10  Comparison  of  the  calculated  and  measured  temperature -time 
histories  of  selected  thin  skins,  Shot  Diablo. 
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110 


Figure  F.12  Comparison  of  the  calculated  and  measured 
temperature -time  histories  of  selected  honeycomb 
faceplates,  Shot  Diablo. 
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Figure  F.13  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree  calorime¬ 
ters,  Shot  Kepler. 


Figure  F.14  Temperature  versus  time  data  for  the  thin  skins  on  the  fuselage  and 
right  outer  wing  panel,  Shot  Kepler. 
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Figure  F.15  Temperature  versus  time  data  for  the  lower  faceplates  of  the 
honeycomb  control  surfaces.  Shot  Kepler. 
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Figure  F.  16  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree  calorime¬ 
ters  and  integrated  180  degree  radiometer,  Shot  Shasta. 
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Figure  F.17  Comparison  of  the  calculated  and  measured  generalized  thermal 
pulse  taken  from  180  degree  radiometer.  Shot  Shasta. 
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Figure  F.18  Temperature  versus  time  data  of  the  thin  skins  on  the  fuselage 
and  right  outer  wing  panel,  Shot  Shasta. 
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Figure  F.  19  Comparison  of  the  calculated  and  measured  temperature-time 
histories  of  selected  thin  skins,  Shot  Shasta. 
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Figure  F.20  Temperature  versus  time  data  of  the  lower  faceplates  e 
honeycomb  control  surfaces,  Shot  Shasta. 
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Figure  F.21  Comparison  of  the  calculated  and  measured 
temperature-time  histories  of  selected  honeycomb  face¬ 
plates,  Shot  Shasta. 
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Figure  F.22  Measured  radiant  exposure  versus  time  taken  from  90  and  180  degree 
calorimeters  and  Integrated  180  degree  radiometer,  Shot  Doppler. 
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Figure  F.23  Comparison  of  the  calculated  and  measured  generalized  thermal 
pulse  taken  from  180  degree  radiometer,  Shot  Doppler. 
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Figure  F.24  Temperature  versus  time  data  of  the  thin  skins  on  the  fuse¬ 
lage  and  right  outer  wing  panel,  Shot  Doppler. 
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Figure  F.25  Comparison  of  the  calculated  and  measured  temperature-time 
histories  of  selected  thin  skins,  Shot  Doppler. 
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Figure  F.26  Temperature  versus  time  data  of  the  lower  faceplates  of  the 
honeycomb  control  surfaces,  Shot  Doppler. 
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Figure  F.27  Comparison  of  the  calculated  and  measured  temperature-time 
histories  of  selected  honeycomb  faceplates,  Shot  Doppler. 
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Figure  G.l  Comparison  of  calculated  and  measured  gust  induced  perturba¬ 
tion  of  normal  load  factor  versus  time.  Shot  Boltzmann. 
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Figure  G.2  Comparison  of  calculated  and  measured  gust  induced  perturba¬ 
tion  of  pitch  angle  versus  time,  Shot  Boltzmann. 
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figure  G.3  Comparison  of  calculated  and  measured  pitch  rate  versus  time, 
Shot  Boltzmann. 
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Figure  G.4  Measured  gust  Induced  perturbation  of  horizontal  stabilizer 
position  versus  time,  Shot  Boltzmann. 


141 

CONFIDENTIAL 


1 


figure  G.5  Comparison  of  calculated  and  measured  gust  induced  incre¬ 
mental  wing  bending  moment  stress  level  at  right  Wing  Station  17.5  versus 
time.  Shot  Boltzmann. 


Figure  G.6  Comparison  of  calculated  and  measured  gust  induced  incre¬ 
mental  wing  bending  moment  stress  level  at  left  Wing  Station  17.5  versus 
time.  Shot  Boltzmann. 
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Figure  G.7  Measured  gust  induced  incremental  wing  bending  moment  stress 
level  at  Rear  Spar  Station  173  versus  time,  Shot  Boltzmann. 


Figure  G.8  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  50.987  versus  time,  Shot  Boltzmann. 
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Figure  G.9  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  63.612  versus  time,  Shot  Boltzmann. 


Figure  G.10  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  115.125  versus  time.  Shot  Boltzmann. 
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Figure  G.ll  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  143.25  versus  time.  Shot  Boltzmann. 
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Figure  G.12  Comparison  of  calculated  and  measured  gust  induced  perturbs* 
tion  of  normal  load  factor  versus  time,  Shot  Hood. 
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Figure  G.13  Comparison  of  calculated  and  measured  gust  induced  perturba¬ 
tion  of  pitch  angle  versus  time.  Shot  Hood. 


Figure  G.14  Comparison  of  calculated  and  measured  pitch  rate  versus  time, 
Shot  Hood. 
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Incremental  Horizontal  Stabilizer  Position,  degrees  Incremental  Angle  of  Attack,  degree* 
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Figure  G.15  Measured  gust  induced  perturbation  of  angle  of  attack  versus 
time.  Shot  Hood. 
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Figure  G.16  Measured  gust  induced  perturbation  of  horizontal  stabiliser 
position  versus  time,  Shot  Hood. 
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Figure  G.17  Comparison  of  calculated  and  measured  gust  induced  incre¬ 
mental  ving  bending  moment  stress  lerel  at  right  Wing  Station  17.5  versus 
time.  Shot  Hood. 


Figure  G.18  Comparison  of  calculated  and  measured  gust  Induced  incre¬ 
mental  ving  bending  moment  stress  level  at  left  Wing  Station  17.5  versus 
time.  Shot  Hood. 
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Figure  G.19  Measured  gust  induced  Incremental  wing  bending  moment 
stress  level  at  Rear  Spar  Station  173  versus  time.  Shot  Hood. 


Figure  G.20  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  50.987  versus  time,  Shot  Hood. 
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Figure  G.21  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  63.612  versus  time.  Shot  Hood. 
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Figure  G.22  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  115.125  versus  time,  Shot  Hood. 
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Figure  G.23  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Bear  Spar  Station  143.25  versus  time.  Shot  Hood. 


Figure  G.24  Comparison  of  calculated  and  measured  gust  induced  perturba¬ 
tion  of  normal  load  factor  versus  time,  Shot  Diablo. 
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Figure  G.25  Comparison  of  calculated  and  measured  gust  induced  perturba¬ 
tion  of  pitch  angle  versus  time.  Shot  Diablo. 
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Figure  G.26  Comparison  ot  calculated  and  measured  pitch  rate  versus  time, 
Shot  Diablo. 
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Figure  G.27  Measured  gust  induced  perturbation  of  angle  of  attack  versus 
time.  Shot  Diablo. 
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Figure  G.28  Measured  gust  Induced  perturbation  of  horizontal  stabilizer 
position  versus  time,  Shot  Diablo. 
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Figure  G.29  Comparison  of  calculated  and  measured  gust  induced  incre¬ 
mental  wing  bending  moment  stress  level  at  right  Wing  Station  17.5  versus 
time,  9iot  Diablo. 


Figure  G.30  Comparison  of  calculated  and  measured  gust  induced  incre¬ 
mental  wing  bending  moment  stress  level  at  left  Wing  Station  17.5  versus 
time,  Shot  Diablo. 
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figure  G.S1  Measured  gust  Induced  incremental  wing  bending  moment 
stress  level  at  Rear  Spar  Station  173  versus  time.  Shot  Diablo. 


Figure  G.32  Measured  gust  Induced  Incremental  wing  shear  stress  level  at 
Rear  Spar  Station  50.987  versus  time,  Shot  Diablo. 
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Figure  G.3S  Measured  gust  Induced  Incremental  wing  shear  stress  level  at 
Rear  Spar  Station  63.613  versus  time.  Shot  Diablo. 
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Figure  G.34  Measured  gust  induced  Incremental  wing  shear  stress  level  at 
Rear  Spar  Station  115.125  versus  time,  Shot  Diablo. 
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Figure  G.35  Measured  gust  induced  incremental  wing  shear  stress  level  at 
Rear  Spar  Station  143.25  versus  time.  Shot  Diablo. 


Figure  G.36  Comparison  of  calculated  and  measured  gust  induced  perturbation 
of  normal  load  factor  versus  time,  Shot  Kepler. 
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Figure  G.37  Comparison  of  calculated  and  measured  gust  induced  perturbation 
of  pitch  angle  versus  time.  Shot  Kepler. 


Figure  G.38  Comparison  of  calculated  and  measured  pitch  rate  versus 
time.  Shot  Kepler. 
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Figure  G.39  Measured  gust  induced  perturbation  of  angle  of  attack  versus 
time.  Shot  Kepler. 
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Figure  G.40  Measured  gust  Induced  perturbation  of  horizontal  stabilizer 
position  versus  time.  Shot  Kepler. 
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Figure  G.41  Comparison  of  calculated  and  measured  gust  induced  incre¬ 
mental  wing  bending  moment  stress  level  at  right  Wing  Station  17.5  versus 
time,  Shot  Kepler. 


Figure  G.42  Comparison  of  calculated  and  measured  gust  Induced  incre¬ 
mental  wing  bending  moment  stress  level  at  left  Wing  Station  17.5  versus 
time,  Shot  Kepler. 
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Figure  G.44  Measured  gust  induced  Incremental  wing  shear  stress  level 
at  Rear  Spar  Station  50.987  versus  time,  Shot  Kepler. 
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Figure  G.4S  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  63.612  versus  time.  Shot  Kepler. 
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Figure  G.46  Measured  gust  Induced  incremental  wing  shear  stress  level 
at  Rear  Spar  station  115.125  versus  time,  Shot  Kepler. 
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Inrr«m#ntal  Normal  Load  Facto* ,  gravltUt  Incramantof  5»ftw,  ptl 


Igure  G.4?  Measured  gust  induced  incremental  wing  shear  stress  level 
Lt  Rear  Spar  Station  143.25  versus  time.  Shot  Kepler. 


Figure  C.48  Comparison  of  calculated  and  measured  gust  Induced  pertur¬ 
bation  of  normal  load  factor  versus  time,  Shot  Doppler. 
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Figure  G.49  Comparison  of  calculated  and  measured  gust  induced  pertur¬ 
bation  of  pitch  angle  versus  time.  Shot  Doppler. 


Figure  G.50  Comparison  of  calculated  and  measured  pitch  rate  versus 
time,  Shot  Doppler. 
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Figure  G.S1  Measured  gust  Induced  perturbation  of  angle  of  attack  versus 
time.  Shot  Doppler. 
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Figure  G.52  Measured  gust  Induced  perturbation  of  horizontal  stabilizer 
position  versus  time.  Shot  Doppler. 
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Figure  G.53  Comparison  of  calculated  and  measured  gust  Induced  incremental 
ving  bending  moment  stress  level  at  right  Wing  Station  17.5  versus  time.  Shot 
Doppler. 
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Figure  G.54  Comparison  of  calculated  and  measured  gust  induced  incremental 
wing  bending  moment  stress  level  at  left  Wing  Station  17.5  versus  time,  Shot 
Doppler. 
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Figure  G.55  Measured  gust  induced  incremental  wing  bending  moment 
stress  level  at  Rear  Spar  Station  173  versus  time.  Shot  Doppler. 


figure  G.56  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  50.987  versus  time,  Shot  Doppler. 
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Figure  G.57  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  63.612  versus  time.  Shot  Doppler. 
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Figure  G.58  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  115.125  versus  time,  Shot  Doppler. 
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Figure  G.S9  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Bear  Spar  Station  143.25  versus  time,  Shot  Doppler. 
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Figure  G.6C  Comparison  of  calculated  and  measured  gust  induced  pertur¬ 
bation  of  normal  load  factor  versus  time.  Shot  Smoky. 
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Figure  G.61  Comparison  of  calculated  and  measured  gust  induced  pertur¬ 
bation  of  pitch  angle  versus  time.  Shot  Smoky. 
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Figure  G.62  Comparison  of  calculated  and  measured  pitch  rate  versus 
time.  Shot  Smoky. 
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Figure  G.63  Measured  gust  induced  perturbation  of  horizontal  stabilizer 
position  versus  time,  Shot  Smoky. 


Figure  G.64  Comparison  of  calculated  and  measured  gust  induced  incremental 
wing  bending  moment  stress  level  at  left  Wing  Station  17.5  versus  time,  Shot 
Smoky. 
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Figure  G.65  Measured  gust  induced  incremental  wing  bending  moment 
stress  level  at  Rear  Spar  Station  173  versus  time,  Shot  Smoky. 
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Figure  G.66  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  50.987  versus  time,  Shot  Smoky. 
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Figure  G.67  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  63.612  versus  time.  Shot  Smoky. 
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Figure  G.68  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  115.125  versus  time,  Shot  Smoky. 

173 

CONFIDENTIAL 


0  0.2  0.4  0.6  0.8  1.0  1.2 


Time  After  Shock  Arrival,  seconds 

Figure  G.69  Measured  gust  induced  incremental  wing  shear  stress  level 
at  Rear  Spar  Station  143.25  versus  time,  Shot  Smoky. 
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Figure  H.l  Measured  section  lift  at  Wing  Station  175.75  versus  time.  Shot  Hood. 
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Figure  H.2  Measured  section  center  of  pressure  at  Wing  Station  175.75  versus 
time,  Shot  Hood. 
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Figure  H.3  Measured  section  lift  at  Wing  Station  175.75  versus  time.  Shot  Diablo. 
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Figure  H.5  Measured  section  lift  at  Wing  Station  175.75  versus  time.  Shot  Kepler. 


Figure  H.6  Measured  section  center  of  pressure  at  Wing  Station  175.75  versus 
time,  Shot  Kepler. 
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Figure  H.7  Measured  section  lift  at  Wing  Station  175.75  versus  time,  Shot  Doppler. 


Figure  H.8  Measured  section  center  of  pressure  at  Wing  Station  175.75  versus 
time,  Shot  Doppler. 
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Figure  1.1  Calibrated  airspeed,  engine  speed,  and  fuel  flow  time  histories,  Shot  Hood. 
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Figure  1.2  Calibrated  airspeed,  engine  speed,  and  fuel  flow  time  histories,  Shot  Diablo. 
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Figure  L3  Calibrated  airspeed,  engine  speed,  and  fuel  flow  time  histories,  Shot  Kepler. 
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Figure  L4  Calibrated  airspeed,  engine  speed,  and  fuel  flow  time  histories,  Shot  Shasta. 
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Figure  1.5  Calibrated  airspeed,  engine  speed,  and  fuel  flow  time  histories,  Shot  Doppler. 
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Figure  1.6  Tailpipe  total  pressure,  tailpipe  total  temperature,  and  compressor  inlet  total 
pressure  time  histories,  Shot  Doppler. 
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Figure  1.8  Tailpipe  total  pressure,  tailpipe  total  temperature,  and  compressor  inlet  total 
pressure  time  histories,  Shot  Smoky. 
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